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ABSTRACT 


An electrical signal obtained during the plastic deformation of sodium 
chloride crystals is described. This effect is attributed to the movement of 
charged dislocations. During a first series of experiments the specimens were 
submitted to a cyclic deformation. The magnitude of the electrical effect 
was measured in function of the amplitude of deformation. The influence of 
annealing was also investigated. In a second series of experiments the 
erystals were deformed by a series of stress pulses. A tentative explanation 
for the shape of the obtained signal is given. i 


$1. INTRODUCTION 


By extending Lehovec’s theory (Lehovec 1953) of a space-charge layer at 
the surface of an ionic crystal to a dislocation, Eshelby ei al. (1958) con- 
cluded that dislocations in an ionic crystal should be electrostatically 
charged. To test this assumption, they suggested that plastic deforma- 
tion of an ionic crystal would probably give some kind of an electrical 
_ effect. 

In two previous notes (Amelinckx ef al. 1959, Remaut et al. 1960) it was 
shown that an electrical effect is produced when a NaCl crystal is 
submitted to a cyclic stress. We then suggested that this effect was 
probably due to the motion of charged dislocations during the plastic 
deformation. In order to verify this hypothesis, a series of experiments 
was performed of which some results have already been described in the 


previous notes. We are now able to give some more results which 
support this hypothesis. 


+ Communicated by Professor Dr. W. Dekeyser. 
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§ 2. EXPERIMENTAL 
2.1. Deformation by Cyclic Bending 


In a first series of experiments practically the same apparatus was used 
as described before. The NaCl specimens were rectangular slabs of 
about 7mm x3mmx45 mm cut from crystals grown in our laboratory. 
The specimens were then bent to introduce an excess of dislocations of the 
same mechanical sign. By varying the radius of curvature we could 
change the dislocation density. One end of the bent specimen was 
rigidly mounted, the other was fastened to the moving part of the vibrator 
(fig. 1). The amplitude of vibration was measured by a vibrating capacitor 
mounted on the crystal. The calibration of the capacitor was performed 
by measuring the amplitude directly with a microscope. 


Bio. 


to electrometer 
tube 


Bent specimen AB with end A rigidly mounted and end B fastened to the 
vibrator. 


The electrical effect was detected by evaporating a gold electrode on 
the crystal and connecting it by means of a gold wire (0-005 mm diameter) 
to an electrometer with a very high input impedance and further to one 
of the channels of a dual-trace oscilloscope. By connecting the vibrating 
capacitor to the second channel of the oscilloscope it was possible to. 
compare the signal coming from the specimen with its deformation. 

A photo of the electrical effect obtained during cyclic bending is given in 
fig. 24. The shape of this signal, except for the relative length of the 
horizontal step, stayed practically the same in the whole frequency range 
considered (10-3000 c/s). In fig. 3 the signal coming from the NaCl 
crystal is compared with that of the capacitor. Here, one can clearly see 
the shift of the two curves with respect to a common time axis. This 
shift is independent of the frequency in the range considered. 


+ Figures 2, 3, 7 and 10 are shown as plates. 
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Fig. 4 
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Relationship between the amplitude of the electrical signal and the amplitude 
of vibration. 


Fig. 5 


to electrometer 
tube 


[eS Se hes ates 


Bent specimen AB rigidly mounted at the end B and fastened to the vibrator 
at the end A. 
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In fig. 4, curve A gives the relationship for a particular specimen 
between the amplitude of the electrical signal and that of the vibration. 
Specimens with the same dimensions and curvature, cut from the same 
mother crystal, gave practically identical curves. For specimens of 
different origin, the value of amplitude of the signal, for a same mechanical 
amplitude, depends very much on the history of the specimen. 

Curve B of fig. 4 gives an analogous relationship to curve A but now for a 
specimen that was first bent and then annealed at 600°c for 48 hours. 
The electrical effect during cyclic stressing is now much smaller. A possible 
explanation for this is the pinning of the dislocations during annealing. 


Fig. 6 


(a) (b) 


Movement of dislocations during an extension and compression cycle: (a) 
represents the movement during extension of the specimen; (b) corres- 
ponds to compression. 


2.2. Specimen Submitted to an Extension and Compression Cycle 


In a second series of experiments the procedure for stressing the 
specimens was changed. The specimens which had been given a small 
curvature to have an excess of dislocations of the same mechanical 
sign, were now mounted in order to subject them to an extension and 
compression cycle (fig. 5). 

The electrical signal obtained this way had exactly the same shape as 
that obtained during the first series of experiments. The displacement 
between the electrical and mechanical signal was also the same. But 
contrary to what we found during cyclic bending, the signals coming from 
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opposite faces of the specimen were now in opposition. This is to. be 
expected, because when the specimen is compressed or extended, dis- 
locations will always leave one face of the specimen while at the same tiine 
others will arrive at the opposite face. This is visualized in fig. 6. 

With this experimental procedure it was also possible to determine the 
shear stress at which the electrical signal begins, i.e. the shear stress at 
which the dislocations begin to move in great numbers. The measured 
value for different specimens lay in the vicinity of 75 g/mm? which is 
very close to the value given by Schmid and Boas (1950) for the yield stress 
of NaCl. We can thus conclude that the electrical effect begins at the 
moment we get plastic deformation. 


Fig. 8 
G, 


A possible explanation for the shape of the signal corresponding to alternating 
pulses. 


2.3. Pulse Stressing 


In order to investigate the effect of a stress increasing linearly, with 
time, specimens were deformed by feeding a series of sawtooth current 
pulses to the vibrator on which the crystal was mounted. The circuit: 
was arranged in such a way that successive pulses could have the same or 
the opposite direction. In the latter case the specimen was alternately 
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extended and compressed while in the former the specimen was successively 
compressed. This type of experiments gave some interesting results. 

The electrical signal corresponding to alternated pulses is seen in fig. 7. 
A possible explanation for the shape of the signal obtained is the following. 
During the period (t, —¢,) [fig. 8 (a)], the stress rises linearly from zero to a 
maximum value. It then drops instantaneously to zero and stays so 
until ¢,, where it starts increasing linearly again, but now in the reversed 
direction. This continues till ¢;, where it drops back to zero again. The 
specimens hysteresis loop which should correspond to this cycle of 
deformation is shown in [fig. 8 (b)]. From this, one can deduce the form 
of the electrical signal [fig. 8 (c)]._ At ¢) the specimen starts to deform but 
the signal stays zero till a critical stress is reached where we then obtain a 
rapid increase of the signal. At ¢,, the stress drops abruptly to zero, and 
we get an elastic recoil of the specimen which causes a slight decrease of 
the observed signal. During the period (t,—t,) there is no movement of 
dislocations and the electrical signal does not vary. At ¢, the saw-tooth 
pulse starts in the opposite direction and the cycle repeats itself. 


Fig. 9 
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A schematic picture showing how the electrical signal changes during continued 
pulsing. 


During sawtooth cycling the signal changes slightly in amplitude and 
in form. In fig. 9 the full line represents the signal corresponding to the 
first pulses; the dotted curve represents the signal after equilibrium has 
been reached. As can be seen the signal decreases. We suppose that 
this results from the damage done in the glide-plane during continued 
deformation; probably the rate of deformation will also have a big 
influence. This should be further investigated. 

The result of imposing on the specimen a series of pulses of the same 
polarity, can be seen in fig. 10. As was to be expected during the first 
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pulse (fo, ¢,) [fig. 11 (c)] the signal changed from zero to a maximum value, 
At ¢, where the stress dropped to zero, the electrical effect decreased 
slightly as a result of the elastic recoil, but it then stayed constant till t,. 
At ¢, a new pulse started and a slight increase in the signal was obtained 
owing to the small elastic and plastic deformation which accompanied. 
this pulse. This continued for some time till the signal reached a certain 


Fig. 11 


(a) 


ae imen, the elongation (e), 
hip between the stress (c) applied to the specimen, 
a the electrical signal (#) and the time (t). 


maximum value where it then stayed constant probably as the result of 
the piling up of dislocations. By slightly raising the stress the Siac te 
were made to break through these barriers and a sudden ee in signa. 
was observed. If this was repeated several times one remarked that the 


increase in signal which accompanied the increase in stress, diminished. 


8 On an Electrical Effect obtained during Deformation 


One finally arrived at the situation where an increase in stress had praeti- 
cally no influence any more on the electrical effect. At this moment a 
slight increase in stress usually resulted in the fracture of the specimen. 


§ 3. CONCLUSION 


From these experiments one can conclude that the electrical effect 
studied is closely related to the movement of dislocations. As a result of 
this, this method of investigation should be a promising tool for the study 
of such phenomenon as plastic flow, work-hardening, Peierls-Nabarro 
forces and fatigue of ionic crystals. 


ACKNOWLEDGMENT 


This work is part of a research scheme supported by the Institut pour 
Encouragement de la Recherche Scientifique dans I|’Industrie et 
V Agriculture (C.E.S.). We wish to thank Professor Dr. W. Dekeyser for 
his continuous interest. 


REFERENCES 


AMELINCKX, 8., VENNIK, J., and Kemavt, G., 1959, J. Phys. Chem. Solids, 11, 
170. 

EsHetsy, J. D., Newry, C. W. A., Prarr, P. L., and Liprarp, A. B., 1958, 
Phil. Mag., 3, 76. 

Lrnovec, K., 1953, J. chem. Phys., 21, 1123. 

Remavt, G., Vennik, J., and AMELINCKX, S., 1960, J. Phys. Chem. Solids (to 
be published). 

See and Boas, W., 1950, Plasticity of Crystals (London: F. A. Hughes & 

o. Ltd.). 


Effect of Slip Distribution on the Fracture Behaviour of 
Magnesium Oxide Single Crystals} 


By R. J. Stoxss, T. L. JOHNSTON, and C. H. Li 
Honeywell Research Center, Hopkins, Minnesota, U.S.A. 


[Received June 23, 1960] 


ABSTRACT 


The tensile deformation of chemically polished magnesium oxide erystals 
has been correlated with the distribution of slip at the onset of plastic flow. If 
two slip bands, generated on orthogonal {110} planes, happen to intersect 
to nucleate a crack before other slip bands develop nearby, then the crack so. 
formed is unstable and the crystal completely brittle. If, however, there 
is a higher density of slip sources cracks can become stabilized by adjacent 
slip bands and the crystals continue to deform. Sometimes slip is confined to. 
a single slip band which expands laterally to fill the entire gauge length, these 
crystals are extremely ductile. The fracture behaviour of a given crystal 
depends critically upon the relative orientation, number, thickness and 
spacing between slip bands. 

The density of slip may be increased artificially when disiocation sources: 
are injected into the crystal surface by sprinkling with carborundum before: 
loading. Such crystals are always ductile in tension. A similar treatment 
for crystals to be bent leads to a profound change in their fracture behaviour:. 


§ 1. INTRODUCTION 


WHEN the mechanical properties of different magnesium oxide crystals 
are compared with one another there is an alarming lack of reproducibility 
both in their yield strength and ductility. Studies of variation in the yield 
strength have led to the general conclusion that a torm of precipitation 
hardening is involved (Kronberg et al. 1960, Thompson and Roberts 
1960, Gorum et al. 1960) so that the rate of cooling from an annealing 
temperature (1000°cC and above) is particularly critical. Indeed it is 
possible to halve the yield strength, by quenching rather than slowly cooling 
the crystals. In the present work we have equalized the yield strength of 
a number of crystals with a standard annealing treatment and then studied 
some of the factors controlling their ductility. 

In some cases the variation in ductility may be attributed to the prepara- 
tion of the specimen surface since cleavage and subsequent handling are 
now known to introduce surface defects which can propagate to cause 
premature failure under stress (Stokes et al. 1960, Clarke and Sambell 
1960). However, when the precaution is taken to eliminate defects by 
chemically polishing and carefully handling the crystals the scatter in 
ductility still persists. It will be shown below that this innate scatter 


+ Communicated by the Authors. 
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may be interpreted in terms of the dependence of the fracture behaviour 
upon certain slip parameters. 

An earlier paper (Stokes et al. 1959a) has described the relationship 
between slip and the cleavage fracture of polished magnesium oxide 
crystals deformed at a slow strain rate under bending. For the sake of 
continuity and the need for reference the details of this fracture mechanism 
will be briefly reviewed at this point. It was concluded from the previous 
observations that a cleavage crack nucleated at the end of a (110) slip 
plane by the coalescence of edge dislocations which became held up against 
an obstacle provided by dislocations concentrated in an orthogonal 
intersecting (110) slip band. The crack nucleus lay parallel to the (110) 
plane of the orthogonal slip band and propagated by cleavage over this 
plane until halted by another (110) slip band parallel and adjacent to 
the original (110) slip band. Although the crack was prevented from 
propagating further in the [110] direction by the slip band, it remained free 
to develop parallel to the [001] line of intersection. In this way a stable 
internal (110) [001] slit was formed along the channel between the two slip 
bands. For complete fracture of the specimen the slit switched out of the 
(110) plane into the (100) plane normal to the applied stress. All crystals 
fracturing by this mechanism contained a characteristic (110) [001] step 
across their fracture surface, referred to descriptively as the ‘line’ source. 

There are two factors which are most important in this fracture 
mechanism. First, slip must be activated on two orthogonal planes 
which intersect each other before the crack can nucleate. Second, the 
effectiveness of the adjacent band in halting and then stabilizing the 
growing crack nucleus depends on its thickness and its distance from the 
[001] line of nucleation. The distance from the line of nucleation (or slip 
band separation) is particularly critical since it determines (a) the velocity 
attained by the accelerating crack front before it hits the slip band and (bd) 
the stabilized crack length, which in turn determines its associated stress 
concentration. An additional factor which could be important is that the 
slip bands must be sufficiently thick to act as an obstacle to an edge 
dislocation avalanche for crack nucleation. 

With this mechanism it is possible to appreciate how different. slip 
band configurations can modify the fracture process and thereby the 
ductility of magnesium oxide crystals. For example, there is the situation 
where a crystal contains just two orthogonal slip bands which happen to 
intersect each other to nucleate a crack before any other slip bands develop. 
Under these conditions crack propagation is uninhibited and the crystal 
extremely brittle. However, if there is a higher density of intersecting 
slip bands the cracks that form become stabilized and a certain amount of 
ductility is to be expected. In addition, there is the configuration where 
slip is confined to a single plane only, in which case no intersections occur, 
no pile-ups form and the conditions for crack nucleation by the above 
mechanism do not exist within the crystal. Such crystals would be 
expected to be extremely ductile. The fracture behaviour of a given 
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‘crystal, therefore, depends critically upon the relative orientation, number, 
thickness and spacing between slip bands operative from the onset of 
plastic flow. 

In § 3 of this paper it will be shown that the innate scatter in the ductility 
of chemically polished magnesium oxide crystals when tested in tension 
may be correlated with their slip behaviour in the above manner. In 
§§ 4 and 5 experiments are described in which the slip distribution at the 
onset of plastic flow was controlled artificially, a procedure which made 
it possible to increase the ductility reproducibly and in some instances to 
modify the fracture behaviour quite profoundly. 


§ 2. EXPERIMENTAL PROCEDURE 


As emphasized elsewhere (Stokes et al. 1960) a detailed study of the 
mechanical properties of rock-salt type ionic solids is possible only if 
exceptional care is taken in the surface preparation and handling of the 
crystals. In all of the work reported here magnesium oxide crystals were 
initially chemically polished and then untouched within their gauge 
length (apart from the special procedures of §§4 and 5) so that the 
observations refer to their innate mechanical behaviour. 

The specimens to be tested in tension were first cleaved from bulk — 
crystals (supplied by Norton Co.) to dimensions approximately 0-1 in. 
square and # in. in length. They were next given a standardized annealing 
treatment in which the temperature was raised to 2000°c in 15min, 
held there for 20min and then slowly cooled to room temperature over 
a period of 30min. Asa consequence all of the crystals yielded at approxi- 
mately the same resolved shear stress, 5000 p.s.i. After polishing in hot 
orthophosphorie acid they were carefully mounted in split tensile grips 
with an epoxy resin cement to leave an approximate 3in. gauge length. 
The tests were conducted in a hard tensile machine at a low strain rate of 
1:5x10->sec. The alignment of each crystal was checked visually in 
the early stage of the tension test with a polarized light technique. Any 
non-uniform stress distribution was revealed by a colour gradient across the 
‘specimen in which case the specimen was unloaded and re-aligned. In 
general the stress distribution throughout the specimen was fairly uniform 
although constraints in the vicinity of the grips (the specimens were not 
machined or polished to give a reduced gauge section) gave rise to local 
‘stress concentrations. ; 

Some of the experimental observations described in this paper were 
‘recorded on motion-picture film. The crystals used for these experiments 
were illuminated with polarized light and the pictures taken through an 
-analyser crossed with the polarizer for extinction under zero stress. 
Because of the high stresses involved in plastic flow and the thickness of 
the crystal the birefringent pattern produced was highly coloured. To 
record this colour 16mm Kodachrome film was used. The photographs 
reproduced in this paper are direct, and therefore negative, enlargements 


-from the film. 
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To change the distribution of slip bands, artificial slip sources were 
introduced into the specimen in two days, first by cleavage and second 
by sprinkling with silicon carbide powder. Dislocations produced along: 
cleavage tear lines may be expanded under an applied tensile stress (Stokes 
et al. 1959b, Gilman 1959). Accordingly, crystals were prepared: whose 
dimensions were such that successive division into two equal parts yielded 
a specimen of the appropriate, 0-lin. x 0-lin. cross section. The large 
crystals were first chemically polished. Special precautions were then. 
taken to avoid the introduction of undesirable surface microcracks during 
the act of cleavage. To do this the crystals were cleaved while standing 
upright on a piece of filter paper and were given a sharp tap with the 
chisel to propagate a fast crack. In this way specimens were prepared. 
which had either one surface or two opposite surfaces in the as-cleaved 
condition. 

A more uniform distribution of slip sources was obtained by sprinkling 
chemically polished and mounted tensile specimens liberally on either one, 
two opposite, or all four {100} surfaces with fine particles of 200 mesh 
silicon carbide powder from a height of 3in. The impact produced by 
these particles was just sufficient to inject dislocation half loops into the 
surface of the crystals, yet insufficient to form microcracks (Stokes et al. 
1959b). In a few instances, however, sprinkled specimens did fracture. 
prematurely due to a surface defect introduced by sprinkling; such 
tests were omitted from the present discussion. 

Crystal beams to be tested in bending were cleaved to a } in. x} in. 
cross section and 1 in. length, chemically polished and then when desired. 
thoroughly sprinkled on all four surfaces with 200 mesh silicon carbide. 
After being loaded in air to the stress level at which yielding could just. 
be detected, the crystal was removed from the machine and repolished. 
The purpose of the second polish was to remove any surface flaws which. 
may have been introduced in sprinkling, a procedure which obviously had. 
to wait until the dislocation loops had been expanded deep enough into. 
the crystal or they too would have been eliminated. (A similar repolishing: 
procedure was not possible with the tensile specimens because they were. 
mounted in grips.) The deformation was then continued to fracture. 

A number of specimens which had been sprinkled and pre-strained. 
slightly by bending to introduce dense slip were further deformed by- 
impact. A small Charpy-type impact machine was used with a capacity 
of 10in./lb and a hammer velocity such that the outer fibre was strained. 
at a rate of 40sec—!, an increase of 108 over the low strain rate tests. The. 
direction of impact was in the same direction as the pre-strain. 


§ 3. [ytRinstc Sire SourcES 
The plastic behaviour of polished magnesium oxide crystals under: 


tension may be broadly classified into three categories, each of which will 
be shown to be characterized by a certain distribution of slip. As can be. 
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seen from fig. 1, the deformation before fracture increased from category 
I to IT to II, even though the yield stress remained consistently in the 
range 4000 to 5000p.s.i. For category I the crystals deformed less than 
0-2°% strain; in some instances plastic deformation was hardly detectable 
and macroscopically the specimens appeared completely brittle. All 
crystals of category I fractured within the gauge length (as indicated by 
the vertical arrow). Crystals of category II were more ductile and fractured 
within the gauge length after approximately 5% strain. For category III 
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Tensile stress-strain curves for chemically polished magnesium oxide crystals 
tested at room temperature. 


there was a distinct change in the yielding behaviour in that the stress-strain 

curves were initially flat, exhibiting in somes cases a yield point. Such 

crystals were extremely ductile and never fractured within the gauge 

length (indicated by the horizontal arrow), they either pulled loose from 
e grips. 

¥ peel abe ene turn the details of the slip band distribution and 

fracture for each of these three situations. 
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3.1. Category I 


Figure 2} shows the slip band distribution, as revealed by etch pits,. 
for the specimen represented by curve (a) in fig. 1. It can be seen that 
along the whole gauge length there were no more than ten individual 
bands of slip, and, most important, two of these bands intersected along a 
line coincident with the fracture surface. The appearance of the fracture 
surface for this particular specimen is reproduced in fig. 3 and shows the 
characteristic (110) [001] line source. 

Figure 4 is a composite photomicrograph taken on the two fractured 
halves of the specimen represented by curve (6) in fig. 1. The individual 
pieces were cleaved down their middle after fracture to expose surfaces. 
suitable for etching, hence the cleavage lines in the photographs. On the 
basis of the stress-strain curve this crystal appeared completely brittle 
but after etching was found to contain just two slip bands within the 
gauge length, the two illustrated in the figure. Again, cleavage fracture 
originated from the (110) [001] slit formed along their line of intersection. 

In summary, it was concluded that crystals of category I were brittle 
because cleavage cracks nucleated at slip band intersections before other 
slip bands had developed which could have stabilized them. 


3.2. Category II 


Crystals of this category had a slip distribution essentially the same as. 
category I except that there was now a sufficient density. of slip sources 
operating from the onset of flow to stabilize the cracks as they formed. 
Thus, although the crystals contained cracks at relatively small strains 
they were able to continue deforming under tension. 

Interesting direct observations of the development of slip and the forma- 
tion of stabilized cracks in a crystal of this category have been recorded on 
motion film. Figure 5 reproduces selected frames from the sequence. 
There were two regions of interest. The first was the group of slip bands 
represented by A on the right-hand side of the figure. Between frames (q) 
and (h), which were approximately 1/16 second apart, a stable (110) 
microcrack could be seen to have formed which was bound at each end by 
aslip band. The crack grew in steps represented by frames (i), ( j), and (hk), 
being arrested as it crossed each (110) slip band in its path. In the frame 
next after (k) the crystal had completely fractured. 

The other region of interest was that associated with the slip band B. 
In frames (a) to (e) the stress concentration associated with this 
band increased in intensity, while in (f) it decreased and in (g) had dis- 
appeared altogether. At the same time at the lower right extremity of the 
band a tiny crack formed, as in (f) (just above arrow) and subsequent 
frames. This crack was completely stable and did not increase its dimen- 
sions throughout the rest of the deformation. 


} Figures 2-6, 9 and 10, 12-15 are shown as plates. 
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There can be no doubt from these observations that slip bands have the 
ability to halt and stabilize cleavage crack nuclei and that this is an impor- 
tant parameter in the ductility of magnesium oxide crystals. 


3.3. Category III 


It was possible in some crystals of magnesium oxide for a single slip 
band to grow uninterrupted in width until it had covered the gauge length 
from end to end. In general, these slip bands were observed to originate 
at the stress concentration near one of the grips and spread from one end 
to the other. However, examples were seen where a single slip band 
nucleated at the centre of the gauge length and spread in both directions 
towards the grips, also where a number of slip bands active on the same 
plane merged together to form a single band. 

The slip bands could be seen to expand through a low-power microscope 
while the crystal was being deformed. The tilt in the crystal surface where 
each edge of the slip band emerged (drawn on an exaggerated scale for 
band Y in fig. 16 (c)) appeared as a fine line which moved along the specimen. 
This was analogous to the appearance of moving low-angle boundaries. 
as observed in zine by Li et al. (1953). If the orientation of the slip band 
was such that a step was not produced on the surface under examination 
(i.e. when viewed in a direction perpendicular to the Burgers vector) then 
the profile of the slip band could be detected with polarized illumination. 
The stress concentration at the edge of the band then produced a marked. 
stress birefringence. Records of expanding slip bands in magnesium 
oxide for both orientations have been obtained on motion-picture film. 
The distribution of slip in a specimen which had started to deform in this 
way but which parted prematurely from the grips is shown in fig. 6. The. 
density of etch pits within the band on the right-hand side of the photo- 
graph was so high that the surface appeared to be uniformly stained. 

On the left-hand side, however, a low-angle boundary and isolated dis- 
location etch pits could be distinguished in the region which had not yet 
deformed. 
As indicated in fig. 1, crystals of category III deformed by the greatest 
amount. Unfortunately, in all instances they parted from the grips when 
‘the gauge length was filled by the slip band. It was possible that the 
stress concentration in the grips or the surface tilt to accommodate the 
slip caused the epoxy resin to relax. It would be interesting to know just 
how much these crystals would elongate and by what mechanism they 
would continue to deform and fracture within the gauge section. The 
reason for the increased ductility may be understood in terms of the 
fracture mechanism outlined in §1. When shp was confined to a single 
plane no intersections developed and one of the conditions for crack 
nucleation by this mechanism did not exist within the gauge length. 

The reason why the number of intrinsic slip sources operative from the 
onset of plastic flow should be different for crystals of the different categories 
is not understood. It is possibly related to the impurity amount and. 


16 R. J. Stokes et al. on 


distribution, i.e. whether in solid solution or precipitate, although no 
correlation between impurity content and characteristic behaviour has 
been found. In general it was true that specimens originally prepared 
from the same bulk crystal fell into the same category and by far the 
greatest number behaved as in category I. 


§ 4. Stipe Sources INTRODUCED BY CLEAVAGE 


We have compared the mechanical behaviour of crystals having carefully 
cleaved surfaces with the behaviour of companion crystals having chemi- 
cally polished surfaces. Tensile stress-strain curves for the cleaved 
crystals are reproduced in fig. 7; their innate mechanical behaviour fell 
in the category indicated. Since the slip sources introduced by cleavage 
were limited in number by the distribution and nature of the tear lines, 
this treatment essentially put all crystals, irrespective of their origin, into 
category II. In this respect, the most interesting change in behaviour 
was shown by the crystal of curve 1. Its ductility should have been 
approximately zero (category I) but this has now been increased to 6%. 

Optical examination of cleaved specimens after fracture revealed that 
many stable (110) [100] slits.had formed during deformation and that 
one of these had propagated for fracture to give the characteristic line 
source of fig. 3. 

Fig. 7 


CATEGORY OF 
SPECIMEN COMPANION CRYSTAL 
I I 
2,3,4 0 


2.0 -3.0 4.0 
ELONGATION ~PERCENT 


Tensile stress strain curves for magnesium oxide crystals having carefully 
cleaved faces. 


Crystal 1; companion of erystal of category I, fig. 1. 
Cleaved on one face only. 


Crystal 2, 3, 4; companion of crystal of category IT, fig. 1. 
Cleaved on two opposite faces. 


§ 5. Strep Sources InrropucEeD By SPRINKLING 
The observations of §$3 and 4 suggested that if the number of slip 
sources could be increased many-fold by sprinkling the surface uniformly 
with fine silicon carbide powder, then the distribution of slip bands would 
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be made so dense as to fill the whole crystal and suppress the formation 
of cracks by the conventional mechanism described in $1 altogether. 
This was found to be the case for crystals deformed in bending, but in 
tension the same procedure led to rather unexpected results. For 
continuity the tension tests will be presented first. 


5.1. Tension Tests 


Crystals which had been sprinkled were found to be consistently ductile 
in tension irrespective of their origin. Figure 8 reproduces the stress— 
strain curves for a number of specimens and shows how the enhancement 
of ductility was approximately the same whether one surface, two opposing 
or all four surfaces were impacted. 


Fig. 8 
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Tensile stress-strain curves for magnesium oxide crystals which have been 
sprinkled with silicon carbide on (a) one surface, (b) two opposite 
surfaces, and (c) all four surfaces. 

Crystal 1, 2; companion of category ITI, fig. 1. 
Crystal 3, 4, 5,6; companion of category L, fig. 1. 
Crystal 7; companion of category II, fig. 1. 


The rather surprising feature of the deformation was the manner in 
which the slip was distributed along the crystal. The gauge length became 
sub-divided into a number of blocks within which slip was very dense and 
confined to a single slip plane. In the individual blocks slip could take 
place on one of the four {110} planes available, although it was generally 
found along the gauge length that slip alternated between the two ortho- 
gonal planes of a conjugate pair. At the boundary between the blocks 
(across which the slip plane changed) the slip was less dense and the two 
systems were able to interpenetrate one another. This is illustrated for 
the two blocks in fig. 9(6) and (c). Only a few isolated bands were able to 
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penetrate completely across the boundary into the region of dense slip 
on the orthogonal system on the other side. Significantly, slits were only 
to be found in those regions of less dense interpenetrating orthogonal 
slip located at the block boundaries as shown in fig. 9 (c). 

While a few specimens fractured within the gauge length by the 
propagation of a slit formed in a block boundary, more typically fracture 
occurred in the vicinity of the grips where stress inhomogeneities did not 
favour the activation of dense single slip. As may be expected, the 
ductility and microscopic aspects of slip varied somewhat for different 
crystals, depending upon the number of blocks developed and the relative 
orientation of slip on either side of a block boundary. The number of 
blocks formed was independent of the number of surfaces sprinkled so that 
one can readily appreciate why sprinkling one surface was just as effective 
as sprinkling all four in promoting ductility (fig. 8). 

Block formation is a well-known phenomenon in sodium chloride single 
crystals when compressed (Taylor and Pratt 1958) or deformed in tension 
(Webb 1960). It is not properly understood why a crystal which contains 
many sources for slip on all four equally stressed planes should prefer to 
confine slip to one plane within a given block, although this could be the 
consequence of a non-uniform stress distribution due to mis-alignment or 
variation in cross section along the gauge length. 


5.2. Bend Tests 


In contrast to the spectrum of mechanical properties observed in 
tension (§3), polished magnesium oxide crystals behaved more consis- 
tently when deformed by bending. Such crystals always slipped on the 
two {110} planes having the bending axis as their zone axis. They were 
fairly ductile, formed stable slits and fractured cleanly in two from a 
‘line’ source (Stokes et al. 1959 a). Sprinkling these specimens to increase 
the density of slip resulted in a three to five-fold improvement in ductility 
but in addition, and more significant, there was a profound change in 
the fracture behaviour. 

In the first place, sprinkling did not result in the formation of slip 
blocks as was the case under tension. Instead slip on the two intersecting 
{110} planes was made very dense in the deformed zone at the centre of 
the specimen as shown in fig. 10. Second, stable (110) [001] slits were 
never found, at least with the resolution of the optical microscope, either 
in the crystals during deformation or in the fractured pieces. Finally, 
when the crystals fractured they shattered to leave a very rough and 
complicated fracture surface. 

In order to understand this fracture behaviour we should consider the 
change in shape which occurs when a crystalline beam is severely bent. 
In a truly elastic body the tension surface becomes curved about an axis 
parallel to the length of the beam as shown in fig. 11(a). This curvature 
is due to the lateral contraction of the beam. It is termed anticlastic 
curvature (Timoshenko 1934) and is proportional to Poisson’s ratio. In 
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plastic crystalline solids, at least those of the rock-salt structure, this 
curvature is accommodated by a kinking mechanism, which will be 
referred to as anticlastic kinking. This is illustrated diagrammatically in 
fig. 11(6). The kinks produced in silver chloride, sodium chloride and 
lithium fluoride bent at room temperature and above are so pronounced 
they can be seen with the naked eye and felt with the finger. Their forma- 
tion will be described elsewhere. 

In magnesium oxide anticlastic kinking at room temperature was a 
more difficult process and by resisting lateral contraction the transverse 
stress at the tension surface was not relaxed. Consequently when cracks 
nucleated in the severely bent sprinkled crystals (by a mechanism not yet 
understood) they propagated parallel with the length of the crystal. 
Figure 12 is a photomacrograph of the tension surface of a specimen which 
bent to the point where anticlastic kinking could in fact just be detected 
in certain regions along the edge of the specimen (fig. 13). It will be 
noticed that the edges of this specimen have for the most part shattered 
instead of kinking and that longitudinal cracks have formed towards the 
centre of the specimen. 

Fig. 11 


TENSION SURFACE 


A 


NEUTRAL AXIS 
(a) (b) 


(a) Anticlastic curvature in an elastic solid. 
(6) Anticlastic kinking in a plastic crystalline solid. 


More typically longitudinal cracks developed before any anticlastic 
kinking occurred whatsoever. Figure 14 shows a longitudinal crack at a 
high magnification. The crack surface was very complex ; it contained 
many forked offshoot cracks which themselves were likewise branched 
many-fold. Focusing into the crystal beneath the tension surface it was 
found that the offshoot cracks visible on the surface did not penetrate 
very deep but were replaced by other branched offshoots from the main 
crack surface. The longitudinal crack surface was neither flat nor charact- 
erized by any readily identifiable crystallographic planes, although the 
branched offshoot fingers tended to point in either [100] or [110] directions. 
Attempts were made to observe the formation of a longitudinal crack 
directly, by loading sprinkled crystals in a bending jig ea NE! or the 
specimen stage of the metallograph. It proved impossible, however, to 
detect them at any stage of development before fracture and it was pre- 
sumed that they formed more or less instantaneously. Slip traces parallel 
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to the tension direction were detected in some cases, indicative of the 
existence of the transverse stress component. Examples may be seen in 
the region of the crack in fig. 14. 

It will be noted that there was a step at the point where the longitudinal 
crack of fig. 14 intersected the specimen fracture surface. HExamination 
of the fracture surface, fig. 15, confirmed that the step did in fact correspond 
with the origin of fracture. In other words these observations taken 
together indicated that the longitudinal crack was formed first in the 
erystal parallel to the tension direction, that it generated offshoot cracks 
which lay in some cases perpendicular to the tension direction and that 
one of these offshoot cracks was eventually able to propagate through the 
crystal for fracture. The appearance of the fracture surface which was 
very rough, and the high number of offshoot cracks which did not propagate 
for fracture, strongly suggest that crack propagation through the densely 
slipped material was a very difficult process. 


§ 6. BEHAVIOUR UNDER Impact 


It was of interest to determine whether densification of slip would 
produce a similar enhancement.of mechanical properties in a test conducted 
at a higher strain rate, especially in view of the strong strain rate sensitivity 
of ionic solids. Consequently, a comparison was made between the energy 
absorbed in impact by chemically polished crystals and crystals which 
had been sprinkled, pre-strained a slight amount, and repolished. It was 
necessary to preload the crystals at four different points along the beam 
to ensure that the whole volume to be loaded under impact was pre- 
strained uniformly. In fact, fig. 10 was taken on a specimen prestrained 
for use in the impact study. The table includes a comparison of the energy 
absorbed for these two states. It can be seen that dense slip generally 
increased the energy absorbed in fracture under impact. Indeed, one 
of the specimens tested actually bent further without fracture and was, 
therefore, able to absorb all of the impact energy (10 in. lb) without fracture. 
This was very surprising for it has been shown that polished specimens 
normally have to be heated to 1900°c to bend under impact (Johnston 
et al. 1959). 


Kffect of dense slip on energy absorbed in impact 


Knergy absorbed (in. Ib) 


Polished-sprinkled 


IN 7 1a - 
8 polished Preloaded—repolished - 


0-31 
0-42 0-78 
0-60 0-81 
0-61 3°88 
5:26 


~ 10-0 bent without fracture 
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§ 7. THe Use or Stress BIREFRINGENCE TO DETERMINE SLIP 
DISTRIBUTION IN Maanestum OxipE 


At this point it is appropriate to comment upon the interpretation of 
stress birefringence patterns with respect to the distribution of slip in 
magnesium oxide. It was found that stress birefringence was produced 
only by the stress concentration at the edge of a slip band. Thus, in 
crystals of category IIT (§3) where slip was confined to a single wide band, 
the birefringence pattern delineated the two edges of the band only, 
within the band extinction was practically complete. Similarly, in fig. 9 (a) 
stress birefringence was most pronounced at the block boundary where 
narrow slip bands intersected one another and along the isolated bands 
(indicated by arrows) which penetrated across the block boundary. In the 
region on the left-hand side of fig. 9(a) where slip was most dense but 
confined to a single plane there was no stress birefringence pattern. When 
the slip was dense and on intersecting planes, as in fig. 10, then the stress 
birefringence pattern gave a reasonable indication of the slip distribution. 

Thus the generally accepted correlation between the birefringent 
pattern and the lines of dislocation etch pits established by Kear and 
Pratt (1958) for lithium fluoride and by Shaskol’skaya and Tui-fang 
(1960) for sodium chloride did not apply for magnesium oxide unless the 
slip bands were sufficiently narrow. Indeed in certain instances, fig. 9 (a) 
for example, stress birefringence gave almost the inverse pattern to that 
expected, being strongest in the regions which had deformed the least. 
It is probable that the reason for the lack of stress birefringence within 
a slip band is due to the cancellation of the stresses of the positive and 
negative edge dislocations retained there. 


§ 8. DIscussION 


While the location and orientation of slits described in §§3, 4 and 5.1 
of the present work are consistent with results of the previous paper 
(Stokes et al. 1959a), there are certain additional observations which 
cannot be accounted for in detail by the fracture mechanism as proposed 
there and outlined in §1. It will be remembered that there are two 
stages in slit formation, first the nucleation and second the growth 
and stabilization stage. With regard to stabilization, it 1s considered 
that the original interpretation is correct. The new evidence supports 
the original contention that a crack is held up by a slip band, probably 
because of the sudden increase in surface energy when it encounters the 
high density of screw dislocations retained in the band. The details of 
the nucleation mechanism, however, are not so clear and the following 
additional observations must be accounted for by the model. 

There are frequent occasions where a slip band bursts across the crystal 
to be held up by an intersecting band within the gauge length. An example 
is indicated by the arrow in fig. 5 (d). Even though this event : equivalent 
to a sudden avalanche of edge dislocations piling into an obstacle a crack 
does not form. Instead the slip band eventually breaks through the 


22 R. J. Stokes et al. on 


intersection, fig. 5(e), and crosses the crystal. Indeed cracks are only to 
be found at intersections where two slip bands have completely cut through 
one another, and are generally associated with slip bands formed in the 
very early stages of plastic flow (fig. 5). The reason why only these slip 
band intersections nucleate a crack remains a topic for more detailed 
investigation. 

The situation of the slit in fig. 4 warrants further discussion. It may be 
interpreted with the aid of fig. 16. Consider two independent slip bands 
as in fig. 16 (a), growing in width to meet and intersect one another. From 
the sign and direction of motion of the edge dislocations it can be deduced 
that a crack can nucleate at A (fig. 16(c)) by the piling up of dislocations 
moving in band Y against the barrier provided by band X. It will be 


The formation of a crack at the intersection of two expanding slip bands to 
give the distribution of slip shown in fig. 4. 


noted that this line of intersection does not exist (fig. 16 (b)) until the edge 
of band Y has cut completely through X. Once the crack has nucleated 
it can propagate along the edge of (in the [110] direction) and across (in 
the [100] direction) slip band X, but at the same time slip band Y continues 
to grow in width until it is finally cut off by the fracture path. This leaves 
the distribution of slip shown in fig. 16(d) which corresponds with the 
distribution of etch pits shown in fig. 4. This mechanism explains why 
the line source of fracture is so close to the edge for crystals in this category 
(fig. 3). However, there is a weakness in the interpretation in the necessary 
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assumption that dislocations moving at the edge of slip band Y pile up 
against slip band X. These dislocations move only a short distance to the 
intersection at A, and slip band X, the barrier, is the narrower slip band. 
of the two. 

The complete suppression of slit formation in bent erystals for which 
intersecting slip has been densified (§ 5.2) may be influenced by two factors. 
First, it is possible that the change in behaviour is a consequence of the 
stabilization mechanism. If the spacing between slip bands is reduced 
below microscopic dimensions by sprinkling, then the slits, even though 
they may form, will not be visible and, being so small, will not exceed the 
critical size for catastrophic propagation across the crystal under the 
applied tensile stress. As fig. 14 and the behaviour under impact (§ 6) 
indicate, the propagation of a crack through a uniformly strained crystal 
is a difficult process. 

Second, it is possible that the crack nucleation process is upset when 
intersecting slip occurs on a very fine scale. Dislocation sources introduced 
by sprinkling are able to expand and multiply to form short slip line 
segments uniformly throughout the crystal at quite low stresses compared 
to the macroscopic yield point (Stokes et al. 19596). As plastic deformation 
proceeds the slip bands develop gradually until slip more or less saturates 
the crystal. For polished crystals, on the other hand, the number of 
intrinsic sources is limited and most of the plastic strain at any instant is 
confined to the edge of a few isolated slip bands. The stress concentration 
in this region is more conducive to the dislocation avalanche which may 
be involved in crack nucleation. The effect of sprinkling then is to change 
the manner in which slip bands develop in the crystal. We may regard 
the high number of slip sources as effectively homogenizing the strain and, 
therefore, the stress through the crystal. 

The observations of this paper indicate an additional factor which 
should be considered in discussion of the apparent difference in ductility 
between crystals strained in bending and in tension. It is frequently 
possible with a polished (unsprinkled) crystal to obtain as much as 10% 
surface strain under bending while an equivalent strain has never been 
produced in tension. This difference might be connected with the distri- 
bution of slip in the two instances. Under tension there is a large volume 
of material under stress and the chances of a crack forming at an inter- 
section before sufficient stabilizing slip has been introduced are greater 
than under three point bending where the highest stress is confined to the 
surface of a narrow region at the centre of the crystal. It is then possible 
for slip from many sources on orthogonal systems to develop, so that 
when cracks form they are immediately stabilized and the crystal continues 
to bend. 

An important aspect of the dislocation theory of fracture is concerned 
with understanding: how much plastic deformation precedes brittle 
cleavage fracture and whether it is the nucleation or the propagation of 
fracture which is the critical stage. On the basis of these observations on 
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magnesium oxide it is obvious that the amount of plastic deformation 
preceding cleavage fracture depends primarily upon the density and 
distribution of slip sources at the onset of plastic flow. It depends also 
upon the chance that a particular slip band intersection will nucleate a 
crack, which, until more is understood of the nucleation mechanism, 
should be regarded as statistical in nature. Whether crack nucleation or 
propagation is the more critical stage depends again upon the density of 
slip sources and when the crack forms. If it nucleates before other slip is 
present to stabilize it then obviously nucleation is the critical stage. 
Whereas if a crack is formed and stabilized then it is the propagation 
through the surrounding slip bands which is the critical stage of fracture. 
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SEGREGATION of impurities to grain boundaries has been described for a 
number of crystalline solids of different structures (Amelinckx and 
Dekyser 1959). Etching, autoradiography, ‘ decoration’ and electron 
microscopy are examples of methods used in the investigation of such 
segregations. In a phosphor, if segregation of an activating impurity 
occurs, it should be visible under the microscope, the impurity acting, in 
effect, as an optical tracer. Segregation of silver ions to dislocation walls 
in a strained silver-activated rocksalt crystal has been demonstrated in 
this way (Van der Vorst and Dekyser 1956), luminescence photographs 
showing diffuse lines corresponding to the dislocations. This note 
describes segregation of impurity to grain boundaries of tin-activated 
potassium bromide during growth from the melt. 

The impurity absorption band for this phosphor has its maximum 
intensity at 2650 A with the corresponding emission band lying completely 
in the visible region. Polycrystalline specimens were grown from a melt 
of potassium bromide doped with 1% by weight of stannous bromide and 
held in a quartz beaker. A conventional tube furnace was used, growth 
being in air at atmospheric pressure. A brass rod, held in a water-cooled 
brass chuck, was lowered into the melt which was at a temperature only 
just above its freezing point. The furnace temperature was then allowed 
to drop to below the melting point. As the furnace cooled, the surface of 
the melt cooling faster than the bulk of the melt, growth suddenly pro- 
ceeded outwards from the brass rod along the melt surface. This ‘ skin ’ 
was lifted off just before it reached the sides of the beaker, total growth 
times being less than one minute. Grain dimensions were about 0-5 mm 
in width with thicknesses from 0:5 mm up to 2 mm, each grain extending 
through the whole thickness of the polycrystalline sheet. Quite regular 
hexagonal networks were often formed, the grains usually having a (100) 
direction close to the horizontal plane. ; 

Figure 1{ is a luminescence photograph of a fragment of polycrystalline 
sheet grown as described above and immediately cooled to room tempera- 
ture. Unfiltered radiation from a mercury low-pressure discharge tube 
stimulates visible emission. Figure 2 is shown for comparison purposes, 
the ultra-violet radiations of the lamp being removed by placing @ glass 


+ Communicated by Professor J. Sayers. 
+ Figures 1, 2 and 3 are shown as a plate. 
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plate between lamp and specimen. Figure 2 shows no sign of the slightly 
diffuse luminescent network of fig. 1. Exposure times were 5 min using 
H.P.S. roll film, exposure and development times being the same for the 
two photographs. Figure 3 shows the surface structure revealed by 
transmitted light, the grain boundaries forming two similar networks 
where they meet the two surfaces of the specimen. Because of the small 
magnification and large depth of focus, both surfaces are well defined in 
fig. 3. Comparison of fig. 1 with fig. 3 shows that the luminescent net- 
work corresponds to where grain boundaries cut the top surface of the 
specimen as it was when removed from the furnace. The good contrast 
~ of fig. 1 indicates that the light emitted from near grain boundaries cannot 
be due to preferential scattering of light emitted from the bulk of the 
grains, but is, in fact, due to a non-uniform distribution of activating 
impurity. Annealing at temperature greater than 300°c followed either 
by quenching or slow cooling to room temperature caused the luminescent 
network to disperse. 

These specimens were prepared while work was being carried out on the 
growth of single crystals of tin-activated potassium bromide from the 
melt. In the latter case the segregation coefficient of tin in potassium 
bromide was found to be too small for it to be possible to grow mixed 
erystals containing any significant quantities of tin. It is, therefore, pro- 
bable that as a polycrystalline sheet grew quickly over the melt surface, 
impurities rejected by the advancing growth surface formed a concentra- 
tion gradient into the melt at the solid—liquid interface. Right at the 
surface of the melt where the degree of super-cooling was greatest, im- 
purities at the growth surface, while not accepted by the bulk crystal, 
could enter the solid phase at disordered regions near the boundaries. 
Conditions for the occurrence of the segregation described in this note were 
found to be very critical, only a few of the examined specimens showing 
luminescence associated with the boundaries. Closer control of the growth 
conditions will be necessary before the precise mechanisms of solidification 
and luminescence can be understood and reproducible specimens prepared 
for study. 
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ABSTRACT 


A discussion is given of the assumptions involved in recent theories of the 
dielectric constants of the alkali halide crystals with particular reference to the 
contributions arising from short-range polarization effects. A theory is 
developed on a semi-empirical basis which makes it possible to divide this 
effect unambiguously between the positive and negative ions, This theory is 
consistent with the experimental data and should provide a means of 
including this effect in calculations on the lattice dynamies of these crystals. 


§ 1. INTRODUCTION 


In two papers on the dielectric properties of ionic crystals Szigeti (1949, 
1950) deduced, from a study of long-wave optical vibrations in: these 
crystals, the following pair of equations: 


Meo,2(€y+2)|(€m +2) = 6R,/K, eee nee 11) 
Ey — Eco = (€x + 2)?47Ne?/9Mw,?. yb) ere tee 0) 


These relate the empirical static and optical frequency dielectric constants, 
€, and ec, respectively, to the corresponding compressibilities and 
restrahlen frequencies K and w,). The other parameters are Ry the 
equilibrium nearest-neighbour separation, M the reduced ionic mass 
and e the effective ionic charge. 

The derivation of (1) and (2) (Born and Huang 1954, denoted by A 
hereafter) leads one to expect that |e|=1 electronic charge. In the table 
(after Szigeti) the values of K and |e| (A* and e*) calculated from eqns. 
(1) and (2) are compared with the actual values. It is at once clear that 
the first relation is obeyed rather better than the second, particularly by 
the lighter salts. Furthermore there appears to be little correlation 
between violations of the two relations. This suggests very strongly 
that at least two distinct effects are present and that any theory should 
allow for this fact. However, the data used in calculating the results in 
the table are all room-temperature data, while the various theories of 
these effects are strictly only true for crystals at absolute zero. It is 
thus a matter of some importance for measurements of the relevant data 
to be made at low temperatures. This applies particularly to the 
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restrahlen frequencies wy and the static dielectric constants €«. (There 
are two sets of the latter available (A, p. 85, Haussiihl 1957), both 
measured at room temperature but which differ by as much as 10%.) 
Thus the data do not make it possible to test any theory very rigorously. 
What is certain however is that the theory which led to eqns. (1) 


Salt ELK e*, e| 
Lif 1-0 0-87 
NaF 0-83 0-93 
NaCl 0-99 0-74 
NaBr eile: 0-69 
Nal 1-05 0-71 
KCl 0-96 0-80 
KBr 0-95 0-76 
KI 0-99 0-69 
RbCl 0-89 0-84 
RbBr 0-83 0-82 
RbI 0-66 0-89 


and (2) must be modified in some way. Such a modification will affect 
not only the long-wavelength lattice vibrations, but the whole spectrum, 
and the problem which faces us is that of attempting to predict this 
effect. We want to do this, if possible, in such a way that the results are 
not too strongly dependent on the detailed mechanisms responsible 
for the failure of the Szigeti relations. Im this way we may hope that 
the results of the analysis will remain valid as more exact theories 
are developed. 


§ 2. Bastc ASSUMPTIONS 

The analysis (A,p. 100) which leads to the two Szigeti relations is 
based on the following assumptions: 

In the case of the first relation: 

(a) That the effective field acting on an ion in a uniformly polarized 
crystal of polarization P and in a field E is 

E+ $nP. 
(This only applies to crystals with cubic or tetrahedral symmetry.) 

(>) That the repulsive interactions due to ionic overlap act only between 
nearest neighbours and depend only on the absolute magnitude of their 
separation. 

(c) That there are no non-central interactions. 

The second relation requires only (a) to be valid, but involves one further 
assumption : 

(d) That the dipole moment at an ion site when the ion in question is 
displaced by a small distance — is +e€ (apart from the electronic 
polarization). The failure of the second relation implies that either 
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(a) or (d) is invalid. Tessman et al. (1953) have carried out an analysis 
of the optical frequency refractive indices of these crystals and have shown 
that it is possible to derive a set of polarizabilities for all the ions such 
that the value for a given ion does not depend on its environment. This 
work was also based on assumption (a), and they found, when they tested 
this by varying the Lorentz factor 47/3, that any other value produced a 
less consistent set of results. It thus seems fairly certain that assumption 
(d) is false. It is therefore necessary to explain this breakdown, if possible, 
in a way which does not imply any violation of (6) or (c). Szigeti (1950) 
himself suggested that the mechanism in question is a distortion of the 
ions associated with the overlap repulsion. Born and Huang (A, p. 112) 
have extended this idea and shown that if each nearest-neighbour bond 
has associated with it a dipole moment depending only on the absolute 
separation of the two ions then the effective charge 


ee lel an tay 4 tey)j ley 2 te Meee ne a(S) 


where m(r) is the ‘deformation’ dipole associated with each nearest- 
neighbour bond of length r, m’(r) =[dm(r)/dr], and R, is the equilibrium 
separation of the ions. This derivation is quite general and independent 
of the details of the mechanism involved. Furthermore, it does not 
automatically imply the failure of any of the other assumptions. 
However, before one can set up a complete phenomenological theory, 
based on the empirical data, one needs to know both m(r) and m’‘(r) 
separately and where the dipoles should be located within the unit cell. 
The theory leading to relations (1) and (2) gives us neither piece of 
information; it only determines e* and, for the homogeneous polarizations 
involved, the exact location of the dipoles is unimportant. It is thus 
necessary to have a more specific understanding of the effect. 


§ 3. THEORY OF THE DIELECTRIC CONSTANTS 


There have recently been two treatments of this problem (Dick and 
Overhauser 1958, denoted by B hereafter, and Hanlon and Lawson 1959) 
which have proposed detailed models with ‘short-range’ polarization 
mechanisms incorporated into them. Of these the first treatment is the 
more complete and will be our main concern here, and this paper follows 
their notation. However, both treatments are open to the same objection ; 
they both employ a ‘shell’ model in which the ions are regarded as 
massive cores isotropically coupled to rigid massless shells consisting of 
the valence electrons. Ionic polarization then appears as a relative 
movement of core and shell and, in antiphase motion, the positive and 
negative ions, by virtue of the overlap repulsion between shells, induce 
dipoles on each other which do not in general cancel and are the main 
cause (in the Hanlon—Lawson theory the only cause) of the deformation 
dipole moment. This feature implies that assumption (6) of the Szigeti 
theory is invalid since the repulsive potential between two ions depends 
not only on the internuclear distance but also on the relative core-shell 
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displacements. This leads one to expect not only the failure of the first 
Szigeti relation, but also that K*/K will always be greater than unity by 
an amount which should depend almost linearly on the sum of the positive 
and negative ion polarizabilities. Both these results are obviously not 
true of the data in the table. It is rather surprising that K*/K is 
generally <1 since one would expect short-range interactions due to 
overlap between other than nearest neighbours, which are the most 
obvious cause of deviations from unity, to be repulsive, in which case 
K*/K would be >1. The fact that this is not so seems to imply that the 
Van der Waals attraction between nearest like neighbours is the more 
important term. 

The theory of Dick and Overhauser does, however, include another 
short range polarization mechanism; the ‘exchange charge’ polarization. 
This effect arises naturally when the charge density in the crystal is 
derived using a Slater determinant of free ion orbitals. These are 
non-orthogonal in the crystal and this fact is responsible for an additional 
term in the charge density; an ‘exchange’ term: 


= 2|e|DS,,U,(r)U;(r) 
Wy 
where U; and U; denote orbitals centred on nearest neighbour ions and 
Sy =JU(r)U,(r)dr. 


The total localized charge q,, is exactly compensated by a spherically 
symmetric increase in the free ion charge densities whose total magnitude 
is —q,,/2 for each ion. In a uniformly polarized crystal there is a 
redistribution of exchange charge which gives rise to a net dipole moment. 
Furthermore, in the unpolarized crystal, this charge distribution has cubic, 
rather than spherical, symmetry about each ion. Dick and Overhauser 
have estimated the contribution of the exchange charge polarization to 
the value of |e*/e|. A lower limit to |e*/e| is given by assuming this charge 
to lie at the centre of the positive ion, but, even when this is done, most 
of the calculated values are still about 10% too large and about half the 
calculated deviation from unity is due to the relative motion of cores and 
shells. If this last effect is ignored, in view of the way in which the 
observed data satisfy the first Szigeti relation, then (1—|e*/e|) is probably 
<10% for most of the alkali halides. (An earlier, but possibly more 
exact, calculation by Lundqvist (1955) gave (1—|e*/e])=8% for NaCl.) 


§ 4. CONTRIBUTION OF IonIc DEFORMATION TO THE EFFECTIVE 
CHARGE 


There is one possibility which is not included in the Dick and Overhauser 
treatment; this is a redistribution of charge within an ion owing to a 


marked difference between the outermost parts of free ion and crystal 
wave functions. 
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There are two pieces of evidence for this view : 

(a) The fairly large differences between the polarizabilities of free ions, 
as calculated by Pauling (1927) and the T.K.S. crystal polarizabilities. 
Dick and Overhauser have attempted to account for these in terms of 
the exchange charge, but the predicted values are too small by about a 
factor of five. 

(b) The contour maps of electron density as measured by Witte and 
Wolfel (1955) for LiF and NaCl. These show a pronounced cubic 
symmetry, but a theoretical calculation by Howland (1957) for KCl 
using free ion functions fails to reproduce the distinct ‘bulges’ of electron 
density along the (110) directions which occur in the experimental data 
for NaCl. He explains these as dynamical effects arising from ionic 
thermal motion along preferred directions. However, in LiF the 
‘bulges’ lie along (100) directions and it seems more likely that they 
result from a genuine distortion of the ions, which is not represented by 
this theory, owing to the change in the outer parts of the wave functions. 

These facts provide support for the deformation hypothesis. If this 
hypothesis is true it would seem that most of the deviation of |e*/e| from 
unity is due to this effect, which does not imply any violation of the first 
Szigeti relation. 


§ 5. Iontc DEFORMATION AND LATTICE VIBRATIONS 


In the theory of lattice vibrations one has to be able to deal with 
lattice distortions which change markedly over a few interionic distances. 
This was first done for an ionic lattice by Kellermann (1940) for NaCl. 
He deduced a lattice vibration frequency spectrum allowing for the 
long-range Coulomb interactions, but regarding the ions as rigid point 
charges. This work has recently been repeated and extended to other 
alkali halides by Karo (1959, 1960) and the agreement with experimental 
heat capacity data (e.g. Berg and Morrison 1957) is very good. However, | 
in another recent paper (Hardy 1959) it was shown that the effect of ionic 
polarization, when it is included in the NaCl theory, is to destroy this 
agreement. If deformation dipole effects are included, the agreement 
can be restored even though the resulting frequency spectrum is very 
different from that predicted by the simple theory. 

In this particular case the entire deformation dipole effect was assumed 
to’ be localized on the Cl- ion owing to its large ‘size’ relative to that 
of the Na+ ion. However, in order to construct a general theory, 
applicable whenever the first Szigeti relation holds, it is necessary to 
split up the deformation dipole moment m(r) associated with each 
nearest-neighbour bond into two parts: m,(r) at the centre of the 
positive ion and m_(r) at the centre of the negative ion, each moment 
directed outwards from the centre of the ion in question. This division 
is implicit in the Dick-Overhauser theory, which has been used by 
Cochran (1959) to calculate frequency-wave number dispersion curves 
for lattice waves propagating along <100) directions. The results are 
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in good agreement with those obtained during the previous work (Hardy 
and Karo 1960) but this might not be so if the positive ion were larger. 

In §3 we referred to the exchange charge q,, defined by Dick and 
Overhauser. To estimate the strength of the dipole associated with each 
nearest-neighbour bond they assumed that the total exchange charge 
between each ion pair is concentrated at the point of tangency of spheres 
centred on each ion, and having the appropriate crystal radii r, and r_ 
defined by Zachariasen (see Kittel 1953, p. 40). Thus 


f= = lel Sgt oe eae ieee 
2 5 
One would certainly expect the maximum exchange charge density to lie 
in this region and to be fairly strongly localized there. 

In the absence of a set of true crystal wave functions it would seem 
reasonable to assume that the change in charge density, due to ionic 
distortion, is proportional to that predicted from the overlap of the free 
ion functions. However, in B it is shown that: 


- S,,? CAR, exp (— Ry /p) 
Yy 
where R,=r,+r_ is the equilibrium nearest-neighbour distance and p 


and 2 are the constants in the overlap repulsion. 


Thus we may write, quite generally, that 
m..=dr,h exp (— Ro/p) (5) 


where d, are disposable parameters. Hence, from eqns. (3) and (5) 


1—letJe|=28,\| “—F |[d_r_—d,r,Jexp(—Rolp) 
p Ro 


or ; 
as (d,—d_) Rae > 
1 jet/e|=2a_| r_?—7,2— aa rR, || -— FZ, \exP (~ Bale) (6) 
If our initial assumption regarding the proportionality of the total 
displaced charge and the exchange charge is correct, then in a given 


compound d,=d_. Let us now make the further assumption that d=d_ 
is the same for all these salts, then from eqn. (6) we have ; 
1 —le*/eloo nr =7.9(2 - a (— R,/p). 
p ky 

In the figure the resulting points are those without*. The K+ and Rb+ 
halides (except RbI) lie on a straight line through the origin of slope 
d, =d_=23-5 x 10'8(ergem)~. Thus the two initial assumptions provide 
us with a self-consistent theory of these salts. There remain the problems 
presented by Lif, RbI and the Na salts. The first presents no real 
difficulty since we assumed, following B, that >i5;;° is proportional to 
the overlap repulsion, but we should reduce the constant of proportionality 
in the ratio 3:1 for LiF to remain consistent with B. When we do 
this the resultant point, LiF** lies close to the K+Rb+ line. RbI seems 
to be anomalous and probably cannot be treated on the same basis as the 
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rest of these saltst. Thus, the sodium salts present the only real problem. 
Naf’, NaCl and NaBr lie on a straight line, but it no longer runs through 


(Na Br* also) 


0-30 


dyd_=23.5 x 10'° 
Cerg cm) for 
Kt and Rb* salts 


Jide = dia 26-0000 4 
OBSERVED) / (erg.cm)! for 
/ Na*™ salts 
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= d_ 1 aay 
1071 xA [Aes eS rR | E -= | exp (— R,/p) 


i ted out errors in 
R tly Havinga (1960, Phys. Rev., 119, 1193) has poin . 
Srpeti (1 950) eee of (1—|e*/e|) for LiF and RbI. The new values are 0-17 


and 0-21 respectively. If one uses these the Lif’ data agree less well with 


the present theory, but most of the large RbI anomaly disappears. 


C 
P.M. 
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the origin. If we now assume that the effective centre of the displaced 
charge lies at a distance r,+6r from the positive ion, then we have 


m.(r) =d(r, + 6r) RA exp (— B,/p), 
d here is assumed to have the same value as for the K+ and Rb* salts. 
Then 
2 
eS e*lel) = 2dr) -— Z| [r2—r,2—2R,6r] exp (— Ry/p); 
p fy 
which is exactly equivalent to eqn. (6) if 


d=d_ and a ae ee d,—d_(1+2™). 


igh ft , 


If we assume that Sr is the same for all the Na* salts, and then choose 
dr/r, =0-114 to fit the resulting point NaCl* to the K*, Rb* line in 
accordance with the initial assumption that d=d_ = 23-5 x 108 (erg em) t. 
The NaF and NaBr points then lie very close to this line. Nal is still 
anomalous; but to accord with this theory (1 —l|e*/e|) would have to be 
~ 0-35 for this salt, and as the theory probably breaks down when (1 — |e*/e]) 
becomes large, the-discrepancy seems to indicate that Nal is a case in 
this region. Values for the heavier Li* halides would be of interest in 
order to test this hypothesis. 

We have now found a self-consistent phenomenological theory of 
effective charge in the alkali halides which involves only two disposable 
parameters (d_=d, and dr) and describes at least nine of these salts 
well. This fact would seem to support the basic assumptions fairly 
strongly. 

The next stage of this work is to incorporate the values of d. which 
we have deduced into calculations on the lattice dynamics of those alkali 
halides which have large positive ions. The important result of this work 
is that we have now obtained a means of dividing the ‘deformation’ 
polarization unambiguously between the two types of ion, and one may 
hope, as more exact inelastic neutron scattering data become available, 
to test the theory by comparing the predicted and observed phonon 
dispersion curves, particularly in the short wave-length region. 
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ABSTRACT 


Single crystals of <100>, <111> and <110) axial orientation, and poly- 
erystals, both of 99-99% pure aluminium, were deformed in tension at 
temperatures of 27°c, 200°c, 400°c and 600°c. The experiments on single 
crystals reveal that as the temperature is increased, the shear stress—shear 
strain curves for the three orientations become more and more similar until at 
600°c they are identical. On the other hand, the surface appearance of the 
erystals strained at 600°c indicates that the mode of deformation depends 
drastically on orientation. 

Taylor’s criterion for the comparison of the stress-strain curves of single 
and polycrystals was found to be valid at 27°c. As the temperature is 
raised, however, the criterion ceases to be applicable. 

From the results of these experiments it is proposed that at the lower 
temperatures in both single crystals and polycrystals the principal hardening 
mechanism is the formation of Lomer—Cottrell barriers that act as obstacles 
against which glide dislocations can pile up. As the temperature is raised 
these barriers become ineffective as obstacles to slip, whereas grain boundaries 
in polyerystals continue to be effective in causing dislocation pile-ups. 


§ 1. INTRODUCTION 
1.1. Relation Between the Deformation of Single and of Polycrystals 


METALS with the face-centred cubic structure deform plastically by slip on 
one or more of twelve geometrically equivalent slip systems. Each 
system consists of one of the three (110) directions in each of the four {111} 
planes. . 

In a single crystal the number of acting slip systems depends on the 
magnitude of the resolved shear stresses in the various slip systems. If 
the resolved shear stress is initially greatest in one slip system, slip will 
start in that system. If the resolved shear stresses are equal on more than 
one slip system then slip occurs on some or all of these slip systems. In 
all cases the shear stresses 7 and shear strains y are related to the uniaxial 
tensile stress o and the tensile strain « by the familiar expressions 
7=0/M and y=Me where M is the ‘resolving’ factor}. 

+ Communicated by the Authors. This work was sponsored by the Office 
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ef tape 1/cos A cos ¢, A=angle between slip direction and tensile axis; = angle 


between normal to slip plane and tensile axis. 
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The situation in a polycrystal is somewhat different. A grain of a poly- 
crystal is not free to deform by single slip, for it is subject to constraints 
imposed by the surrounding grains. In order to maintain continuity 
across the boundary, a grain must deform in a particular way. Since an 
arbitrary strain is described by the six independent components of the 
strain tensor, in the most general case six independent conditions are 
necessary to produce such a strain. If one assumes that the volume 
remains constant during deformation, then slip on five independent systems 
will satisfy the remaining five conditions necessary for the grains to undergo 
any arbitrary deformation. These considerations must be taken into 
account in deriving the factor M which relates the applied tensile stress 
and the resolved shear stresses in various slip systems of an inividual grain 
in the polycrystal. 

Among many attempts to compute M for a polycrystal, the work of 
Taylor (1938) is most noteworthy. Taylor's method is based on the 
principle of virtual work, and can be summarized as follows. 

Consider a set of n shears that satisfy the continuity conditions at the 
grain boundary. The energy expended in producing the 7 shears is 


n 

> tsdys where dys is an element of shear strain in the glide system s and rg 
s=1 
is the shear stress to activate it. This energy must equal the work done 


by the external uniaxial tensile stress in producing an elongationde. Thus 
n 
ode= > Tsdys. Sea, ost on cine eek 
s=1 
In the restricted case where 7; has the same value in every slip system, 


gic = al Se ee eee ee 


This may be rewritten as 


= 2 as Reese er cee) 


In analogy with the single crystal, M is defined by M=o/r. 

Taylor argued that M can be computed by finding the combination 
of slips which produces the imposed external strain and gives the 
minimum value of >'|dy|. All other combinations of slips which satisfy the 
conditions of continuity at the grain boundaries would require a higher 
value of o than (7/de) times the minimum value of >|dy| to produce the 
required deformation. 

If one assumes that at any stage of deformation the hardening value is 
the same throughout the polycrstal, the problem reduces to finding the 
numerical value of M. For this purpose Taylor assumed that the strain 
in each grain is the same as that of the whole specimen. He further postu- 
lated that exactly five slip systems will operate in each grain and that this 
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set of systems will be the one out of 96 irreducible sets in a face-centred 
cubic crystalt which gives the lowest value of 


dy| 
Si eS >| Y i 
de 


Taylor calculated these minimum values of M for 44 orientations of the 
tensile axis with respect to the crystal orientation. Averaging these 
values for a quasi-isotropic, face-centred cubic polycrystal in tension 
gave M=3-061. 

Taylor first tried to compare the stress-strain behaviour of a polycrystal 
with a single crystal of ‘average’ orientation (whose axis was situated 
inside the stereographic triangle). Unfortunately such a comparison 
may not be meaningful, since a single crystal so oriented will deform pre- 
dominately by slip on one slip system, whereas an average grain in a poly- 
crystal will slip on more than one slip system. Kocks (1958) pointed out 
that one should compare a polycrystal with a polyslip orientation single 
crystal (one having five or more slip systems with equal resolved shear 
stresses), for such a single crystal is found to deform (by slip line exami- 
nation) in a manner similar to that of a grain in a polyecrystal. The (100) 
and (111) axial orientation single crystals satisfy this requirement. For 
reasons to be discussed later, the (111) orientation with six slip systems 
having equal resolved shear stresses is most suitable for this purpose. The 
comparison is carried out between the actual tensile stress—tensile strain 
curve obtained for a polycrystal and a fictitious tensile stress—tensile 
strain curve for a polycrystal, obtained from the stress-strain curve of the 
single crystal; this curve is derived by converting the true shear stress— 
shear strain curve of the single crystal into a tensile stress—tensile strain 
curve using the resolving factor M =a/r=3-06. 


§ 2. EXPERIMENTAL PROCEDURE 


Single crystals for this investigation were grown from the melt by a 
modified Bridgman technique developed by Chalmers (1953), permitting 
one to obtain specimens of any desired orientation. The material used 
was aluminium 99-99% pure. Specimens of (111), (100) and <110) 
axial orientation, three millimetres in diameter and about ten inches long, 
were produced by the above method. They were subsequently cut into 
two-inch long pieces with an oxy-acetylene torch to prevent deformation ; 
small spheres of aluminium were melted on the ends of the samples to 
provide a method of gripping the crystals in special slotted holders in the 
tensile machine. This method prevents slipping in the grips and pre- 
mature deformation. The polycrystals studied were prepared from hard 
drawn aluminium wire (99-99% pure) three millimetres in diameter, by 


+ In his original work Taylor took into account only 24 systems; for the 


tion see Taylor (1955). ; 
itor recent a on this subject by Bishop and Hill (1951), who used more 


elaborate assumptions, gave the same result as Taylor’s calculations, in the case 
of uniaxial stress. 
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annealing for five minutes at 620°C. This treatment produced grains 
about 0-3-0-4mm in diameter, corresponding to about 50 grains in the 
cross section. Subsequently these specimens were prepared in the same 
way as the single crystals. An x-ray examination of these polycrystals, 
using the Laue transmission method, showed complete recrystallization 
with no evidence of preferred orientation. 

After the single crystals were grown and cut into individual tensile 
specimens they were annealed for twenty-four hours at 625°c. All 
specimens were deformed in a ‘hard’ type (Instron) tensile machine. The 
various test temperatures above room temperature were achieved by 
means of a low-temperature gradient furnace mounted on the tensile 
machine. The temperature for any test was held constant to within + 1°. 

The single and polycrystals were deformed at room temperature, 200°C, 
400°c and 600°C at an elongation rate of 4x 10-°sect. In addition, 
single crystals of the (111) axial orinetation at 600°c, and polycrystals 
at all the above indicated temperatures were also extended at a rate of 
4x10-%sec-1, To determine the reproducibility of the curves, two single 
crystals of the same orientation were extended at each temperature ; only 
several percent variation was noted between two such specimens. The - 
stress-strain curves were obtained from continuous recordings of load 
versus elongation, which were replotted to show resolved shear stress— 
shear strain diagrams. 


§ 3. EXPERIMENTAL RESULTS 
3.1. Deformation of (100), <110) and (111) Single Crystals 


At room temperature (figs. 1, 2 and 3) all curves show a steeply rising 
initial portion of the shear stress—-shear strain curve. In the (100) crystal 
this is followed by a sharp bend, the stress then rising only slightly. Kocks 
(1958) observed a similar phenomenon; in both investigations the bend 
occurred at about 900g/mm?. The slope of the curve for (111) crystals 
tested at room temperature gradually decreases after the steep rise to a 
much higher level of stress; the same is observed for (110) crystals, but 
the decrease of slope begins at a lower stress level. 

At 200°C the (100) curve still bends over rather sharply and then 
becomes practically linear with a small slope. The {110) curve rises 
from the start with a gradually decreasing slope. The <111) curve becomes 
almost horizontal in contrast to the curve at room temperature. 

At the higher temperatures the curves for the three orientations become 
horizontal after about 10% shear strain ; at 600°c the levels of stress become 
almost identical for all three orientations. The rate of work hardening 
decreases markedly above room temperature in all three cases, becoming 
essentially zero at 400°c and above. 

Surface examination of the three types of crystals strained at the higher 
temperatures reveals interesting features of the deformation pattern. 
At 400°c, the markings (figs. 4, 5 and 6, Pl. 9) on all three crystals are 
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similar, showing evidence of multiple slip ; the slip lines are closely spaced. 
At 600°C (figs. 7,8 and 9, Pl. 9) the surfaces are strikingly different. 
The (110) crystal shows widely spaced, sharp-edged slip bands, similar to 
the bands that occur on centre orientation crystals deformed at 600°C; 
this is strong evidence of single slip. The cross section of the crystal was 
elliptical after deformation. 

The (100) crystal also shows widely spaced slip bands with these 
differences: the edges of the bands are rounded, and in different regions of 
the crystal, deformation took place on different {111} planes. In any one 
region of the crystal the cross section is elliptical. The rounded edges are 
probably due to a combination of slip on many adjacent planes in the 
neighbourhood of the slip band, and to cross slip. 


Fig. 3 
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The (111) crystal remained circular in cross section after deformation 
at all temperatures. The surface displayed many short, intersecting slip 
lines and deformation appeared to occur fairly homogeneously over the 
entire length of the crystal. 


3.2. Comparison of (111) Crystal and Polycrystal 


A comparison between the tensile stress-tensile strain curves of the 
actual polycrystal and the one derived from the deformation of the CITI 
single crystal was made using Taylor’s criterion (see Introduction). 
Figures 10, 11, 12 and 13 show the results of this comparison. A measure 
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of the similarity, that is, the ratio of the stresses at a given strain, 
Spoly/% 411), 18 given in the table. The strain rate is 4x 10-5 sec—. Also 
included in the figures are stress-strain curves of polycrystals deformed at 
a strain rate of 4x 10-3sec1. The ratio (Srast/Fslow) for polycrystals at 
15% strain is also given in the table. 


Temperature 27°C 200°C 400°C 600°C 
Spoly/F<111) 1-00 (e=11%) 1:16 1-46 
(Gfast/ slow) Tpoly 1-11 1-63 2-70 


tT Ofast—at 4 x 10-3 per sec. ogiow—at 4 x 10-5 per sec. 


§ 4. Discussion 
4.1. Orientation Dependence of Single Crystal Deformation 


Ifa uniaxial stress is applied along a (100) or (111) direction (‘ polyslip ’+ 
orientations) or along a (110) direction (‘multislip’} orientation) of a 
f.c.c. single crystal, the resolved shear stresses are equal on eight, six and 
four slip systems respectively, and are nil on the rest of the twelve available 
systems. In these cases plastic deformation occurs, from the onset, by 
slip on at least several of the equally stressed systems. Furthermore, if 
slip occurred on all the equally stressed systems, in the (111) orientation 
three sets of Lomer—Cottrell dislocations would form along stressed (110) 
directions, in the (100) orientation two such sets would form along un- 
stressed (110) directions and in the (110) orientation only one set would 
form along an unstressed (110) direction. 

The differences in the stress-strain relationships for these three orien- 
tations, found at room temperaure, can be rationalized in terms of the 
above differences in the numbers and distribution of Lomer—Cottrell 
dislocations (Hosford et al. 1960). 

The results of the present experiments show that as the temperature of 
testing is increased, the shear-stress shear strain curves for the three axial 
orientations ({111), <100) and <110)) become more and more similar, 
until at 600°c they are identical. On the other hand, the surface appear- 
ance of the crystals strained at 600°c indicates that the mode of deformation 
depends drastically on orientation. It should also be pointed out that for 
deformation at 600°c, the shear stress-shear strain curves of crystals 
oriented for single slip (‘centre’ orientations) are the same as for the three 
‘corner’ orientations discussed above. One is thus led to conclude from 
the similarity of the shear stress-shear strain curves for all these cases, 
that the mechanism of strain hardening at high temperature is the same 


+ We are defining these terms as follows: polyslip orientation—one having 
five or more slip systems with equal resolved shear stresses; multislip orienta- 
tion—one having more than one but less than five slip systems with equal 


resolved shear stresses. 
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regardless of orientation. This happens in spite of the fact that the 
numbers of Lomer—Cottrell dislocation sets expected for form are not the 
same for each orientation, and that the number of both active and equally 
favoured slip systems is also different in each case. It is apparent, there- 
fore, that in aluminium at 600°C Lomer—Cottrell barriers become in- 
effective as obstacles to slip. 


4.2. The Relation Between Single Crystal and Polycrystal Deformation 


To satisfy continuity conditions at the grain boundaries, a grain in a 
polycrystal deforms plastically by slip on several slip systems (five systems 
in the most general case). Hardening by the formation of Lomer— 
Cottrell barriers will also take place within a grain of a polycrystal, the 
only difference between the single crystal and the polycrystal being the 
presence of grain boundaries in the latter. In the case of aluminium 
deformed at room temperature, it was shown by work on multicrystals 
(Elbaum 1960) that the only effect of the grain boundaries is to force addi- 
tional slip systems into operation in order to satisfy the continuity con- 
ditions; the boundary contributes no intrinsic strength of its own to the 
polycrystal. 

The results of the present study on polyecrystal deformation show 
that, at room temperature, the stress-strain curve of a polycrystal can be 
derived from the deformation of a (111) single crystal using Taylor’s 
criterion. Equally good agreement is obtained in the case of (100) crystals 
for the first few percent of plastic strain. For higher strains the <100) 
crystals undergo a characteristic change from deformation by slip on 
‘many’ slip systems to deformation consistent with slip on two systems 
only. ‘This difference between the two orientations is thought to originate 
in an inherent instability of the (100) orientation. In particular, in the 
«100 case, it can be shown (Hauser and Jackson 1960, Kocks 1959) that if 
deformation occurs to any great extent on a particular slip plane, it will 
continue on that plane with a resulting gradual change in orientation. 
For the specific purpose of the present discussion it remains significant 
that good agreement is found between the behaviour of polycrystals and of 
(111) or (100) single crystals, as long as the latter deform by slip on several 
systems. ‘These results suggest that the energies consumed in the defor- 
mation of a polycrystal, or of a single crystal oriented for polyslip, consist 
in both cases of the energy for polyslip and are essentially the same. More 
specifically, these results indicate that at room temperature the work 
hardening mechanisms in polycrystals and in single crystals oriented for 
polyslip are also the same. It will be assumed that work hardening is due 
to the formation of dislocation pile-ups caused by Lomer—Cottrell barriers. 
This mechanism is not sensitive to the presence of grain boundaries in a 
polycrystal because the spacing between effective barriers to slip in each 
grain is small compared to the size of the grain. 

At elevated temperatures the stress-strain curves for polycrystals, 
derived from the deformation of (111) single crystals, no longer coincide 
with the actual curves for polycrystals. It should also be pointed out 
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that at 600°C, after the initial deformation of 2 or 3%, no work hardening 
is detected in either single crystals or polycrystals (deformation takes place 
under a constant stress). This implies that a steady state is reached, 
where the dislocation configuration remains the same in the crystal during 
deformation. In particular, the number of annihilated active dislocations 
is equal to the number created; however, the steady state is reached at 
different stress values in the single crystal than in the polycrystals. Also, 
both single crystals and polycrystals are strain rate sensitive, but at 600°C 
the strain rate sensitivity is higher for polycrystals than for single crystals. 
These observations indicate that the predominant relief mechanisms are 
different in single crystals and in polycrystals, although both mechanisms 
lead to a steady state dislocation distribution. Of the two mechanisms, 
at least one, or both, must depend on the nature of the obstacles to slip. 
If both mechanisms were independent of the type of obstacle, or if the 
obstacles were the same in the two cases, one would not expect to find any 
difference between single crystal and polycrystal deformation. (This 
happens at room temperature where Lomer-—Cottrell dislocations are 
thought to be responsible for work-hardening in both cases.) Since no 
orientation dependence was found in the stress-strain relationship of 
single crystals deformed at 600°C, it was postulated that Lomer—Cottrell 
dislocations become ineffective as obstacles to slip at this temperature. 
The only likely obstacles to slip that remain in a polycrystal and which 
are not present in a single crystal are grain boundaries. It is therefore 
suggested that at 600°c in single crystals and in the interior of a grain in a 
polyerystal, break-up of Lomer—Cottrell dislocations occurs, whereas a 
grain boundary in a polycrystal is capable of supporting larger dislocation 
pile-ups. Relief from these pile-ups at a grain boundary may occur 
through climb, cross slip, or escape of dislocations into the grain boundary. 
No clear distinction can be made between these relief mechanisms on the 
basis of the available results. 

It should be pointed out that grain boundary sliding is not important at 
600°C even at the slowest strain rates used in these experiments. If grain 
boundary sliding required a lower stress than slip within the grains, the 
stress-strain curve for a polycrystal would be expected to lie below the 
stress-strain curve derived from (111) single crystal deformation, whereas 
the opposite is found to occur. 


4.3. The Validity of the Taylor Criterion 


Using Taylor’s criterion 


the stress-strain curve for a polycrystal can be derived from the stress— 
strain curve of a single crystal only when the microscopic deformation 
mechanisms, and particularly the work hardening mechanisms, are similar 
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in both cases. This condition is satisfied at room temperature, if the com- 
parison is carried out between polycrystals and single crystals deforming 
by polyslip. Here work hardening is in both cases due to obstacles to slip 
(Lomer—Cottrell dislocations) formed by similar dislocation interactions 
within the single crystals or within the grains of a polycrystal. This 
mechanism is not sensitive to the presence of grain boundaries, because 
the spacing between effective barriers to slip in each grain is small com- 
pared to the size of the grain (or spacing between grain boundaries). 

At elevated temperatures, where grain boundaries become more effective 
barriers to slip than Lomer—Cottrell dislocations, polyslip deformation in 
single crystals and in polycrystals is not the same. In this case the stress— 
strain curve for a polycrystal can no longer be derived from the stress- 
strain curve of a single crystal, using Taylor’s criterion. 


§ 5. CONCLUSIONS 


It may be concluded from the experiments described above that: (1) 
Taylor’s criterion for deducing the stress-strain behaviour of a polycrystal 
from that of a single crystal ceases to be valid (in the case of aluminium) as 
the temperature rises much above room temperature, and (2) the principal 
work hardening mechanism that operates at lower temperatures (namely, 
the formation of Lomer—Cottrell dislocations that become obstacles 
against which other dislocations can pile up) ceases to function at elevated 
temperatures. That is, the Lomer—Cottrell locks are no longer effective 
obstacles to slip at higher temperatures in both single crystals and poly- 
crystals. The only barriers that remain effective are grain boundaries 
present in polycrystals. 

The experiments dealt with in this investigation and in previous studies 
of room temperature plastic deformation of multicrystals (Hosford e¢ al. 
1960) were carried out on aluminium only. It is, therefore, not possible 
to extend all the conclusions of these investigations to other metals with 
the f.c.c. structure. Similar investigations in this area, on other metals, 
particularly on f.c.c. metals that are characterized by a low stacking- 
fault energy (for example copper), are necessary for a more complete 
understanding of the relations between single and polycrystal deformation. 
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ABSTRACT 


Cross-slip now appears to be an important factor in the initiation of fatigue 
cracks, and this lends support to the Mott screw-cycling model. The 
reason for progressive unidirectional cycling of the required kind under 
alternating stress is, however, not easily understandable. Two possible 
gating mechanisms are here proposed to account for this effect, both of 
which involve interactions with Lomer—Cottrell barriers. One of these 
mechanisms can provide an explanation for the very rapid appearance of 
extrusions in a few cycles of stress, and does not require any initial internal 
void to exist, as does Mott’s original model. 


§ 1. INTRODUCTION 


Or the mechanisms proposed to account for the phenomenon of fatigue- 
induced intrusion and extrusion in slip bands, that which involves the cross- 
slip of screw dislocations (Mott 1958) appears to accord most satisfactorily 
with the experimental observations. The double-slip system of Cottrell 
and Hull (1957) and the mechanism of Thompson (1959) are both open to 
certain objections; the former in requiring intrusions and extrusions to 
be necessarily associated, and to occur in different slip bands, and the 
latter in its limitation to the case of extrusions only.. The experiments of 
McEvily and Machlin (1959) on certain single crystals (LiF, AgCl, NaCl) 
were designed specifically to test some of the features of the different 
theories, and the general outcome was that extrusions and intrusions, and 
susceptibility to fatigue failure, are strongly associated with ease and pre- 
dominance of cross-slip. McEvily and Machlin go on to examine possible 
dislocation configurations which could lead to the occurrence of intrusion 
and extrusion in the same slip band under an alternating stress. 

One of the difficulties of the cross-slip hypothesis is that systematic uni- 
directional cycling of screw dislocations is required, of such a type that 
when the stress is reversed the screw dislocations do not simply retrace 
their forward paths; in other words, the dislocation movement is not 
reversible. May (1960) takes the view that no specific mechanism need 
be invoked to account for the observations, which could arise as the con- 
sequence of a kind of random cross-slipping of dislocations from one slip 
plane to another. This model is developed statistically by May, who 
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assumes that as the surface roughens, the valleys so created become pre- 
ferred regions for slip, and that these deepen until a crack may be considered. 
to exist. With certain assumptions, for example that groups of about 20 
dislocations are the operative elements, the model leads to values for V 
(cycles to failure) of the right order of magnitude. The model does not 
attempt to examine how the random cross-slipping might be initiated, 
nor can it explain how relatively large extrusions, up to ~5y in length, 
can develop in a relatively few cycles. 

’ In view of the support that the cross-slip hypothesis of Mott has received 
from recent experiment it appears worthwhile to attempt some closer 
examination of the possible gating mechanisms which might lead to 
progressive cycling of the required type. 


Fig. 1 
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The impeding action of stacking faults associated with Lomer—Cottrell dis- 
locations, showing in the sequence I, II, III the operations required for 
the extrusion or intrusion of the shaded area. 


§ 2. Gatine By Srackine Faunts associarep with LomER—CorrrRELL 
DISLOCATIONS 


Consider the situation depicted in fig. 1. Screw dislocations approaching 
the Lomer—Cottrell (L-C) dislocation could be impeded by the stacking 
fault OPP’O’. If the sign of the partial PP’ is such that it moves 
towards the sessile dislocation OO’ during this stress cycle, then the screw 
dislocations could pass freely over the stacking fault. On the negative 
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half cycle, OPP’O’ extends, and the returning dislocations are held up by 
the stacking fault. Cross-slip is thus initiated, the screw dislocations 
slipping onto some other parallel plane, as shown. If during their slip on 
this plane they negotiate an L-—C barrier in the open state, then the 
elements exist for some kind of cycling, because in the next positive half 
cycle the dislocations traverse the path indicated in fig. 1. It is not 


Fig. 2 
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Dislocation system required to promote cross-slip by the cyclic blocking of 
screw dislocations by L-C barriers. 


essential for such an ideal two-barrier unit to form in order to introduce 
a bias into the behaviour of dislocations. The principle could be applied 
to a random distribution of such barriers, each operating as a directional 
device. In such a case, the barriers might well constitute the physical 
elements required by May’s statistical theory. One weakness of such a 
model is that the barrier offered by a single stacking fault would appear 
to be too slight to promote extensive cross-slip: this difficulty might Be 
resolved by considering the gating elements as closely spaced groups o 
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similar L—C barriers. However, it does seem more reasonable to seek a 
mechanism which utilizes the blocking action of the L—C barrier itself, 
and this leads to the next model. 


§ 3. Gatine By LomeR—CoTTRELL DISLOCATIONS 


Consider the situation shown in fig. 2. The stressing system is such that 
screw dislocations are generated on ABCD and edge dislocations generated 
on EFGH. In the positive half cycle screw dislocations (say positive) 
traverse the plane ABCD from left to right. Edge dislocations traverse 
the EFGH plane and shift the L-C barrier away from ABCD. In the 
negative half cycle of stress, edge dislocations of opposite signs are generated 


Fig. 3 
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The consequence of the operation of the system indicated in fig. 2. 


the L-C barrier is translated into the ABCD plane, and thus screw dis- 
locations (negative) are held up, and are forced to cross-slip as shown. 
The essential elements for progressive extrusion or intrusion now exist. 
Similar events at other points on these slip planes could result in segments 
of extrusion and intrusion occurring in the same slip band. Very rapid 
development of these features can occur: in fact in a single stress cycle, 
hundreds or thousands of screw dislocations could be generated and 
large extrusions/intrusions result. It is thus not essential for a given group 
of dislocations to be taken around a closed path (although this could indeed 
occur by such a mechanism) but merely necessary that there should exist 
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a path discrimination between positive and negative dislocations in the 
way indicated. Such a model suggests that finite extrusions might well 
be seen after a very few cycles. 

The model does not require any initial internal voids, as does the un- 
modified Mott model, nor does extrusion imply that such voids must be 
formed. The passage of the screw dislocations will leave a concentration 
of positive and negative edge dislocations (see fig. 3) and this is a sufficient 
consequence of extrusion or intrusion. However, high dislocation con- 
centrations could lead to an actual void forming, in precisely the manner 
which Mott (1958) envisaged. This modified model therefore accounts 
for both of Mott’s original features by a single basic mechanism. 
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ABSTRACT 


It is proposed that the change in crystal structure on alloying the noble 
metals with elements of the B sub-groups (e.g. copper—zinc) is due to the 
Fermi surface touching the cube faces of the Brillouin zone of the face- 
centered cubic structure rather than the octahedral faces which are already 
touched in the pure metals. A reasonable band structure for copper is con- 
sistant with this suggestion. 


It has been established empirically (Hume-Rothery et al. 1934) that, 
when the atomic size factors are favourable, the maximum solid solubility 
of many B sub-group elements in copper and silver occurs at an electron 
concentration of approximately 1:4. The explanation of this advanced 
by Jones (1936) was that, in so far as electronic energy is concerned, if 
the density of states or N(E) curve is of the general form of fig. 1, there 


Fig. 1 


NCE) 


will be a tendency for a phase to become unstable in the range A-B. For, 
if all other possible crystal structures are considered, there will probably 
be one or more of them for which the N(E) curve continues to rise after 
the electron concentration of A, and which will, therefore, accommodate 
the electrons with a lower energy. rae 

For the face-centred cubic structure, the first Brillouin zone is the 
well-known truncated octahedron. If the Fermi surface is assumed to be 
spherical, it touches the octahedral faces of the zone at an electron 
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concentration of 1-36, in remarkable agreement with the observed values 
of the maximum solubilities. In later work Jones (1937) calculated the 
N(E) curve for the face-centred cubic structure more accurately by 
assuming the energy gap across the octahedral faces of the zone to be that 
indicated by the optical properties of copper. The peak at A then 
corresponded to an electron concentration of about 1-0, and thus indicated 
a marked distortion of the Fermi surface. This was later confirmed by 
the work of Pippard (1957) on the anomalous skin effect in copper which 
indicated that the Fermi surface has just touched the octahedral face 
of the zone, and is in agreement with more recent work by Shoenberg 
(1960) on the de Haas—van Alphen effect which shows that the zone face 
is touched by the Fermi surface in copper, silver and gold. 

In the calculations of Jones referred to above it was assumed that the 
same energy gap would exist across the dodecahedral faces of the Brillouin 
zone for the body-centred cubic structure. Marked distortion of the Fermi 
surface was again found, and the energy relations were such that in the 
equilibrium of the a/f8 brass type (between face-centred and body-centred 
cubic structures), the «-solubility limit was again in general agreement 
with experiment. 

The calculations of Jones were thus satisfactory in predicting the 
marked distortion of the Fermi surface, and in indicating the approximate 
position of the «/«+ 8 solubility limits, but their acceptance meant the 
dropping of the correspondence between the position of a peak on an N(E) 
curve, and the electron concentration at the limit of a solid solution, 
which had been such a remarkable feature of the early theory. Cohen and 
Heine (1958) pointed out that with solutes such as zinc, the s valency 
levels were lower than the p valency levels by an amount which was 
considerably greater than for copper, and for this metal they concluded 
that, at the centre of the octahedral face, states inside the zone would be 
purely p-like, and those outside purely s-like. They, therefore, argued 
that the solution of zine in copper would tend to lower s states relatively 
to p states, and would thus decrease the energy gap at the zone face, 
and make the Fermi surface more spherical, so that with increasing 
zinc-content the surface would at first shrink away from the zone-face, 
and then swell out again in a more spherical form and so satisfy the 
conditions of the original simple theory. For silver, Cohen and Heine 
concluded that the Fermi surface was nearly spherical, so that the 
conditions of the simple theory applied. For gold they concluded that, 
at the centres of the octahedral faces, the states inside the zone were 
purely s-like whilst those outside were p-like. They thus argued that 
solutions of zinc, and similar elements would increase the energy gap and 
make the Fermi surface more distorted, so that contact would occur at 
a lower electron concentration, and in this way they explained the lower 
solubilities of elements in gold than in silver and copper. 

For copper, the explanation of Cohen and Heine can still apply, but 
the recent discoveries that the Fermi surface has touched the zone in silver 
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and gold disproves their interpretation of solid solutions in these elements, 
and the following alternative explanation is suggested. 

The complete N(E) curve for the face-centred cubic structure is of the 
form of fig. 2, when the peaks at A and B correspond to contact of the 
Fermi surface with the octahedral and cube faces of the zone. For a 
spherical Fermi surface the peak A corresponds to an electron concentration 
of 1-36 and if the surface continues to expand whilst retaining a spherical 
shape on its free surfaces, the peak B corresponds to an electron 
concentration of 1-88. If, therefore, in the actual structure of copper 
the peak A has been shifted back from a value of 1-36 to one of less than 1-0, 
it does not seem unreasonable to suggest that the peak B has been moved 
back from a value of 1-88 to ca. 1-4. After the point B, two sets of zone 
faces are reducing the number of available energy states and we may, 
therefore, expect the N(E) curve to fall more steeply after B than after A, 
and so to intensity the effect proposed by Jones. We suggest, therefore, 
that the solubility limit in copper alloys corresponds with the fall in the 
N(E) curve after the second, rather than after the first peak. 


Fig. 2 


a 


Unfortunately, there is no complete calculation of the band structure 
of copper from which one could easily verify this suggestion. However, 
Garcia Moliner (1958) showed that the surface obtained by Pippard could 
be adequately described by the Slater and Koster (1954) interpolation 


formula 
1 1 
«=C,(3—cos 4ak, cos $ak, — cos 40k, cos 34k, — COs tak, cos ak.) 
+ C,(3 — cos ak, — cos ak, — cos aks) 
+ C,(3 — cos ak, cos aks — Cos ak, cos ak, — cos ak, cos ak) 


+ 0,(3 —cos 2ak, — cos 2ak, — cos 2aks) 


Ai ay ekg eerie 


with ¢/C,=3-6301, C,/C,;=0-0995 (C, and higher coefficients zero). 
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The member of this family of surfaces which just touches the (200) face 
of the zone has <¢/C,;=4. Sampling in 27 directions (representing 794 points 
by symmetry) we find that this surface contains 1-63 electrons per atom. 
This is rather higher than the observed ratio at the phase boundary. 

The general form of energy contour at the face of a Brillouin zone is 
such that the Fermi surface will probably be pulled out as it approaches 
the (200) zone face. This is confirmed by Shoenberg’s work which suggests 
that the Fermi surface of pure copper is pulled out in the (200) directions 
(as the area of a central cross section in the (200) direction exceeds that 
in a (111) direction by 2-7%). A surface described by (1) with 

e/C, = 3°84, C,/C,=0-119; C,/C, = 0-064 eee 2) 
contains one electron per atom and has the correct area of contact on the 
(111) faces and the correct ratio of cross-sectional areas (though these 
are taken 8°% below Shoenberg’s areas which exceed the area of a central 
section of a free electron sphere). The member of the family with (e/C,=4) 
which touches the (200) face contains 1-21 electrons per atom. This is 
rather too small but it is probable that the true Fermi Surface is not as 
depressed in the (110) directions as (2) suggests. The addition of suitable 
harmonics to allow for this might well give a surface touching the cube 
faces at a higher electron concentration. For example 

e/C,=3'823, C,/C,=0-118, C,/C,=0-045, C,/C,=0-01 

fits Shoenberg’s surface as well as (2) and the member of the family in 
contact contains 1-25 electrons per atom. These calculations show that 
the position of the (200) peak on the N(E) curve is very sensitive to the 
exact form of the Fermi surface. The calculations are tentative and 
qualitative, but suggest, that the identification of the a-solubility limit 
with the fall in the density of states after the (200) peak is at least 
reasonable. 

The a/a+f type of equilibrium is found only in alloys of copper and 
silver. In the gold alloys, the more electro-negative nature of gold results 
in an increased electrochemical factor, and in some systems (e.g. Au—Cd, 
Au-Zn) the random body-centred cubic 8-phase is replaced by an ordered 
Cs—Cl type of structure with a maximum freezing point in the equilibrium 
diagram. In these systems the smaller solubility of Cd and Zn in gold 
follows from the usual free energy principle of the effect of a stable 
intermediate phase on the limits of a primary solid solution. Heine 
(private communication) suggests that, since in the free atom of gold 
the s level is so much lower than the p level, it is still possible that the 
mechanism of Cohen and Heine is responsible for the lower solubility 
limit in the gold alloys. The solubility limits of the face-centred cubic 


structure in equilibrium with a close-packed hexagonal structure are not 
yet understood. 
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ABSTRACT 


A linear integral equation is derived for the correlation function in the 
presence of a magnetic field, and solved and the errors estimated. 


§ 1. LyTRopUCcTION 


In a fully ionized plasma whose charge density in p(x, ¢), a fundamental 
quantity in the description of the plasma in say thermal equilibrium is 
the correlation function p(x, t)p(x’,t’)), where the brackets mean that 
the quantity is averaged over the thermal equilibrium state. This quantity 
is a directly observable quantity by probe measurements, and plays a 
key position in the theory of transport phenomena of states slightly removed 
from equilibrium since it appears in the kernel of the transport equation. 
It will be studied in this paper in a simple approximation which will allow 
the inclusion of an external magnetic field. A linear integral equation 
for the correlation function will be deduced and solved and the corrections 
briefly discussed. A more general approach which can lead to an analysis 
of the complete solution will be discussed in a later paper since it requires 
a much more elaborate mathematical apparatus. ‘The application to 
transport theory is also deferred to a subsequent paper. 

To simplify the analysis the gas will be considered.-to consist of only one 
kind of particle with a neutralizing background. Before considering the 
whole problem some simplified versions which can be solved exactly will 
be considered. If x, are the coordinates of the electrons, the density is 


o(x)= > 8(x—x,)-N/V Seiad al) 
* which has Fourier components 
pe= > exp (tk ory). Hos agtck. Wray aaa 1; 2) 


The overall neutrality has the effect of making py zero, so that the equation 
of motion of an electron becomes 


G. 
ee ee cap iiKe(x. —3G))t = X,~ By) 4) 2 (18) 
. V io e 


+ Communicated by the Author. 
+ Vacation Consultant. 


62 S. F. Edwards on the 


where k runs over the values (27/R)(1, %2, 23), & being the side of the box, 
V = R%, and nj, no, nz are integers. 
The correlation function C(k,t) is defined by 


Ok, t) = (pg(t)p*(0)) ms iss 
so that C(k,t) is the Fourier transform with respect to x—x’ of 
(p(x, t)p(x’, t’)>. MS) fi RS | 


Consider firstly the case of free particles. Then 


Mook 1.6 
x, (t) =X, +V,E. \ eee 


Averaging against the Maxwell distribution 


Sa 0)) 
x | ex| - 2 dmv 2/KT +ik.v,t+ik(x, —X,) ] [[ dv, dx, 


=5 ap 
(1.7) 
where 
C= [exp (— > $ mv2/KT) [[ dv, dx,. ~ Se ERS) 
Therefore 
C(k, t)=N exp [— (kT )t?/2m]. wee es he) 


Now consider the case with Coulomb interactions and long wavelength, 
i.e. small k. This means that the first term, that of order k? in 


Pr= >, —(k.x,)? exp (tk.x,) +7 > (k.x,) exp (tk. x,) (1.10) 
can be ignored so that 


_ 4ne 
V 


pel ats 5 p(i(k-+j)x,—ijx,). - . (1.1) 


When this is averaged in nee equilibrium only the term k+j=0 
survives, giving NV, so that 


Pull +( Fo) Pe Ojoig ie aes ace ee i, 


Vm 
and 


(Pr(t)pr* (9) ) = (py (O)p,*(0) ) Cos cot - +. (1.18) 


where w is the plasma frequency (47Ne?/mV)"2. 
To evaluate (p,(0)p,*(0) ) allowing for the Coulomb interaction requires 
a solution of the equilibrium properties of the plasma which will be quoted 
here. In first approximation 
N 


Ojo." (0) \ =e 
(Px Pr ( )> 1+ (k2A2)-2 
where A=(VkT'/47Ne?)"?, the Debye-Huckel length. (Higher approxima- 


tions to oe can be obtained for example by the methods given by Edwards 
(1959).) 


(1.14) 
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Next consider the effect of the external magnetic field on these two 
limiting cases. For the free particle the equation of motion is 


MX, = (e/c)x, x B. ee. ea By) 


Taking B in the Z direction and Q=e|B|/mc the solution is 


v(t) =%, + O4(u, sin Qt+v,(1—cos Nt), v1 (1-16) 
y At) =y, + Q41(v, sin QF +u,(1— cos ME)), een LAL) 
z,(t)=z,+w,t (ee ee, eee eee te (11S) 


Proceeding as in (1.7)—(1.9), one obtains 
(Px(t)py (0) >) =N exp [—k,?«T 2m 
—(k,? +k,?)(1 — cos Qt)KcT/mQ?] 5 Te) 
or in invariant form 
N exp [ —(k. B)?#«7'/2m.B? — (k x B)?(1— cos Ot)«T'/OQ2B2m]. (1.20) 


For the interacting case, neglecting as before terms in k?xT'/m compared 
with w? and now also 92, the simplest way to close the equations is to go 
to fourth order 


ee ces So ye Beep Gk xo) ee ere (de 21) 
ey a ule k. [(k B 
Pe PK Sa [(k x B) x B]p, 
es. k.(((v,x B)xB)xB)  . . . (1.22) 
11C 
Lae = ¥ k. (v, x B) exp (ik. x,) 
+ wQ?[(k2B? —(k. B)?)/(k2B2)]p,;  . . - (1.23) 
i.e. 
2 k2B2—(k. By 
(5 +o?) yh=— (5 +a") 026, + QF? es (1.24) 
1.€. i 
Ue 2) p= Riera 25 
a 
where 


y= cos (sin (k. B)/kB). 
This equation corresponds to oscillations of frequency q where 
(q? — w?)(q? — Q?) — pQ?w? = 90 ; see ages (1.26) 


= @=Hor+ OF + V/[lo®— 232+ 402904). a3) 


0, giving the rocts g?= 
The two simple cases are k and B parallel, i.e. w= 
and g* = 0?, ne k and B perpendicular giving q?=w?+ Q? or0. As p ae 
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from 0 to 1, the solution w goes to »/(w?+ Q?) and the eet Q goes to 
zero. In peneral there will be two solutions +w(k) and +¢(k) say, and 
the analogue of 

Cog(t)puc*(0)) = <px(0)pye*(0) ) C08 wt 
becomes 


Smet =“ cost + ; 008 pt | (1.28) 


The equilibrium distribution depends upon exp ( - rs. and since the 
kinetic energy is }(p—eA)*, ie. still }4mv*, the magnetic field has no 
effect upon <p,(0)p,*(0)) which is the same as in (1.13). 

In general the correlation function C(k,t) will have collective and 
individual particle motion aspects and in the next section integral equations 
will be derived containing the two extreme cases as limits. 


§ 2. Tue InrEGRAL EQUATION FOR THE CORRELATION FUNCTION 


Consider at first the case when the magnetic field is absent. Then the 
equation of motion 


mi (82 Sy es 


can be formally integrated into 


e* r7 . X,(s)—x,(s) 
x (=x, +Vv,f+ — I d ds : 2.2 
m Jom Joe xe) —ag0)P a8 
X(t) = Xp Vib a(t) ceay, a 9 a ee ee ee ee 


where x,, Vv, will be used to mean x,(0), v,(0). The correlation function 
O(k, t) is ne 


= 2 Serpe x,)]> so. feat So aed 
“Shasiact «2a « (2.5) 


Now assume that the motion of any one particle deviates slightly from the 
rectilinear path, and that only the collective effects of the many slight 
deviations matter, so that exp [ik.a,(t)] can be expanded in a series. 
A more general disoussion of this expansion will be given in $4. This gives 


C(k,t)= S <exp [ik(x,—x,) +ik.v,t]) 


ay 


47re (k. j) 
+ —— ds eth 
Ve > bor sf rf’ ce 


x exp [ik(x,—%, +¥,t) +i]0%,()—X,(8)) (2.6) 
In the same spirit the x,(s) can also be expanded giving 
C(k,t)= » (exp [7k(x, —x,) +ik.v,t]) 
i 
ae es >» >) an asx! 


1) nhs 


x (exp [i(k e j)x, +o(ké+ ee — 1jx,(s) + vkx,]). (2.7) 
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Now in equilibrium the average value of >. exp [i(k+j)x,] is zero 
unless k+j is zero, when it gives VN. Assuming that this averaging will 
not interfere with the averaging of exp [—ijx,(s)+ikx,] the equation 
becomes, putting g(k,t)C(k, 0) for the first term on the right 

rt il 
O(k, t)= Ck, 0) 9(k, t)—w | i | Hon sea\C las) (2.8) 
0 0 
or integrating by parts: 
rt 
O(k, t) = C(k, 0)g(k, t) — w? | (Gg (ee rye a) ar. (2.9) 
0 


This equation is readily solved by the Laplace transform: 
If 


Ok,p)=p| erp eCUc dL ee ee (210) 
0 
qk, p)=p | exp (—pt)g(k, t) dt 
Jo 
I 
= sit exp (p?/4a) erf (+) See CLAS 
where «= k?xT'/2m, and 
G(k, p)=p | exp (—pt)ig(k, t) dt 
0 
eee raat (2.12) 
Pp ap p > 
then 
O(k,p) _ Gk, pyClk, 0) _ w® Gh, p)C(k, p) (2.13) 
Pp p Pp 
A Jk, p)O(k, 0) 
Y k =) eS 
mveP) 1+w?G(k, p)/p 
and 


O(k,0) pr exp(pt) kp) 
= S| io, P t+eO(k, pip 
(2.14) 


This readily reproduces the limiting cases w=0 trivially, and for 
x= k2xT/2m=0, Jk, p)=1 and G(k, p) =1/p® so that 
re C(k, 0) ay Ox) (wt) dp : ‘ 2.15 
C(k, t)= [iin pr+o ( ) 


27 
= O(k, 0) cos wt “lia et en Oe 9) 

as in (1.13). 
Though the integral (2.14) cannot be evaluated without further 
approximation, an important use of it is in transport theory where its 


E 
0 P.M. 
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time integral is required in the kernel of the Boltzmann integral equation. 
This can be obtained analytically, for let it be called A,,: 


A= | Ok, t) dt. Ree ee rE 
O 
Then 
Ay= PueE) (2.18) 
Pp p=0 
_ Gil Dy pyClk, > | Ree en) 
1+0G(k,p)/p \p-0 
1+ w?/2a 
an (2m/k2kT)¥2 ae 
fale Oi2 eM 1+ w2m/k2«T ot 


But w2?=47e2N/mV and \?= VkT'/4re?2N, therefore 


A= de (Se) a ae. Sod RBS) 


The dynamical effect of the Coulomb interaction upon A, is to multiply 
it by a factor (1+(k?A?)-1)-4, so using the fact that C(k, 0) is N times this 
factor 


1/2 
Ay= Jszlcn) (L+(R22)4)%. |. (2.28) 
ZK\K 


Approximate evaluations of (2.14) can be made near the limits w=0 
and a=0. For the latter one can compare the discussion given by Landau 
(1946), for the problem of plasma oscillations, which, though discussing 
a slightly different problem and using a different derivation, reduces to 
very similar analysis to that given above. 

Another way to approximate to O(k,t) has been given (Edwards 
1958) by fitting the exact Taylor expansion of C(k,t) to the form 
cos +/ (at) exp (— bf”), and this checks with a similar approximate evaluation 
of (2,14). 

The correlation function for the current can be deduced immediately 
for the above in the absence of a magnetic field since the only vector 
available in the problem is k. They are 


Ci (px (0) ) = k(k*) *0C/0t, » + + (2.24) 
vhs k kg @C 
CHa) I eg()) = oe ok (2.25) 


These results can also of course be obtained by setting up integral 
equations for them and solving directly. 
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§ 3. THE CorRELATION FUNCTION IN THE PRESENCE OF AN EXTERNAL 
FIELD 


The presence of a constant external magnetic field (0,0, B) alters the 
equations of motion to 


x(t) —x,(t) an 
MX , me DT wes tatcn = > & od ICH 


which has the formal solution 


@(t) =x, + Qu, sin QE + v,(1 — cos Qt)) 


Ar 
ef se dr {cos (Q(7 — s))(x, (7) —x,(7)) 
+sin Q(7—s)(y,(7 (r))}[x,.(7 i) | et een S32) 


y(t) =y, + Q4(v, sin Qt —u,(1 — cos )) 


ee. _ 2 ie Ie dr {cos [Q(s — 7) ](y,(7) —y,(7)) 


+sin [Q(7—s sale (r) }}lx, —x,(7)|-8, ee os) 
z,(t) =z, +w t+ aes ds pure ee ve (3.4), 


ev 


One proceeds exactly as in the last section, (2.5)-(2.9), the final equation: 
being 


t 8 
C(k, t) = 9(k, t)C(k, 0) — w? | ds | drQ(s—r)g(k,t—7)C(k,t) (3.5), 
0 0 
where g,(t) is the appropriate function in the magnetic field 
exp [—(k. B)??«T/2mB? — (k x B)?(1 — cos Ot)eT/O?B?m] (3.6) 
and 
Q(s—t) = (kB?) [(k. B)? + (k x B)? cos Q(s —7)]. (3.7) 


In this form (3.5) clearly reduces to (2.8) when Q goes to 0. 
By integration by parts (3.5) becomes 


t 
O(k, t) =9(k, t)O(k, 0) —w? | {aan | (k. B)2(t—7) +(k x B) 


x g(k, t —7)C(k, 7) ar » (3.8) 


sin Q(t —7) 
oo | 


This equation can now be solved by the Laplace transform, since 


O(k, p) p) _ Gk, p)C(k, 9) 2 ep) (ke. B)? “an »P) 


= 


p He v 
(k x B)? G(k, p —1Q) e glk, p+) kp+t =) | C k, 3.9 


E2 
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and therefore 


By ls €}+t0 pl 
okt) =) | exp (pt) 


270 €—-10 p 


2 ° 7 k, — 792) 
(3.10) 


x 


the symbol . meaning the imaginary part. 
As before the time integral of C takes a simple form, but not so simple 
that it can be evaluated explicitly 


C(k, 0) | g(k, t) dt 


\ en 
er | Q(t)g(k, t) dt 
0 


(3.11) 


§ 4. HigHER APPROXIMATIONS 


To investigate the approximation used above, before analysing the 
complete problem, consider the next approximation to the integral 
equation; for simplicity work without the magnetic field though there 
is no reason in principle why it should not be included. It turns out 
that to the next order the new terms coming from the expansion of the 
exponential and from the further expansion of w,(f) give 


ne =I I. dr [; doy (i do i doy a a “I 
“3 ak io) ete doy A} Roney) | (4.1) 
where e 0 | 


Rit, oy, 0,)= > <exp [ix,(k + j+1) + tv, (kt + joy + ley) 
aBys . 
—j.x,(o,) —1.x,(o,) —tk.x,]). (4.2) 


These terms can be estimated from the original integral eqn. (2.8) and 
are smaller by a power of the density, except in as far as there will be some 
terms of the expansion in which pairs of the indices «, 8, y, 5 are equal. 
These correspond to the fact that the expansion of the collision of two 
particles as a power series in e? has nothing to do with the density but only 
_ with the smallness of e and the analyticity of the scattering cross section 
in terms of e®. For the Coulomb case this is strictly speaking a difficulty 
for the potential is unbounded near the origin, and one has the same kind 
of difficulty that arises in say the Fokker—Planck equation. This in 
practice is not very serious and can be avoided crudely by cutting the 
potential off at small distances, i.e. large k, or by using the correct 
description of two-body scattering after the manner of the recent 
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developments in nuclear structure (see e.g. Watson 1953). A more general * 
discussion of the higher approximations to the correlation function will 


be given in another paper since they involve more powerful mathematical 
techniques than are involved here. 


ACKNOWLEDGMENTS 
The author would like to thank Professor B. H. Flowers and Dr. W. B. 
Thompson for their discussions of this work and Mr. T. Gaskell for checking 
the calculations. Part of the work was done whilst the author was 
Vacation Consultant at A.E.R.E., in the summer of 1958 and the author 


would like to thank the head of the Theoretical Physics Division for his 
hospitality. 


REFERENCES 


Epwakps, S. F., 1958, Phil. Mag., 3, 302 ; 1959, Lbid., 4, 1171. 
Lanpav, L. D., 1946, J. Phys., Moscow., 10, 25. 
Watson, K. M., 1953, Phys. Rev., 89, 575. 


a 
es 
A eryew aS ° 
: s.r ee Bite 
z f, i tye * Ee 
: ; - 
; fos cue iys he 
= pe hb 
<a 
& <i 
a <i 
iS a ee 
= 
ms pa a 
_ i : 
“~ - — ss" 
: — 
a , 7 
be 
i ~ 7 e a: 
er ro. 
4 , ae S a 


; Phe a eeu 
ake ae to eg ae 
ee 


A New Approach to Transport Problems in F ully 
Ionized Plasmast 


By S. F. Epwarps and J. J. SanpERson 


Theoretical Physics Department, University of Manchester, and 
Theoretical Physics Division, A.E.R.E., Harwell} 


[Received July 11, 1960} 


ABSTRAOT 


Exact formulae for transport coefficients are given in terms of one and two- 
body equilibrium Green functions in velocity space. Equations are derived 
for these Green functions which apply under plasma conditions and which 
automatically include shielding and plasma oscillation effects. The method 
is illustrated by simple examples. 


§ 1. IntTRODUCTION 


THE problem of obtaining the transport coefficients of a fully ionized 
plasma concerns the fact that one cannot describe the physical consequences 
of the long-range coulomb potential in terms of the two-body encounters 
of the usual Boltzmann equation of gas dynamics. The many-body 
effects resulting from the long-range forces alter the static and dynamic 
properties of the gas, producing the shielding of the forces and the ability 
of the plasma to sustain plasma oscillations. The object of this paper is 
to derive a method of calculating transport properties when the long- 
range nature of the forces is properly taken into account. The method 
depends essentially on the fact that the equilibrium properties of the plasma 
are completely given by the Gibbs distribution exp(—H/«7'). If one 
considers situations not far removed from equilibrium one can derive 
transport equations based on the static averages derived from the Gibbs 
function, and on the dynamic correlation functions which are derived from 
the Gibbs function and the equations of motion. In this way the Debye— 
Hiickel shielding and the effects of plasma oscillations should appear 
naturally in the transport equations and do not have to be inserted in the 
ad hoc manner adopted for example, when one uses a shielded potential in 
the Boltzmann or Fokker Planck equations. ; 

Any approach to the transport problem must be approximate in some 
sense to be useful, and it is important therefore to be able to understand 
the nature and magnitude of higher approximations, and a considerable 
part of this paper is involved with this point. The actual formulation of 
the problem is in terms of the Green functions for the propagation of one 
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and two particles which reproduce mathematically a faithful analogue of 
the actual physical processes taking place. In this paper the simplest 
transport coefficient, the electical conductivity, is treated in detail and the 
equations for the appropriate Green functions derived. The method is 
illustrated by applying it to the Lorentz gas, detailed calculations of the 
conductivity and other transport coefficients being deferred to a subsequent 


paper. 
§2. THE Derivation oF Exact FoRMULAE FOR THE TRANSPORT 
COEFFICIENTS 


The simplest transport problem is that of electrical conductivity, in 


which the distribution function f(r,...r,, U,..-. U,, t) may be taken in the 
form 
f=foth OR Qh Yee Oe eee) 
fp being the partition function 
fo=exp (=A (nee. pris UG ey UE i ee 


and f, the change resulting from the presence of the external electric field 
& which may be time dependent. Liouville’s equation governs f: 


7] 0 F(r.—r t) e C] 
a+ Du, “ + 2, SMacnie) . a) i +> me é. = =. (2.3) 

Since fy satisfies this equation but for the term in 6, and 
H= ¥im,u2— > O(r,—r,), dpe Rancmmeniun Aue 
where ® is the coulomb potential e,e,/|r,—r,|, the equation for f,, ignoring 


for the moment the dependence of f, upon 6, is 
7) 0 Fos aee = 
(5; + 2M. ar, + 2 ar . xa) i= 2eatta- Efo/KT’. - (2.5) 


Consider the Green function of the left-hand side defined by 


7] d 1 0 ; 
(5; oe or, 3 Jat, Fog. Far) Seas sass Yao ae 
XUy’,..,U,’,t, t’,)=T]d(r,—r,’)S(u,—u,’)8(t—-t’). 2. (2.6) 


a 


This function formally solves (2.5) giving 


Flt F thy Unt) | vo [OU yey Mayes Bahay 


, , ’ ex , g , , t ¢ 
MUG 30 tae ents 2 ap Ma Oto Uyekwesr tte tee 
x [][d*r’ da’, dt’. Vai OS tie, ee By 
Hence the current 


i(r, ad bee D> u,e,d(r—r,)[]@u,dr, . . . (2.8) 
B 


is, under uniform conditions, with no contribution from To 


t 
es ’ e ; ; ; 
i=| dt [-- [Zeus Su, » Ef T[ dud [ ]d*r"’ ,d3u, 


(2.9) 
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where dashed functions refer to functions of the dashed variables. If the 
lons are considered fixed so that the current is carried by electrons, the 
current contains two terms: 

Ne? 


j= ml +. fuyGuy’. Ef) [[d3r’'T [abu] [d3u’ Tar dt’ 
: N(N = Te" r ’ OL 1 3 3,./ 3 14 
+ |... | uGuy’. &fy TdT] a’ T fT ]a'at’ (2.10) 


This formula can be reduced by introducing the Green functions G, and 
G1, defined by 


ed 


, , 
Guy ws t0)= | ne (Gases Cy seed, Fire on Ul peet SU (gan Ua af) 


x exp [(— H+ # + 3mvy"*)/(«T) | [a*r dr ,'T[u,d3u,” 
a a+B 


(2.11) 
and 


tA i / / 
G,.(uU,, u, : Us »U> Ss t =] | G(r, ee ey Gc Pie ree u,, Ciceu 
XU,,U,,...U,’ t,t’) exp [(—H’+ F +4 mu’,? + dmu’,?)/KT][ ][d3r,d°r,’ 
a 


geet 
Xx glib dieted Us (2.12) 

These Green functions have a clear physical meaning: G, is the prob- 
ability that in equilibrium one finds a particle at a time ¢’ with a velocity 
u’ and that subsequently the same particle has velocity u at ¢. Similarly 
Gy refers to the velocities of two particles. In both functions the spatial 
coordinates have been integrated out. Consider the electric field switched 


on at time ¢=0, and define the conductivity by 
EAD) Tee) a: eles ei ty see (2h } 
B 


then using (2.11, 2.12) one obtains 
2 3/2 pt ; 
ate ( Se ) | dt’ [rasGamg exp (— mu’ ,?/2xT')d*v,'d*v, 
0 


08 OP \OnnT 
area Tales 2 3 rt ; : 
ae Ors “) | dW [rs Cts exp [(—mu,’?—mu,’?)/2«T7"] 
K 277K 0 
x Bu, diu,'dund?u,’ PP ee er Ne 


In the same way other transport coefficients can be calculated if one 
assumes that the system is not far removed from a local equilibrium, so 
that if there is a local density n(r), a local velocity V(r) and a local tem- 
perature 7'(r), the starting distribution function will be 


fo= TIn(ri) exp] (— Semele Vera) 22+ SO(re eM) + F/eL) |, 
z (2.15) 


The temperature is assumed not to vary across the effective Tange of 
the force so that it is immaterial whether 7'(r;,) or 7'(r;) is used in (2.15). 
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If this form is substituted into Liouville’s equation it will no longer satisfy 
it and just as before an /, will be required depending upon the gradients of 
n, V, T, and giving rise to fluxes of matter, momentum and energy and 
hence to the coefficients of viscosity and thermal conductivity. Though 
it is hoped to return to explicit calculations of this type elsewhere they 
will be pursued no further here and the rest of the paper will be devoted to 
deriving G, and G4, for the electrical conductivity problem. 


§3. Tur Equations FoR THE AVERAGE PROPAGATION FUNCTIONS 


It is convenient to perform the averaging in (2.11) upon the unprimed 
rather than the primed variables. If S is introduced by the expression 


S=exp ie mu? | 2+ 2 Dag) [at | 
x Gexp es mau? [2+ > 5) [x] vera (es) 


which satisfies 


() 0 1 0 0 F, /m,u, O« ) 
Es moe eS — +—.(——-+—]]8 
(5 af », ts or, o » Mime" SOs ou Ors : m4 ( KT 3 = 
= [[d(r,—r,’)8(u,—u,’)d(t—-#’), 2 ox daa meas a ese 


then the average propagator G', is given by 


Go [--- [exe] - (x0+ S yu) [P+ F |e |S 


TL @r,.@r dud yer er | 5. eee 
a1 
Before considering the general case, which is quite complicated, the 
averaging process will be illustrated by considering the simpler case in 
which all particles but particle one are very heavy and can be considered 
at rest relative to particle one, though distributed in thermal equilibrium. 
Using @ for S in this case, the equation for @ depends dynamically only 
upon the variables r,, v,: 


0 7) ee 0 mu, 
é Yo or, +2 m (sm ei 


) )@=8tr.—n1 980% - uy")8(¢—"’). 
(3.4) 
The averaging is static, so employing brackets to mean averaging value, 


1 
(Q)= eal sisi [ee ee u,, Uy, ; t, t’, R,, Rs, Ssh5 R,,) 
x exp( — ¥ OR,, R,) / «I —YO(R,, r) / D TTER, 
a, B a a 


(3.5) 


4=|...[exp| - > o/er | eR, oes at G} 


a, B, 1 


where 
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and capital letters R,,...R,, are used for the static particle coordinates. 


To perform the averaging consider the equation for Q written by iteration 
in the form 


ae aN ae : , r F, 0 mu 
G of oe) a8 n )8(uy — uy )8(t—t = las % sake 
eae rd) mu 
7 a -{(— saat " " " 
~ ~ | | m (sa 7 or) Ql r, , U,,uU, ,Z, ¢") 
F,.” fa] mv,” ” ’ ” R i ng ry eae 
x a (503 +t Sin ore ,1,/,U,”, uy’, t”, t’)d3r,"d2u," dt 
(3.7) 
where @, satisfies 
0 r) é , 
eto) a JOO, Sh) ae (3.8) 
i.e. 
Ooo, 4; —u,(t—¢))s(u,—u, OCS) . 9) 9 (3.9) 
where 


O¢—-f)=1, t>t, 
=) ee Gat 
It will now be argued that the average of the last term in (3.7) can in a 
good approximation be replaced by 


F, 7 a OL Fie/a? mu,” Se ee 
~.({—— +— pI (rowel Aap Apa yet ERA 
(EE ] m & si KT ) Go m ‘ae ¥ KT )<@) ry Oy 
» (3.8) 


The average of the linear term in (3.7) vanishes, and by referring to (3.5), 
one can use integration by parts to show that 


ta) F,. uy 
Suita ae = 9) mona wiele(3:9) 
which can be averaged approximately by the same methods: 
A Oe ae wv Fi ne. ut) 
(eS 2 ee pe)= mak ye PU i 0oFia- (503 ae sia ) 
OO Ua siae dt (a. L0) 


So the final approximate form of the equation satisfied by (@) is 
0 1 ne) " u u " i pe 
ae Q) = me [QoFen R,)- du, Qo(ta, fi") Up Ur", tH) 


0 nu," ” "oyun 
"— — + — Bu,"d*r," dt 
x DF(r R,)- (sa Te ig ))<@) Uy ary 
= 6(r,—ry’)5(u, — uy’)d(t—¢’) he hce eS ee (Ss 18) 
That this is a good approximation can only be checked by looking at the 
corrections. If the exact equation for (Q) can be written in the form 


560) + | MKQ)=8(r A (— 41 )HE—H), come (Gis Las) 
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then (3.11) can be regarded as giving the first approximation to M. The 


next approximation is then 


coe hear (Gemeente) ow ze 
cr (e+ 28))-(3F- Lease (4h) 


x (DF. a Qo dF. (a+ sa ur): _ eats 


since the terms involving the correlation of three forces give a small con- 
tribution of order of the square of the ratio of the interparticle distance to 
the Debye length. 


The term : F) mu’ <ase 
(a ee ) Ae 


m m 


is at a time ¢”, but the average has a material value only if ¢” ~¢, in fact if 
t” —t is less than a collision time, because the two F’s are strongly de- 
creasing functions. This being the case one can integrate the term 
u" .d/dr" by parts, just as in (3.9), replacing it by —F.mu"/«k7T. Thus (3.13) 
becomes 


Gr Boar (ae « Boar Baar (d+) 


(3h Losni(2 -s2)e(sr sr QdF. (55, + “-)). 


(3.14) 


This term has the same dependence upon the coordinates of the free 
particle in both parts, and its magnitude therefore depends on the size of 
multiple averages which are schematically in the form 


(SFSFSFDF)— (SFYF)(SFYP). app Sate SOF 


In the same way the next term will depend on a combination of the six, 
four, and two force averages. The validity of (3.11) depends then on the ~ 
cancellation to a good approximation of the two terms of (3.15), since it 
can be shown that either of these terms can give rise to corrections to the 
transport coefficients comparable with their values derived from (3.11). 
Writing out the terms of (3.15) in detail they are 


> <F(ry— R,)F(r;’— R,)) 
and " 
> Fr,- R,)F(r,’— R,)F(r,”— R,)F(r,” — R,)). (3.16) 


a,B,y, 6 
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It is simpler to consider these quantities Fourier transformed 
2 CFi(ry— R,)F(r,/— R,) ) 
Hs 


D4 


as — » k,jpexp (ik or, + ij. ry’) S (exp (—ik. R,—ij. R,)), 
(3.17) 
‘ 2 Fn RJFi(ry — Re)Fn(ry" — R,)F,(r” — Ry)) 
“ ma ime EDP mIn XP (i(k. ry +j.ny! +p.ry"+q. r,”")] 
x os {exp [—i(k. R,+j.R,+p.R,+q. R,)]). (3.18) 


The homogeneity of the averaged system requires that 
2 {exp (—ik.R,—ij. R,))=5 (k+j) ( > exp (—ik.(R,— R,))) 
a, a, B 


=6(k+j)C(k?), say, er, arom Cock) 
and 
> <exp[—i(k.R,+j.R,+p.R,+q. R,)]) 
a, B, y, 6 
=8(k+j+p+aq)d(k,j,p.q). . . (3.20) 
A detailed account of the evaluation of the averages will be given in the 


Appendix; to the accuracy of taking the Debye length X much greater 
than the interparticle distance ,, one obtains 


O(k?) = N/(1+k-A-). 5 Disa Ee (38519 


It is now essential for the cancellation in the expression (3.15) that the 
four-particle average splits to a good approximation into 
d(k+j+p+q)d(k, j, p,q) 
~ 3(k + j)3(p + q)C(k*)C(p*) + d(k + p)oj + @)C(k*)C(q") 
+ 0(k + q)d(j + p)C(k?)O(9°). (3.22) 

That this is indeed the case in a plasma is discussed in the Appendix. 
The term with 6(k+j)6(p+q) produces the desired cancellation, and the 
two remaining terms can be shown to produce only small corrections to 
the transport coefficients, of order (u/A)*. It is important to realize that 
this cancellation does not happen in more general systems such as a normal 
liquid, though it is extensively employed there, being equivalent to the 
superposition approximation. Indeed a fully ionized plasma is perhaps 
the only physical system for which the approximation réally works. 
Readers familiar with quantum field theory or many-body theory, will 
recognize the above discussion as being equivalent to saying that the 
single particle mass operator may be calculated by perturbation theory 
i.e. by the first Feynman diagram. 

Returning now to the general problem (3.1) one can see how (3.11) 
will be generalized in physical terms. Equation (3.11) shows how a 
particle moves on an average in velocity space through an equilibrium 
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distribution of heavy charged particles, the kernel of eqn. (3.11, being 
derived from the mean second approximation to the motion of the »article, 
or in quantum mechanical language from the mean second Born 
approximation. Now in fact the particles of the plasma are moving with 
their random velocities and with the thermal distribution of plasma 
oscillations. So one may expect not just the instantaneous correlation 
function to enter, but the time dependent correlation function. If R,(é) 
is regarded as the coordinate of a particle in the Lagrangian sense, and is 
taken to be a solution of the equations of motion having the initial 
conditions averaged over the Gibbs function exp(— H/T’), one expects 


C(k,t—)= ( ¥ exp (ik. (Ra ~ Ralt’)))) aes} 


to play the part previously taken by C(k?). 
Another way to look at this is to recognize }exp(ik.R,) as the kth 


fourier component P,, of the density ez 
P(r)= > d(r—R,) . nk ee ea 
in terms of which ; 
U(My=< PP ey 1) PS aerate Achaea 
and 
O(k, t—t') = (P,(t)P;*(t’)). jeoee tees 26) 


The electrons move under the influence of forces from protons and 
electrons, the force acting upon it, depending upon the density of charge. 
So if P;, is used as above for the density of protons, and p,, for electrons, 
then p;—P;, is the density of negative charge. One would thus expect 
in the general problem to find the charge correlation function, 


O(k, tt) = ((pxlt)— Py O)(o.*(t') — Pi) (3.27) 
playing the part of C(k?). 

To derive the differential equation for G, as defined by (3.1) it is 
convenient to think of all the particles except particle one in the 
Lagrangian description, retaining the Liouville description of particle one. 
The equation for S then becomes, retaining capital R, for the Lagrangian 
coordinate of the ath particle, 


7) 0 1 0 mu ’ 
Gi +u,. or, +> m F(r, — R,) : & + —)) S(ry, ry, Uy, U,, R,(¢), t, t’) 


= 6(¢—?')8(r, —r,')8(u, — u,’) Deter Wa wei ie 
where 
mR= 2 F(R, (¢) — Ra(t)) + F(R,(¢)— Ryd), =. =. (3.29) 
and ~ 
R,(¢) = | rS(r1,0,', Uy, uy’, R, (i), t,t’) dr d8u,. . . (3.30) 


The boundary conditions upon the R,’s are 


R,(’)=r,’ and R,(’)=u,’. he (3.31) 
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This system is of course identical with (3.1, 2). Physically particle one 
moves in the fluctuating field of all the other particles according to 
eqn. (3 29). The other particles interact amongst themselves and are 
influenced by particle one through the coupling of (3.29) with (3.30). 
Hence in (3.28) the field acting on particle one depends upon its own motion 
in the past through the dependence of the R,’s upon R, of (3.29). This 
must be a very small effect, however, only of order (1 (N ), and will be 
ignored, so that the system (3.28) and the modified equations 


m,R,( =A R,(t)) oe (EES 


will be assumed henceforth. Now one can proceed precisely as before, 
leading to 


im p> [(Fees R,(¢). ai ry, Uy, Uy", 6,0")  F(ry"— R,(¢’)) 


. ( a] 6 mu, )) Gita,” a” pe t') du,” dr,” dt” = §(u,— Uo =), 
1 
(3.32) 


ea the correlation function C(k,t) defined in (3.27) this can be written 


” 0 mu,” 
+E (a) & C(k,t—t")k. 5a; Colts ti" pent) Ke Gad Bh 7) 


x G,(uz”, uy’, t”, t’) dt” du,” =8(u, —u,’)S(t—t’). (3.33) 


The function O(k, t) is a property of equilibrium, and can be calculated 
to the accuracy of the other approximations of this paper by various 
means, for example by the method of the preceding paper. If g(4,t) is 
defined by 

g(k, t) =exp (— xT /2m), Fed RO ETOL | 
then C(k, t) satisfies 


O(k, tt’) =4(g(k, tt’) + IC(k, 0) + w? | g(k,t—s)O(k,s)ds, (3.35) 


where w is the plasma frequency given by w%=47Ne?/V. This is solved 
by the Laplace transform giving 


O(k, 0) pt exp(pt) (Gl, P)+)) 9, (3.36) 
C(k,t—t')= ae p 2(1+0G,(p)/P) 
where = 
AUPE | exp(—pi)g(k,t)dt =. . . . (3:37) 
0 
and 


. a9 
an)=-v = (2). aaa al ganesh (3°38) 
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Since Q, is given by (3.9), the kernel of (3.33) contains 
Jexe (ike. ut)O(k, t) dh. 


Moreover it can be assumed that G, varies very slowly in time compared 
with a collision time, so that the time integral can be extended to infinity. 
This gives the kernel directly in terms of the Laplace integral 


P(k, u) = \3). exp (ik. vt — kx /2m) dt + 78(k . u)/2 + i77/2k. up 
0 


a -1 
nm {i+e? | texp (ik. ut — k?#?«T'/2m) a : See ofS 239) 
S 0 
We find 
6G, _ (d?k 0 Onc TAY ee yt 
aI + {Grek O)k. ie u)k. & + a G,(u, u’, t,t’) 


=§(u—u’)8(t—t’). (3.40) 


The plasma oscillations are represented by the denominator in (3.39) 
a fact noted by other authors, e.g. Balescu (1959), Hubbard and Thompson 
(1960). But consider for the moment only the term 6(k.u). 

This is the only surviving term for the case of a Lorentz gas, i.e. a gas 
of stationary scatterers which was used as the model at the start of this 
section. It gives 


0 d(k.u) Qo mu 
3h = 28a BOE Te 
5 [auc uel d le. du u (1+ w?m/ mT *- (Gat or) at 


for the integral in (3.40). Since 


C(k, 0)=NA++kK A+) and atm/cP =x 

this reduces to . 
ad ew eret: saa (3.42) 
(1+ k-A-*) kb u® aQ2 ‘ 


where 0°/0Q? represents the angular part of the Laplacian in velocity 
space. This integral converges at small ’, whereas in more elementary 
treatments, such as the Fokker—Planck equation, it diverges at small k 
and is customarily cut off at k*\?~1 which agrees with the convergent 
form above. It is interesting to notice that the factor (1+k-2\-?)4 
comes from two sources: the static sheilding effect contained in O(k, 0) 
giving one factor (1+~*A~*)-! and the dynamic effects of the motion 
of the charge through a collision represented by the time integration of 
O(k,t) giving another. If one uses the elementary theory the shielded 
potential appears twice in the cross section for scattering and over- 
estimates the effect of static shielding by a factor (1+k-A~)-. But 
the elementary theory ignores the dynamic effect, and in the end reaches, 
rather surprisingly, the same answer (3.42). This agreement however 


Transport Problems in Fully Ionized Plasmas 81 


is not present for any of the other terms of (3.39). The integral (3.42) 
still diverges. at large k, which comes from the inadequacy of using 
perturbation theory for large angle collisions, i.e. close encounters. I+ is 
customary to make a cut-off at k=2/ry where ry=e2(3«7')/2, so as to 
give the same behaviour as the Boltzmann equation which handles ‘hard’ 
collisions exactly. This cut-off could be removed in the present work 
using the exact two-body scattering, instead of the Born approximation, 
but this requires quite complicated analysis which is not matched by any 
physical interest and so will not be discussed further. 

The differential operator 0?/0Q? depends only upon angles because 
with all the ions at rest all collisions are elastic and only affect the direction 
of the velocity of the scattered particle, not its magnitude. The remaining 
terms in (3.39) due to electron—electron scattering are more complicated. 
For small v they contribute a term 


[ace 0)(1+h7A*)k. = exp [| —m(k. u)?/(2k?«7’) | 


a = mu\1 a(KT’) 
= + —]= 3.43 
xk (54+ o)5/ Ge ( ) 
to the integral in (3.40). 


For general u however the expression is complicated and shows resonance 
structure at w?~(k.u)?. 
For large w one can write the denominator 


1+ wf texp [ik. ut — k?#?(«Z'/m) | dt 
0 


1 w? iw?(k.u)m 


Te ee ee 
faye TPN mete amen) ate Ee 


which becomes small when w?~(k. u)?. To get some idea of the effect of 
the resonance one can consider its effect concentrated at the points 
w= +k.uand approximate that part of (3.39) coming from the resonance 


by 


du kT 


fa] m 
(n/16) | O(8 0) d3kk. es [5(w—k.u)+5(o+k.u)]k. ( aa ‘| 
(3.45) 


For large wu this again gives rise to a term like that of the Lorentz gas 


kdk 11 o nts 

oo) [pees A) a8 BOP" a) 
Crude as this approximation is it shows that nothing spectacular arises 
from the plasma resonance, and within the near-equilibrium analysis of 
this paper does not block the runaway nature of the conduction problem 


EF 
P.M. 
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which is implicit in the u? velocity dependence of (3.46) as will be discussed 
further below. The final equation for G, 


(5; + A(uy) Oa (uy tn’) = (uy — wy U= 8) eee yi 


where 
0 mu 


2G é pe 3.48 
Ses on | ET mead tscerh 


alk, u) =| (kw [1+ or aya (1 — Blk, 4) » | 
7m a 2mw? i | mk. >| 
ss (sm Tie) OF THe 


mie 2 
x [1+ as (1 Bik, uw) | 
ama4(k .u)? | mk. =|} (3.49) 
exp| — : 


PATS «Th? 


k k 2m (k. u)(m/2«Tk? 

B( eu) ( »), |e) | dx 
oats ae m(k .u)? 

exp | « aS : 


In the same way the equation for G,, becomes 


and 


(ar + A(u,) + A(u,) + B(uy, ug) + B(ug, u) Giy9(Uy, Uy’, Uy, Uy’, t,t’) 


= 6(u, — u,’)d(u, —u,’)d(t—-t’) 54 at tasoO} 
where 


ST in ek a od mu; 
BOW ed Eee k. Facto u—upk. (a + ae). (3.51) 


It is intended to give detailed calculations based upon these equations 
in a subsequent paper, but they will be illustrated in the next section by 
evaluating the case of the Lorentz gas. 


§4. Toe Lorentz Gas 


In the case of the Lorentz gas, the cross terms B in (3.50) are absent, 
and this equation has the solution 


Gis = GG. ° e ° e . . ° ° (4. 1) 


Hence the conductivity comes solely from the first term in (2.14) because 


[oem exp ( — H/T) du, dus Bu,’ du,’ = 0. . @ (4.2) 
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The equation for G, is 


ro] a 0 , , 1 , 
(5; ate = oh) Gi,(u, u ,t,t )=d(u—u )o(t—t’) 


(4.3) 
where 
y _ _ 2aNe* (70 kdk 
m* {. (kB 4-2)’ hed 
2nNet 
— log (2/r9+1), se le el WRT 
SP (2 e*/mn= log (2A) 7.) see eee (4.0) 


Since there are no dynamic effects in the strict Lorentz gas, only one 
factor (1+k-*A~)-! appears, but this suffices to make things finite. 
Since the differential operator in (4.3) does not involve |u| one can write 


G, = G,'d(|u] — |u’|) bed vo gt ema 4e7) 
and as 


d(u—u’) = ~3(c0s 6 —cos 6’)8(A—¢4')d([u|—|u’|]) . (4.8) 


the equation for G,’, becomes 


fa] Xo 0? ‘| 1 , , 
(5; + SSaqp) Ci’= 7290080 co8 884-49. auido) 


This is solved in terms of spherical harmonics Y,,,(9,¢) which are eigen- 
functions of the differential operator 07/0Q?: 


Gy! = SY exp [=ag(t HL + 1/8] Yin( #)Yin® OP )OUE—H). (4.10) 


Hence o,,, which is clearly isotropic 
On) = a6, 


is given, when the electric field is switched on at (=0, by 


Net / m \s2 fé 1 
sa adie 80088 ¥ — VimYimt w 
er (sen) | seo HOG arte 


x exp [—ap(t—t’ )U(L+ 1)/u3 —mu?/2xT].u? du. (4.11) 


The only values of /, m contributing are /=1, m=0, and as ¢ tends to 
infinity o tends to 


oN E (air), [, exp (—2avilu alder duexp (—mat/2<2 
0 0 


ae 3xT \2rnT 
2Ne*ar Vea Na IP 7 2 E 
—— utexp (—mut/2eT)du 6 ) aint eu (4kd) 
Bak L' (<7) 0 P( / 
=2(2«T |)3?(e2m"? log (22/79), we es ee fh 13) 


which is the customary solution. 
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In terms of the distribution function 


f=foth 
INGE tips Seat Per 
=f,+ — Yom» mm* &XP (— %l(E+ 1)(t—-#)/u?) 
ot Mesa) [  Yin¥act exp (— all +168) 


x u-*8([ul —|u’|) & . u’ exp (— mu?/2xT’) dv 
=f yt (Ne(& .u)u3/(2a9xT")) fo, ». » «6 gee eae 


when t—?t’ tends to infinity. 

It is interesting to note that the knowledge of G, allows one to see 
how f settles down to the value (4.14) when the electric field is switched on 
at the time t=0. This point of view can be used if one wishes to consider 
the runaway problem, for in this formulation a continually increasing 
current is represented by a conductivity which increases with time. The 
dependence of f, upon & was ignored from (2-5) upwards, if it is restored 
the equation for G, for the Lorentz gas becomes 

2 
(5+ Stam + ef 5a) G,=d(u—u')d(t—-t’). =. . (4.15) 

It is now clear that for large enough w the term in & dominates and such 
electrons are accelerated indefinitely giving a correspondingly increasing co. 
It is also the case that if G, were expanded in a series in 6 and the 
conductivity evaluated term by term, each term is finite, but the series 
diverges as ¢ tends to infinity. As has been noted above the plasma 
resonance in the realistic problem does not affect this conclusion within 
the accuracy of this paper. 

The Lorentz gas model also illustrates the modification due to a constant 
magnetic field, which gives rise to an external force e(u x B)/c (and to a 
modification of O(k, t) in the general case (3.40) which is negligible provided 
that the Debye length is less than the Larmor radius) so that 


0 é 0 eg (ch , ‘ 


Choosing B in the Z direction this becomes 


C0. OB ie iin OF , ' 


which has the solution 
Gy = (|u| —|u’|)u-2 2 Yim Y im*® EXP [(— ol (+ 1)/u3 + ienB/me)(t—t’) | 


x @O(t—t'). (4.18) 
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The calculation of o proceeds as before; o,, is unchanged, but o,, 
and o,,,, involving 1=1, n= +1, become 


Ong = Cyy = [ u® duu? exp (— 2apt’ /u®) cos (eBt’ /em) exp (—m U?/ QT) at’ 
J0 
a ie (=) uw exp (— mu?/2xT) as (4.19) 
SiGe) Leia (aine? 32) (tortor yee ee 


This method of handling the magnetic field is important when the 
Larmor radius and mean free path are comparable so that the transport 
coefficients cannot be expanded in terms of either ratio. For example, 
it is hoped to discuss in a future paper the off-diagonal elements of the 
pressure tensor which vary like B™ in strong fields and are difficult. to 
obtain by conventional methods. 


§ 5. CONCLUSION 


This paper has attempted to set up a formalism for transport processes 
in a plasma in which the equations used are direct representations of the 
physical processes involved. Thus eqn. (2.10) states that the conductivity 
can be related to the rate at which a single particle travelling through a 
plasma loses momentum and how the momentum lost by such a particle 
is in part transferred to other particles in such a way that it still contributes 
to the current flowing. The Green function G, represents the probability 
that a particle with velocity u’ at ¢’ will have u at ¢ and similarly for a pair 
of particles in G,,, and can be handled simply because of the immense 
simplification brought about by the predominance of weak collisions in 
a plasma. This is probably the only case in nature in which this condition 
really applies and here only within use of a cut-off at very large momenta. 
Perhaps the greatest simplification in matters of principle is however to 
consider only situations which are near equilibrium, because the use of 
the partition function as the first approximation transfers all questions 
of randomness to questions concerning the equilibrium state which are 
of course all answered by the use of the partition function. The 
assumption of molecular chaos or stosszahlansatz never has to be invoked, 
since such questions are referred back to the equilibrium state, whose 
status is well understood. 

The final eqns. (3.47) and (3.50) are still of course very complicated but 
it is hoped to give a quantitatively reliable evaluation of the transport 


coefficients from these in a subsequent paper. 


APPENDIX 


It is required to evaluate the static four-body correlation function of 
the charge density in thermal equilibrium and to show that the four-body 
correlation function is represented by the combination of products of 
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two-body correlation functions. If the Fourier components of the charge 
density are denoted by o;, 
o,= > exp(ik.x,)— > exp (ik.x,), en eet 
electrons protons 
then one wishes to evaluate 
<o Koy) and (04,0 j5/Fm > 
with the weight factor 
exp (— H/«T) =exp [— (47e?/v) > (0,0;,* — 2N)/k?xT’}. (A 2) 
The o,’s themselves are not good coordinates as they are an infinite set 
where there are N electrons and N protons, but it is possible to introduce 
a set of coordinates ¢; which are completely independent and allow the 
correlation functions to be written in closed form as integrals over the 
¢;'8. The details of this procedure are given in an early paper (Edwards 
1959) and here only the results of a straightforward extension will be 
given. If ¢(r)=> exp (tk.r)¢;, then using for convenience the physically 
equivalent grand canonical distribution rather than the canonical 
ensemble above, one gets 


Re | pe | E sper {2 [eos E ((r) — 4( (s)) |arrars! \ 
eee | (-1/8m2) | (va)? ar |T154. Re 
Wr | oh | exp | (— 1/82) i (va) | Tiss, (A4) 


(ote tee oe | nt ieee x a 
x exp {(2) | cos [e(4(r) — d(s))/KT ] @rars| | 


erp | (- 18m) | (v#)?@r |TT4, connie 


where Z = (277'/m)* exp (v/«T’), v being the chemical potential. 
In the first approximation of writing 


an 


2 | eos (e(r)/K7") dr =ZV(1—e?/V(KT)?) > bi.6_, 
=ZV—$yb_,|kTR  . . . . (AB) 


so that the exponent contains 
exp | — 3 (tla P x(t +k | 
k 
these formulae lead to 
{o;,0;) =N8(k +j)/(1+k-A-), hv ee SS 
To obtain higher approximations one has to consider the coulomb 


potential cut-off at some small radius to avoid an essential singularity 
inherent in the Gibbs distribution, but apart from such terms the 
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corrections to (A7) can be shown to be of order (/A)®. If one now 
proceeds to (c,0,0,0,,) there are terms coming from «=B=y=8 in 


exp [i(k. R,+j.R,+!.R,+m.R,)]) © . . (A8) 
which are static analogues of the failure of the Born approximation 
discussed after (3.42), but the remaining terms are again to order (/A): 


© (84,0 5%1%m) = (7 4%5 Fm) + (74%) (F7Fm) + KOKI m)<5;01)- (AY) 
For a detailed discussion of the handling of the short-range part of the 


potential and the evaluation of the ¢,, integration the reader is referred to 
the reference above. 
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ABSTRACT 


Two composite stacks of alternate sheets of nuclear emulsion have been 
used to investigate electromagnetic cascades in the cosmic radiation. The 
stacks were exposed in BOAC Comet aircraft for a total of 1400 hours at about 
38000 ft. Approximately 600 cascades were observed, and a sample of these 
was classified as arising either from y-rays incident from the atmosphere or as 
a result of nuclear interactions in the stack. The zenith-angle distribution 
of the cascades has been used to determine the attenuation length of the 
radiation producing the cascades. The number of electron tracks in the 
cascades has been determined by a photometric method, and the relation 
between the photometric density and the cascade energy has been established 
for the two types of cascade. 


§ 1. INTRODUCTION 


Dvrxinc the past few years, many papers§ have appeared describing the 
results of studies of energetic nuclear interactions or ‘jets’, produced in 
emulsion by particles in the cosmic radiation with energies in the range 
from 1000 to 100000 Bev. Most of the experiments were carried out with 
small stacks of nuclear emulsion, exposed for periods of only a few hours 
at balloon altitudes. The number of recorded events was therefore small, 
and it was difficult to establish or maintain any rigorous criteria for their 
detection ; the results were correspondingly uncertain. 

The great advantage of work with stacks of pure emulsion is that the 
tracks of individual particles can be followed for great distances, and their 
secondary effects observed in detail. Important advantages can be gained, 
however, if this feature is sacrificed by employing large stacks built up of 
emulsion and layers of heavy elements. Many of the nuclear interactions 
then escape observation, but it is possible to make a much more powerful 
attack on the detailed features of the electromagnetic cascades which result. 


eee ee a eee 
+ Communicated by the Authors. 
t On leave from University of Damascus, U.A.R. 


§ For review article see, for example, Perkins (1960). 
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from the nuclear interactions. Briefly, the advantages thus gained are 
threefold : 

First, in an appropriate assembly, the radiation length can be made so 
short that cascades of comparatively low-energy (a few hundred Bev) 
appear as concentrated bundies of tracks, visible to the unaided eye. 
Secondly, such electromagnetic cascades grow, reach their maximum 
intensity, and are attenuated, within the assembly; the precision with 
which their energy can be determined is thereby increased. Thirdly, 
the depth of the stack can readily be made large compared with the con- 
ventional value of the nuclear ‘interaction length ’—about 150 g/cem?, so 
that some features of development of the individual nuclear cascades can 
be investigated. 

To exploit the above advantages, two very large assemblies of lead and 
emulsion were exposed to the cosmic radiation in BOAC Comet 2E airliners 
for total period of about 1400 hours at an altitude above 30 000 feet during 
the period November 1957 to March 1958. Because of the length and 
intensity of these exposures, they allowed a statistically significant analysis 
to be made of some features of the interactions of particles with energies 
up to about 10’ ev. 

In the present series of papers, we present an account of the methods 
employed together with some of the results of an analysis of the events. 
This is followed by a discussion of the origin of the cascades and of their 
distribution in zenith angle, of the method of determining their energies, and 
of the energy spectrum of the cascades. The work relates both to cascades 
due to single gamma rays or electrons incident upon the stack, and to those 
associated with nuclear interactions occurring within the stack. 


§ 2. EXPERIMENTAL DISPOSITION 


Both assemblies, referred to hereafter as Comet 1 and Comet 2 stacks, 
were made up of vertical sheets of lead and emulsion. The proportions of 
lead and emulsion were varied, however, so as to give three horizontal 
layers of different mean density; see fig. 1. The top 7-5cm, designated 
section 1, consisted of pure emlusion made up of alternate sheets 400 
and 600 thick. The next 12-5cm, section 2, was composed of sheets of 
emulsion 600 1 thick alternating with 400 sheets of lead. The bottom 
layer of 10 cm, section 3, was made up of 600 u emulsions alternating with 
2400 y lead sheets. The sheets of lead were all wrapped in a thin layer of 
terylene, 6 thick, which prevented electro-chemical action between the 
metal and the emulsion. 

In the three sections, the effective values of the ‘radiation-lengths ’ 
were 29cm, 1-O0cm, and 0-6cm, respectively}. Using a value of 32cm, 
for the mean ‘nuclear interaction length’ in pure emulsion, the corres- 
ponding values for sections 2 and 3 are 23 cm, and 18cm, respectively. It 


a ee ee 
t In this paper, values of radiation length have been calculated from the 


formula given by Rossi (1952, p. 54) which assumes the validity of the Born 
approximation. 
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follows that for particles at vertical incidence, each stack presented a 
depth of 32 ‘radiation-lengths’ and 1-4 conventional ‘nuclear interaction 
lengths’. 

This particular-design of the stacks, incorporating several centimetres of 
pure emulsion at the top, was originally chosen to allow the observation 
of high-energy nuclear interactions in pure emulsion, and to ensure a rapid 
‘development, in the lead-emulsion layers, of any electromagnetic cascades 
resulting from them. It was also thought desirable to be able to observe 
the points of conversion into electron pairs of any high-energy y-rays 
incident from outside the stack, and to investigate the change in the cross 
section for bremsstrahlung predicted by Migdal (1956, 1957) and others 
(see Fowler e¢ al. 1959), 


Fig. 1 
ee) os ee 
Section 1 
1-4 nuclear interaction 
; lengths 32 radiation 
Section 2, lengths. 
7 1 | Section 3 
WAWHA |} 
<3.Omme Pb 
“unit cell repeated 85 times [J] Emulsion 


‘Cross section of Comet stack, showing disposition of lead and emulsion sheets. 
Taking the three sections individually, the values of radiation length and 
nuclear interaction length are respectively: section 1, 3-0 cm and 
32cm; section 2, 1cm and 23cm; section 3, 0-6cm and 18cm. 


The Comet 1 stack had a sensitive volume of 20 litres, with approximate 
dimensions 25cm wide, 26cm thick and 30cm deep. The dimensions of 
the Comet 2 stack differed only slightly from those of Comet 1; it was 
30cm wide, 24cm thick and 30cm deep. Both stacks were thermally 
insulated and packed in wooden crates. No cooling arrangement was 
installed for the first stack; the second was kept below 4°c with the aid of 
-solid carbon dioxide. 

In the Comet 1 stack, all the emulsion was of type Ilford G5 with normal 
gelatine content. In the Comet 2 stack, the 400. (even-numbered) 
emulsions of section 1 were x 2 gelatine, either G5 or K5. Section 2 of this 
stack was in two halves, the first incorporating G5 emulsion only, whilst 
in the second, every third emulsion was of type G5, the others being of K5. 
In section 3, all emulsions were of type G5. The reason for employing 
emulsions with different values of the gelatine content was to give the 
possibility of investigating nuclear interactions in light and heavy elements. 


92 J. G. Duthie et al. on 


Emulsions of type K5 were used because they were thought to provide 
better conditions for investigating ‘jets’, since their smaller grain-size 
gave improved optical resolution. Unlike the G5 emulsion, K5 was found 
to be particularly sensitive to surface marking from the terylene employed 
to separate the emulsion and lead sheets; as a result the K5 emulsions. 
could not be used for photometric measurements. 

The most interesting results obtained in the experiments were those 
concerning the energy spectrum of the cascades produced both by incident 
y-rays, and those arising from nuclear interactions occurring in the stacks. 
For such studies, a more uniform distribution of lead and emulsion, and the. 
use of ordinary G5 emulsion throughout, would have been preferable. 

The Comet 1 stack was exposed above 30000 feet for 840 hours, and the 
Comet 2 stack for 580hours. From the typical flight curves shown in 
fig. 2, it may be seen that the effective altitude of exposure was 11 500 metres. 
corresponding to a depth of 220 g/cm? within the atmosphere. Most of the 
total time of the exposures was obtained whilst the aircraft was flying 
between London and Beirut. 


40,000 


LONDON - BEIRUT 


30,000 x—x 1.1.1958 
o---o 13.12.1957 


20,000 


ALTITUDE IN FEET 


10,000 


gm /cm2 OVERLYING ATMOSPHERE 


ie) 1 2 3 4 5 
HOURS 


Height-time curves for two of the Comet flights from London to Beirut.. 
The return flights were typically half an hour longer. About 300 such 


flights were made for the total exposure of 14 i 
30.000 feet. xposure 0 00 hours above 


§ 3. SCANNING 


Cascades occurring in the stacks were detected by viewing the emulsion 
sheets by unaided eye against a suitable illuminated screen. The cascades. 
then appear as lines, which are, in most cases, readily distinguishable from 
scratches and other background marks. All lines thus detected were- 
examined under the microscope to confirm that they were in fact formed by 
cascades. In this way, about 200 events were found in the first, and 400 
in the second Comet stack. The first stack, although of slightly smaller- 
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volume, was given the longer time of exposure, and more cascades should 
have been recorded in it. The apparent discrepancy is due to the lower 
grain density in this stack, which was not cooled during exposure, and to a 
greater degree of surface marking. This increased the energy threshold 
for detecting a cascade by about 30%. The ‘plateau-value’ of the orain 
density for relativistic particles of charge || is 8-5 grains per 50 p in the 
first Comet stack, compared with 11 per 50, in the G5 emulsions of the 
second. The developed grains in the first stack, however, have an area 
about 10% greater than that in the second. 

The efficiency with which cascades are detected depends on their black- 
ness with respect to the background, and on their angles of dip with respect 
to the emulsion plane. Cascades of high energy may be detected with very 
high probability, even when they are very steep. Low-energy cascades are 
more readily visible when they are flat. The biases involved in the scanning 
will be discussed in more detail in a later paragraph (§1, Part Il). By way 
of illustration, fig. 3 Cf shows a contact print (natural size) of an emulsion 
containing seven cascades of various densities, and having various angles 
of dip with respect to the plane of the emulsion. 

After finding a cascade, it was traced back through a series of plates to- 
wards its origin, using the illuminated screen. When it became too faint 
to be distinguished by the unaided eye, it was followed under the microscope 
to its origin on to its point of entry into the stack. ‘Those cascades pro- 
duced in nuclear interactions in the stack, or in matter in its immediate 
vicinity, were classified as ‘nuclearevents’. Those produced by individual 
y-rays or electrons which entered the stack with no associated nuclear- 
active components were classified as ‘electromagnetic’ events. 


§ 4. ORIGINS OF THE HVENTS 
In a total of approximately 600 events detected in the two stacks by 
scanning, 300 events were followed back to their origins. This sample 
included all events with a cascade energy estimated to be greater than 
3000 Bev (see §§ 7 and 8), together with those of lower energy which occurred 
in conditions of favourable geometry and which were used in determining 
the energy spectrum of the cascades. 


4.1. Classification 
The origins of the cascades were classified as : 
(i) ‘ Electromagnetic’ events 
The cascade was followed back to either: 
(a) A single electron pair, originating in the emulsion, or occasionally in 
the lead of sections 2 or 3, with no neighbouring associated tracks. 
(b) The track of a single electron, entering the stack from the outside. 


The identification as an electron depended upon the observed rapid 
development of the cascade without evidence for nuclear interaction. 


+ Figures 3 and 9 are shown as plates. 
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(c) A number of photons and/or electrons, with no measurable angular 
divergences, which entered the stack within ~ 1004 of one another. 

In almost all cases, events of type (a) are due to single y-rays incident 
on the stack from the outside, originating directly from a nuclear interaction 
in the atmosphere; many of them will arise, for example, through the 
creation of 7°-mesons. In events of type (6) and (c) the y-rays will com- 
monly have materialized in the atmosphere before reaching the stack. The 
effects of multiple scattering and deflection due to the earth’s magnetic 
field, over distances of ~1km, which is equal to one radiation length 
(r.1.) at this altitude, can account for the lateral spread of the particles in 
events of type (c). 

In events of all three types, the maximum density of tracks in the cascade, 
as measured photometrically (§ 6), occurs at about 7r.1. from the point of 
entry into the stack. 


(ii) Nuclear events 
In these events, the cascade was followed back to either : 


(a) A nuclear interaction in one of the emulsion layers. 


(o) A nuclear interaction inone of the lead sheets interposed between 
emulsions. Such an event is characterized by a measurable convergence of 
the shower tracks and electron pairs in the first emulsion in which they are 
recorded. In all cases, it was required that the convergence measurements 
should show that the origin was inside the adjacent lead sheet. 

(c) A nuclear interaction in the ‘dural’ end-plates clamping the stack, 
in the wooden containing box, or in the aircraft fuselage. In such events, 
a close bundle of y-rays and shower particles entered the stack, and their 
measured convergences indicated the position of the origin. The nuclear 
origin was confirmed by secondary nuclear interactions in all except one 
of the events of this class, and the events were easily distinguished from 
those in (i) ¢, which arise from bremsstrahlung and pair-production in the 
atmosphere. 


(iii) “Families of cascades’ 

Occasionally, ‘families’ are observed, consisting of a large number of 
cascades, having no measurable divergence and laterally separated by 
distances of the order of a few centimetres. Most of the cascades in a 
family are ‘electromagnetic’ in origin, but a few are ‘nuclear’ events. 
Figure 3 (A and B) shows contact prints of some of the members of a large 
family (total 28) in the Comet 1 stack. Events of this type are attributed 
to very high-energy nuclear collisions in the overlying atmosphere, 
typically 1km above the stack, the individual y-rays having a divergence: 
of the order of 10->radians. They represent large air showers in the very 


early stages of development. The detailed features of such events will be 
discussed in a later paper. 
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The relative numbers of events of the different types are indicated in the: 
following table: 
Origins of all cascades of energy >3000 Bev in both stacks (total 75 


\ 
fo) 


Number observed 


(i) Electromagnetic events 
(a) Pure y’s 26 
(6) Close bundles of photons and/or electrons 11 
(ii) Nuclear events 
(a) In emulsion ll 
(5) In lead 22 
(c) In packing © 5 


4.2. Primaries of Nuclear Events 


In addition to tracing back all cascades of high energy, those of lower 
energy and of favourable geometry, were also followed, and gave an. 
additional 15 events of nuclear origin in the emulsion. Among all the 27 
nuclear events in emulsion, 18 were due to singly-charged primaries, 8 to 
neutral primaries, and one to an «-particle. This restricted sample of 
events does not provide any indication of a variation in the ratio of charged. 
to uncharged primaries with changes in the energy of the asscciated cascade. 
Thus, if the sample is divided into two groups, it is found that for a cascade 
energy 1500<#<3000 Bev, there are ten charged and five neutral 
primaries; for H>3000Bev, eight charged and three neutral primary 
particles. The large excess of charged over neutral particles cannot be 
attributed to a residua! excess of protons over neutrons in the nucleonic: 
component in the atmosphere. As will be shown below, the attenuation 
length of the nucleonic component is ~ 110 g/cm?, whereas the interaction 
length can be taken as ~70g/cm?. The calculated ratio of the numbers 
of true primary protons to the total number of nucleons at a depth of 
220 g/cm?, is then only about 18%. It follows that if the observations are 
interpreted in terms of pions and nucleons only, the observed excess of 
charged primaries is consistent with the view that 20% of all primaries 
are charged pions; the statistical uncertainty in the number of pions is. 
however of the order of 100%. 


§ 5. ZENITH-ANGLE DISTRIBUTION OF CASCADES 


Both stacks were flown with the emulsion sheets in a vertical plane, at an. 
effective atmospheric depth of about 220g/cm*. The zenith-angle distri- 
bution of cascades therefore gives a measure of the attenuation of the 
radiation producing them. Assume that the attenuation in the number of 
cascades of a given class varies exponentially with the mass of overlying 
atmosphere. Denoting the attenuation length by A, the vertical atmos - 
pheric depth by h, and the zenith angle by 6, it follows that the number of 
cascades per unit solid angle at 0 is given by 


dN dN(0) .exp [ —h(sec @— 1)/A]. 


dqaQ= dQ 
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The distribution with zenith angle of the number of the cascades has 
been calculated from the above relation for various values of h/A. From 
this formula, the distribution in the projected zenith angle, Op, in the 
emulsion plane, for different intervals in the angle of dip, 6, relative to the 
emulsion plane can be computed. This procedure is necessary as the 
probability of detection for the weaker cascades falls below unity for the 
larger angles of dip. Over a small range in the values of sec 5, the projected 
angular distribution will be independent of 5; variations in detection 
efficiency with 8 will then not affect the distribution in Op. Although the 
detection efficiency for the low-energy cascades is low, there is no reason to 
suppose that this efficiency will depend on the projected angle Op. It 
follows that for values of 5 < 30°, the projected angular distribution for the 
events can be closely approximated by the expression 


aN 
dbp 


where 1-05 is the effective value of sec 6 averaged over the events considered. 


ocexp [ — 1-05h (sec #,, — 1)/A] 


Fig. 4 
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Distribution of cascades as a function of zenith angle. 
(a) True zenith-angle distribution for energetic events where scanning losses 
were negligible. 


(b 


— 


Distribution of the zenith angle of the cascades, projected into the plane 

of the emulsions, for the more numerous cascades of lower energy 
(£'<3000 Bev). Only cascades with angles of dip below 30° are included. 
Although the scanning efficiency for the events varies appreciably with 
an euele, there is no reason to expect any dependence on projected 
angle. 
_ The two distributions are consistent with an exponential absorption 
in the atmosphere of the radiation producing the cascades, with an 
absorption length A=110g/em%. The figures for A given in the two 
graphs represent the best fits for each distribution. It is of course 
assumed that the attitude of the aircraft was horizontal and constant 
during the high-altitude part of the flight. 
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In fig. 4 (a), the true angular distribution of the 75 events of highest energy 
(£ > 3000 Bev) is plotted, as a function of cos 6. We are confident that in 
these events, the scanning efficiency was unity up to the highest angle of 
dip. A ‘maximum likelihood’ analysis for different assumed values of A 
gives a best fit for A= 105 + 16 g/em?; the computed curve for this value of A 
is indicated. 

In fig. 4 (6), the projected angular distribution is plotted for events of low 
cascade energy (1000 Bev <H<3000Bev), and of dip angle 6<30°. A 
‘maximum likelihood’ analysis gives A= 115 + 10g/cm?. 

It follows from the above results that the attenuation in the number of 
cascades, over the interval 220 g/cm? to about 400 g/cm? of atmosphere, is 
exponential, with an attenuation length of approximately 110 g/cm?. 
This value of A appears to apply over the whole range in the values of the 
cascade energy here considered, and to be valid for both ‘nuclear’ and 
‘electromagnetic’ events. The primary energies involved in the case of 
these cascades are in the range 104 to 10° Bev. 

The above value of is in good agreement; (i) with the measurements of 
Kaplon et al. (1953) on the ratio of the numbers of cascades produced by 
nucleons of primary energy > 1000 Bev at balloon and mountain altitudes ; 
they obtained A=130+15g/cm?; (ii) with the variation in the rate of 
extensive air showers between sea-level and mountain altitudes, and with 
barometric pressure or zenith angle at sea-level. The values thus obtained 
group around 110—-120g/cm? (for a review see Greisen 1956). There is 
indeed a remarkable consistency in the value of A obtained by different 
methods, over the very wide range in primary energy from 10°-10’ Bev. 


§ 6. THE PHOTOMETRIC MEASUREMENT OF CASCADE DENSITY 
6.1. Optical System 


In order to measure the energy of a cascade, the numbers of tracks must 
be determined at some given distance from the origin. When the density of 
tracks is very high, they cannot be individually distinguished and counted. 
Tn such circumstances, however, their number may be estimated by photo- 
metric measurements, which determine the opacity due to the assembly of 
tracks. The densitometer employed consisted of an (RCA 931) photo- 
multiplier tube attached to a modified Cooke M 4000 microscope. A CRO 
displays the signal from the photomultiplier. The optical system is shown 
in fig. 5. The axes of all slits in the system are aligned to the y-motion of 
the microscope stage, as also is the cascade under consideration. An 
image of the lower slit is focused into the plane of the emulsion by means ofa 

x 22 oil-immersion objective (N.A.0-65), used as a condenser. A x 22 
oil-immersion objective (of the same type as used for the condenser) is 
then adjusted to project the image of the lower slit onto the plane of the 
upper slit. This projected image is made to sweep across the upper slit 
by means of a rotating glass cuboid, driven by a small two-phase a.c. motor 
running at 1000r.p.m. ‘Two opposite faces of the cuboid were blackened so 
that at each revolution the image traversed the upper slit twice. 


G 
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The image of the lower slit in the plane of the emulsion had a width ee 
about 15, and a length of 200. The upper slit was adjusted to accep 
light coming from a region in the object plane, 10 pin width. : 

As can be seen from fig. 6, the grains contributing to the attenuation O 
light as registered by the photomultiplier are contained in the region 
ABCDEF. The figure here relates to the unprocessed emulsion, and the 
angle Omax is given by the formula 

N.A. 
ps ; 
where N.A. is the numerical aperture of both condenser and objective 
(0°65); py is the refractive index of the processed emulsion (1-5); and s is 


sim na60= 
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Diagrammatic representation of the optical system employed for photometric 
determination of cascade density. This system is an adaptation of 
that developed by Ney (1955) for measurements on the tracks of 
a-particles and heavier nuclei in the primary cosmic radiation. 


the shrinkage factor (the ratio of thickness of unprocessed to processed 
emulsion, approximately equal to 2). Thus Omax= 13°, represents the 
greatest angle of inclination of light rays to the normal to the emulsion 
plane, given by the condenser and accepted by the objective. 

Suppose now that the axis of a cascade—cylindrically symmetrical in the 
unprocessed emulsion—lies in the centre of the slit BE—the slit itself is 
normal to the plane of this diagram. A ray GH, normal to the emulsion 
plane, at distance « laterally from the centre of the slit, is attenuated 
according to the number of grains per square micron in the emulsion plane, 
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integrated over the emulsion thickness (600 mu). Aray JK, at an angle 0 
to the normal, will be attenuated according to the integrated number of 
grains per square micron at a distance x cos 6 from the axis, over a thickness 
equal to 600 sec @ microns. 

We have so far considered light rays lying in the plane of the diagram in 
fig. 6. Consider now the more general case of a ray making an angle @ to the 
normal to the emulsion, but with azimuthal angle 4 with respect to the 
plane of the diagram. Then one has to project the ray into this plane. 
The value of the angle %, which this projected ray makes to the vertical is 
always less than @. Since the values of cos @, or cos¢ lie between 1-0 and 
0-975 for all permitted light rays, it follows, for our purpose, that it is 

sufficient to consider all light rays as being normally incident on the 
emulsion. Thus, the arrangement employed measures the number of 
grains per square micron projected normal to the emulsion plane; in effect, 
we have to consider the attenuation in a parallel beam of light, 10 4 wide, 
produced by the tracks in the cascade. 


Fig. 6 
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An enlarged view of the distribution in light rays traversing the emulsion 
layer. The thickness has been taken as that of the unprocessed emulsion. 
A cascade with axis perpendicular to the plane of the diagram will have 
a track distribution possessing circular symmetry about its axis. The 
significance of individual rays is explained in the text. 


6.2. Relation between the Cascade Density and the Number of Tracks 


Let A be the cross-sectional area of a grain, and assume that a track 
consists of n grains per unit path length. Let f(x,y) represent the number 
of tracks per unit area in a plane normal to the cascade axis; we shall 
assume at present that this axis lies in the emulsion plane. The coordinates 
of a particular track relative to the cascade axis are designated w and y ; wis 
the coordinate in the plane of the emulsion. Then one has the following 


G.2 
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relation for the attenuation of the light intensity J, J, being the incident 
intensity 


(5) = -tAnsiow) SE Pe ee 
dy 
or 
+300u 
I(x) =I,exp(—An7(x)), where (x)= | f (x, y)dy. (2) 
—300yu 


Here, 7(x) is the integrated track density per micron, at a projected distance 
z. ‘To obtain the total intensity, J, incident on the photomultiplier tube, 
one finally has to integrate over the slit width, 1.e. 


Jady| "exp (—Anx(a))  (de/10), where Jy=J1,).10y. (3) 
—5u 


If one assumes, as a first approximation, that the projected track density 
z(x) does not vary over the slit width, then one obtains the important 
relation 

J=J,exp(—Anr), or In(JjfJ)J=Anz. . . «4 (4) 
Here, the values of A and are constants throughout the stack of emulsions, 
if the development is uniform, J, can be determined from a “ background ° 
reading, and 7 is the projected number of tracks per micron in the slit. 
The quantity In (J,//) will be called the cascade density, D. If variations 
of 7 over the slit are taken into account, it can be shown (see below) that, toa 
high degree of accuracy : 

D=I\n(J,/J)=Ant 


where 7 is the value of + averaged over the slit width. 


6.3. Effect of Angle of Dip of a Cascade 


There are three ways in which the angle of dip, 5, of a cascade, relative to 
the plane of the emulsion, affects the estimation of track density : 

(a) The projected grain-density is proportional to sec 6. 

(b) The integral in eqn. (2) is carried to the limits +300 cos 8 microns. 

(c) Since the slit image in the emulsion is 200 « in length, it is focused 

only on that part of the cascade in the centre of the field of view. 

The effect of (a) is to increase the projected grain density, and hence 
absorption, produced by a given number of tracks in the cascade, whereas 
(6) and(c) reduce it. From the form of the radial distribution—see Pinkau 
(1957)—it is found that the overall effect is to increase D by a factor of 


(0-1+0-9secd). This factor is therefore only slightly smaller than the 
effect of (a) alone. 


6.4, Effects of Multiple Scattering of Light and Saturation of Grain-size and 
Density 

In the derivation of the above formulae, two assumptions have been 

made: (i) That the multiple scattering of light back into the ‘allowed 

cone’ of the optical system may be neglected, and (ii) That the area of the 
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developed grains, A, and the grain density, n, of a track due to a particle 
of minimum ionisation are constants for all values of 7. 

(i) The magnitude of the effect of scattering has been investigated by 
increasing the numerical aperture of the objective, to 1-0, the aperture of the 
condenser remaining unchanged, so that Omax becomes ~ 20°; this increases 
the solid angle accepted by the objective by a factor of ~ 2:4. Even for 
very dense cascades (D~3 as measured with the original objective), the 
change in D amounted to less than 7%. Thus we are confident that, 
except possibly for the three most dense cascades considered in this 
experiment, multiple scattering of light is unimportant with the optical 
system employed. ; 

(ii) There are two possible reasons why the second assumption (A and n 
constant) might be invalidated. First, if the particle density is extremely 
high, the number of developed grains may no longer be proportional to the 
number of electrons, for there is a significant probability of a particular 
halide crystal being activated by more than one electron. Calculation 
shows that such saturation effects are completely negligible in G5 emulsion 
for the particle densities considered in this investigation. In emulsions 
with fewer and larger halide crystals (e.g. x-ray emulsions) such effects 
might become important. , 

A second source of error may arise from a failure to develop sufficiently 
all the activated grains in the core of dense cascades. Experiments carried 
out with intense beams of x-rays and f-rays, show that, with the normal 
development procedure we employ, the blackening produced in prolonged 
exposures is no longer proportional to exposure time. This is due to the 
fact that the amount of developer available in the emulsion is insufficient 
to develop fully the activated grains. Such effects will be important when 
the lateral spread of the beam of tracks extends over a distance of several 
hundred microns; if however the width of the pencil of tracks is small 
compared with the emulsion thickness, developer can diffuse inwards from 
all round the axis and replenish the shortage in the region of the pencil. 
Such replenishment, for example, takes place in the development of the 
tracks of heavy primary nuclei. For an extended beam of tracks, however, 
the effect of ‘developer-starvation’ results mainly in a reduction in the 
degree of physical development, and a corresponding decrease in the size 
of the developed grains. 

The results obtained with x-ray and B-ray exposures indicate however 
that the effect is unimportant for the cascades under discussion. As a 
check on the effect in the actual cascades, a visual count was made of the 
projected grain density (in grains per square micron) at various positions 
in several cascades, together with the results of photometric measurements 
in InJ. The regions of the cascades employed in this measurement were 
such that the projected grain density did not vary significantly over the 
slit width. Applying formula (4) for two readings, p, J, and ps, Je 


respectively : 
a z A(py— pi) =n J4/2 Spare Mote Ben | 
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The results of the measurements are displayed in fig. 7. It may be seen 
that the experimental relation between (p2—p1) and In J,/-/, is linear, the 
straight line passing through the origin; thus eqn. (4) is confirmed. The 
slope of the line yields a value for the effective cross-sectional area of a 
grain, A = 0-32? (for the Comet 2 stack). The contribution by grains on 
5-ray tracks must, of course, be included in the estimate of n, the number 
of grains/100 per cascade track. These extra grains amount to ~ 10% 
of the total. This estimate is obtained readily from measurements in 
other emulsions of the form of tracks of relativistic nuclei with Z ~ 26 
(Powell et al. 1959) and applies to contribution from 6-rays with ranges 
lu<R<300pu. Longer-range 6-rays contribute very little, as they 
leave the emulsion, and grains on shorter 5-rays will be counted as part of 
the grain density of an individual track. 


Fig. 7 


Po = Pp (number of grains per square micron ) 


Comparison between the observed number of developed grains due to a cascade 
and the photometric density in the same position in the Comet 2 stack 
p is defined as the projected number of developed grains per square 
micron through the emulsion thickness, J is the output signal from the 
photomultiplier tube. The projected area covered by developed grains 
throughout the emulsion thickness, per unit area, will be proportional 
to the logarithm of the transmission since double scattering of light 
can be shown to be unimportant. p,, J, and pg, J, refer to two positions 
relative to the cascade axis. The observed linear relation between 
(p2—p,) and In (J,/J,) shows that the mean area, 4, of a developed 
grain is independent of the grain density n. The slope of the line 
corresponds to a value of 4=0-32 py? (see eqn. (5)). In the grain 
counting, grains associated with low-energy electron tracks (mainl 
d-rays) account for 10% of the total number of grains, but grains Pa 
heavily ionizing tracks were not included. 
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The value of A above is in accord with direct measurements of the 
apparent grain size from the gap-length distribution in tracks produced 
by particles with ionization about four times the minimum value. 

A further check was made by determining the central density of a 
number of high-energy cascades (3000 to 6000 Bev), both at 7 r.l. distance 
from the origin, where the density is near its maximum, and at 10r.1., 
where the number of electrons per square micron in the core of the cascade 
is reduced by about a factor of 5. The average difference between the 
values of the cascade energies based on the two different sets of measure- 
ments was only 3%. 

Thus we conclude that systematic errors introduced, either by multiple 
scattering of light, or by variations in number and size of developed grains 
in the cascade tracks, as a function of cascade density, are unimportant, 
over the entire range of densities encountered in this experiment. The 
validity of the simple formula (4) is thus verified and it can be applied 
with confidence. 


6.5. Statistical Errors in Determination of Density D 

It is found in practice that an individual reading, J, on the oscilloscope, 
of the transmission of light at any point in the emulsion is subject to an 
r.m.s. deviation of the order of + 8%. These variations are due, firstly, to 
fluctuations in the number and grain density of background tracks, and 
their inclination to the slit; and secondly, to variations in the degree of 
surface marking. An individual estimate of D is obtained by taking one 
reading, J, on the cascade, and two background readings, J, one on each 
side of the cascade. Thus, the r.m.s. error in D=In J,/J is approximately 
0-10. This figure applies to the Comet 2 stack. In the Comet 1 stack, 
the surface markings on the emulsions were more troublesome, and the 
error in D was typically twice as great. As a result, only events with 
values of Dmax > 0-7 were accepted for the energy-spectrum in stack 1. 

For normal values of dip, ~ 20°, several independent readings of D 
can be taken in each emulsion; if the cascade is steep however—é = 45° 
for example—only one reading is possible, but the cascade will pass 
through more emulsions near the region of maximum density, and the 
statistical error in the determination of Dmax is almost unaltered. 


§ 7. RELATION BETWEEN DENSITY AND ENERGY FOR CASCADES 
INITIATED BY y-RAYS 
7.1. Relation between Maximum Density and Energy 
Having obtained a relationship between D and the projected track 
density, 7, the latter quantity must be related to the cascade energy. It 
is convenient first to consider cascades produced by single y-rays. From 
the curves of Pinkau (1957), which are in good agreement with more recent 
theoretical calculations by Nishimura and Kamata (1958), one can con- 
clude that the electron density f(x, ¥) in the ‘ target diagram’ of a cascade 
must be a function of the form 
f(a, 9) =£y'Gr x Eo) ey ch nw 8 (0), 


104. J. G. Duthie et al. on 


where g(r x Hy) is a function of the product Hy andronly. Here, Hy is the 
cascade energy, and r=(x*+y?)'*. Such a formula follows at once if 
most of the particles responsible for the optical density of the central 
region of the cascade are of high energy, so that ‘Approximation A’ is 
valid. 

From eqn. (2), the relation between 7(0)=7(x = 0) and Hy can now be 
obtained : ¥ 

+300 u PA 
tS { f(0, y) dy= 2 | E°9(r x Hy) .dr 
0 


—300u 


Zo 
= 2H | g(z) dz 
10) 


where z=r x H, and z)= Hy x 300 p. 

If the upper limit to the integral may, to a sufficient degree of approxi- 
mation, be taken as infinity (i.e. if the contribution for r> 300, can be 
neglected), then 7, will be strictly proportional to Hy (again only under 
Approximation A). In practice, only a small proportion of the tracks are 
excluded on account of the finite emulsion thickness, so that tT, is found to 
increase only slightly more rapidly than the primary energy £). 

A second effect is concerned with the finite slit width employed (10,). . 
As already explained, eqn. (4) assumes the track density 7 to be constant 
over the slit. If this is not so, then the light transmitted is a measure of 
the average value of exp(—Anr) taken over the slit—eqn. (3). This 
quantity is nearly equal to exp (— Anz) where 7 is the average value of 7 
taken over the slit. This approximation holds to a high degree of accuracy 
(better than 5%), provided the ratio of the intensities of light transmitted 
at the edges and at the centre of the slit is less than four. Such extreme 
variations of projected track density over short distances will only be 
found in cascades of extremely high energy and in their early stages of 
development. Where the effect begins to operate (H, > 4000 Bev att=5r.l.), 
the value of D as measured photometrically with a slit 10 « wide, results in 
a slight underestimation of 7, and this tends to offset the overestimation 
resulting from the finite emulsion thickness described above. Empirically 
indeed, it is found that the relation between the central density, D,, of the 
cascade at t= 7r.1., where it is close to its maximum, and 4, is practically 
linear. Briefly, H,=4000(D,)°*4Bev over the range 6000 Bev >E,> 
2000 Bev, for section 2 of the stack where the radiation length is 1 cm. 
If, instead of D,, one takes the maximum central density, Dmax, one finds 
E, = 4000D max Bev, no appreciable departure from linearity being expected. 

The general relation between cascade energy and Dyax is: 

31 000X 


Ee EES Dax | BRY, 


where X is the radiation length measured in centimetres, and An represents 
the total area of grains in square microns per cascade electron track per 
hundred microns, and includes contributions from d-rays. Dymax refers 
to the maximum density as determined with a slit of width equal to 


Ey= 
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10 microns. If we substitute X=1lem and An=7-7 square microns per 
1004 per minimum track, we obtain the previous relation applicable to 
section 2 of the Comet 2 stack. In order to determine accurately the 
energy of an individual cascade, we find that the most convenient 
method is in fact to determine the maximum central density, Dmax, since 
this value is the one associated with the smallest statistical errors of 
measurement. 


7.2. Energy Calibration from Examples of Neutral Pion Decay 

So far, a relation between #, and Dmax (the maximum central density) 
has been established on a theoretical basis. The method must now be 
calibrated. For low energies, Pinkau (1956) has determined values for 
HE, by performing scattering measurements on the tracks of the cascade 
electrons. Since his results showed agreement with the theory of 
Nishimura and Kamata (1958), it is clear that the curves in fig. 8 must be 
in agreement with the points obtained from his cascades ‘A’ and ‘D’, 
In fact, this result would be sufficient to calibrate the curves of fig. 7 over 
the entire range. This becomes clear if one remembers that essentially 
only two assumptions are involved in order to extend the distribution of 
cascade theory from low to high energies, provided Approximation A is 
valid in both regions. These are (i) that the cross sections for pair 
production and radiation remain homogeneous; and (ii) that multiple 
Coulomb scattering varies inversely as the particle energy. Effects such 
as the suppression of low-energy bremsstrahlung by high-energy electrons 
are completely negligible in these energy regions, as far as the development 
of the cascade at appreciable depths is concerned. 

A completely independent check on the determination of the cascade 
energies is provided by measuring the angular divergence of the two y-rays 
resulting from the decay of a single neutral pion, and determining the 
density of the subsequent cascades. The angle of divergence, y, of the 
two y-rays is related to the energy of the parent neutral pion by the 
equation 

135 a 
Ee er et) MOV io tree Gees (10) 
x 
where ¢ is the ratio of the energies of the two y-rays. In this formula, 
« need only be known approximately (for a value of ce 16, the value of 
y is only twice that for equipartition). In the material available to us 
there are three nuclear events, designated «8, P20 and 192a, in which 
only one or two neutral pions were emitted. It has been possible in these 
cases to give rather close limits to the total energy in y-rays, from the 
observed angular divergence of the photons. Consider, for example, the 
event 192 in the Comet 2 stack. In this case, a nuclear interaction in 
the ‘dural’ end-plate gave rise to a single energetic neutral pion, but no 
observable shower particles. The two y-rays by chance converted in pure 
emulsion (section 1) within | mm of each other, after passing through 7cem 
of emulsion. Their separation was then ~5y. The angular divergence 
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Fig. 8 
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Relation between the primary energy H, of a cascade initiated by a y-ray, 
and the central photometric density D. D is equal to In (/J9/J), where 
J is the amplitude of the photomultiplier signal with the cascade axis 
focused in the centre of the image of the defining slit (fig. 5), and Jp is 
the signal for a point in the emulsion remote from the cascade. The 
curves indicate the expected variation of D with Hy, calculated from the 
data of Pinkau (1957), at various depths from the point of conversion 
of the parent y-ray (4, 5-5, 7-0, 8-5 and 10-0r.1.). The lower envelope 
of the curves gives the relation between # and Dyax, and it can be seen 
that this envelope is very nearly a straight line ; H)=4000 Dax Bev. 
The relation between D and H, shown here refers to section 2 of the 
Comet 2 stack, where the radiation length was 1 cm, and the mean 
area of developed grains on a minimum track was An=7-7 square 
microns per hundred microns of track length. 

For differing conditions of development, the values of D, for a given 
energy, will be proportional to An. For given values of An and D, 
the true value of Hy is proportional to the radiation length. For example, 
in section 3 of the Comet 2 stack, where An is identical with that in 
section 2, but where the radiation length is 0-6 em, a cascade of maximum 
density 1:0 would be produced by a y-ray of energy 0-6 x 4000 Bev. 

The experimental points shown in the graph are as follows : 

(i) The three lowest points are taken from the cascades ‘A’ and 
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between the two axes could be determined and was found to be 
(741).10-rad: From the appearance of the subsequent cascade, it 
would appear that ¢ was close to unity; thus the energy estimated for the 
parent neutral pion (3800 + 600 Bev) will be independent of small errors 
in e. Photographs of this event are shown in fig. 9. At 7cm distance from 
the point of conversion of the y-rays (71.1. in section 2), the calculated 
separation of the cores is ~10yu. This separation is quite insignificant 
compared with the lateral spread of each cascade at this position ; hence 
the cascade density will be indistinguishable from that obtaining if the 
axes of the cascades had been exactly collinear. Since the density D is 
proportional to the energy, it follows that D will be the. same as that 
produced by one cascade due to a single y-ray of the same total energy 
(3800 Bev). This calibration point is plotted in fig. 7, together with 
events P20 and a8 which occurred in stacks of pure emulsion, and to which 
corrections were applied to allow for the different grain density, grain size 
and radiation length in the different medium. 

The four points at highest energy in fig. 8 refer to events where the 
radial track distribution was determined directly by visual counting on 
the cascades at depths of 3-4r.l. and the energy estimated from Pinkau’s 
curves. 


§ 8. ENerGy DETERMINATION FOR CASCADES 
ASSOCIATED WITH NUCLEAR INTERACTIONS 


8.1. Differences between Nuclear and Electromagnetic Cascades 


Let us now consider cascades produced as a result of nuclear inter- 
actions in the stack assembly. These cascades differ from those produced 
by single y-rays, in that: (i) several y-rays, distributed in energy, are 
involved; (ii) the axes of the individual cascades are laterally separated, 
according to the relative transverse momentum of the y-rays (typically 
200 Mev/c); (iii) secondary interactions may occur which reinforce the 
cascade from the original interaction. 

Both effects (i) and (ii) increase the lateral spread of a ‘nuclear’ cascade 
over that of an ‘electromagnetic’ cascade of the same total energy, if this 


‘D’ of Pinkau (1956) where the energy was estimated from scattering 
measurements on cascade electrons. The stack employed was of pure 
emulsion, and the energies have been reduced as indicated above. 
[] measured at 7r.1. gg measured at 5-511. 

(ii) The next three points refer to cases where the energy could be 
estimated from the divergence of the y-rays in 7°-decay. 

A ; Event 1920, Comet 2 stack (see fig. 9), measured at 7 ell. 

vy and ©; Events P20 (two 7's) and a8 (one 7°) respectively, in pure 
emulsion stacks, at 61.1. 

(iii) In the three highest sample points shown without errors, the 
energy was estimated by determining the total numbers of tracks 
inside 300 » and comparing this with the energy estimate based on the 


central density D. The measurements were performed at 7r.l. on 


three cascades in the Comet 2 stack (compare fig. 10). 
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spread is measured at the same depth ¢ (measured from the nuclear inter- 
action in the one case, and the stack edge in the other). To a first approxi- 
mation, the effect of a source spectrum of y-rays on the lateral spread will 
automatically be allowed for if the maximum density is determined. 
This maximum will occur about Ir.l. earlier however, than that of an 
electromagnetic cascade produced by a single y-ray of the same total 
energy. This result follows since the spectrum of secondary photons 
and electrons about 1r.1. from the point of conversion of a single primary 
y-ray, is similar in form to the spectrum of y-rays produced at the origin 
of a nuclear interaction. In both cases about 30% of the total energy is 
carried by the most energetic y-ray or electron, and the energy content of 
the low energy radiation is small. These results are obtained from the 
study of nuclear interactions occurring in the emulsion stack, and from the 
study of ‘families’ of parallel y-rays originating in nuclear interactions 
in the overlying atmosphere. The experimental evidence will be des- 
cribed in detail in a later paper. 

As a specific example, let us consider a cascade of total energy 4000 Bev. 
The most energetic y-ray will typically have energy 1000-1500 Bey. 
y-rays of energy 100 Bev or less will account for less than 20% of the total 
energy, and, at 6r.1., their axes will deviate laterally by about. 100 1 from 
that of the most energetic y-ray. Because of this source of lateral spread, 
the projected central density of the cascade will be reduced from the value 
obtaining when all the y-rays are collinear. (We at present neglect the 
effect of secondary nuclear collisions.) However, the total number of 
electron tracks within a fairly large radius of the main axis will be un- 
altered. Since the emulsions are of 600, thickness, it is convenient to 
measure the integrated track flux through a square of side 600 u, centred 
on the main cascade axis. (This flux will for all practical purposes be the 
same as that through a circle of radius 350 u centred on the axis.) Inside 
this radius, it can be assumed that the integral number of tracks will be 
equal to within a few per cent to that in a cascade produced by a single 
y-ray of the same total energy at the same effective depth. As explained 
above, this will correspond to a y-ray converting approximately Lr. 
before the star origin. 


8.2. Integrated Lateral Distribution of Nuclear and Electromagnetic 
Cascades 


In order to determine the energy in the soft component associated with 
nuclear events, we have measured the mean cascade density over a lateral 
distance of 300 on each side of the main axis. Using the values of 4 and 
n above (§6.4), this yields the total number N,, of electron tracks within 
350» radius of the axis. Ny was determined typically at a depth ¢=7r.l. 
from the star origin. From the curves of integral number of tracks as a 
function of (# x r) (Pinkau 1957), the value of the energy E of the equivalent 
single y-ray was then determined for r= 350 and t=8r.l. (measured from 
the point of conversion of this equivalent y-ray). We shall call this 
energy H/39), to denote the value obtained from the number of tracks 
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within a 300 box, at 7r.]. from the star origin. Hj) denotes the energy 
determined from the central density at 7r.l. from the star, which is the 
same as that for a y-ray of this energy at 8r.l. from the point of conversion. 
If secondary interactions could be neglected, the value of Ho, thus derived 
should be largely independent of the choice of t. 

The histograms of fig. 10 show the distribution in the ratio 


R= Lzis00 

7/0 
for cascades initiated by interactions in the stack and by single y-rays from 
outside. For cascades due to single y-rays, the mean value of R is con- 
sistent with unity, as it should be. The r.m.s. spread in individual values 
of & is c~0-3; this spread is to be accounted for in terms of fluctuations 
in cascade development, and errors of measurement. The majority of 
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R is the ratio of the cascade energy in individual events as estimated from the 
total number of tracks within 350, of the cascade axis, to that 
determined from the central density. The distribution of Rf is shown 
for a sample of 36 ‘nuclear events ’ and 28 ‘ electromagnetic events °. 


the nuclear events (~ 90%) have values of & similar to those for re 
y-rays; the distribution is however skewed to high R vaiues, and ~ 10 "S 
of the nuclear events give R values greater than unity by more than twice 
deviation, o, of the y-rays. 
oi ile values of # in the Hake events must be attributed to the 
effects of abnormal divergence of the individual y-rays at the point of 
production. Detailed inspection of these events clearly reveals the occur- 
rence of abnormally large values of transverse momenta. There appears 
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to be no significant variation of the mean value of & with energy, over the 
small energy range investigated, and its absolute value is only slightly 
greater than unity (R=1-11+ 0-05). This result means that the expected 
increase in R above unity, due to the effect of relative transverse momenta 
between individual y-rays, is approximately compensated by the increase 
of the central density resulting from the effects of secondary nuclear 
disintegrations, and that the effect on the energy estimate due to the source 
spectrum of the y-rays is correctly allowed for by considering the equi- 
valent y-ray to have converted Ir.l. prior to the nuclear disintegrationy. 

In order to determine the spectrum of the total energy in the soft 
component liberated in nuclear interactions, we have employed the 
maximum central density of the cascade, Dmax, this value having the 
greatest statistical weight. The energy derived from Dmax, treating the 
cascade as initiated by a y-ray—i.e. Emax = 4000 Dmax Bev—will in general 
be greater than either of the other estimates 7/39) or Hz) just considered, 
since the observed maximum of nuclear cascades occurs on average about 
. 9r.l. from the origin. This shift is of course due to the effects of secondary 
nuclear interactions. From the longitudinal development (see § 3, Part I1) 
of high-energy nuclear cascades it is found that the Emax is 110+ 0-05 
times the value H,,. This factor exactly cancels the factor F in the fore- 
going analysis. Thus one obtains the result that Hmax is equal to Lz/g9 
to within a few per cent, and is therefore, on average, a good estimate of 
the total energy liberated in the soft component in the original nuclear 
interaction. As emphasized above, this estimate will be incorrect in a 
few individual cases ; where the value of R is abnormally large, the energy 
Emax will be an underestimate; and where the effects of secondary inter- 
actions are severe, Hmax will clearly be an overestimate. 
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ABSTRACT 


The results of the measurement of the density of electromagnetic cascades 
described in Part I are discussed. For events initiated by single y-rays from 
the atmosphere it is found that the integral spectrum of cascade energy may 
be represented by a power law with exponent —3-0+ 0-20 for the range 
1200-8000 Bey. The corresponding exponent for cascades associated with 
local nuclear interactions is —2-9+ 0-2 in the same energy range. The y-ray 
spectrum may be used to calculate the pion spectrum in the atmosphere and 
this is compared with that deduced from measurements on p-mesons at 
sea-level or underground. Reasons for the resulting discrepancies are 
discussed. The flux of cascades associated with local nuclear events when 
compared with the nuclear flux deduced from measurements on extensive air 
showers yields values for the fractional energy, Ky, radiated as pions in 
collisions of nucleons in the energy range 10 000-100 000 Bev. It is 
found that K, decreases strongly with increasing primary energy. An 
upper limit to K, has also been deduced from observations on the longitudinal 
development of these cascades. 


§ 1. LyrropuctTIon 


TuHE study of high-energy nuclear interactions has often proceeded through 
the detailed investigation of individual events. A more indirect approach 
can, perhaps, yield some valuable complementary information. An 
example of this second approach is the study of the .-meson energy 
spectrum at sea-level. This energy spectrum can be interpreted in terms 
of that of the parent 7-mesons which have been created in the high 
atmosphere. Barrett et al. (1952), Pine et al. (1959), and Ashton e¢ al. 
(1960) have investigated this problem over an extensive energy range. 
The study of high-energy y-rays in the atmosphere is an alternative 
method of establishing the energy spectrum of 7-mesons. The rapid 
absorption of y-rays in the atmosphere limits the scope of experiments 


at sea-level and consequently long exposures at aircraft or higher altitudes 
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are desirable. The emulsion-lead sandwich assembly described in Part I 
is a very good detector of y-rays in the energy region above 1500 Bev. 
Furthermore the study of electromagnetic cascades associated with high- 
energy nuclear interactions, when combined with the primary nucleon 
spectrum, allows one to decide the fraction of the primary energy radiated 
in the form of 7-mesons. 


§ 2. EneRGy SPECTRUM OF CASCADES 


Applying the formula H = 4000Dmax Bev derived in Part I to both the 
electromagnetic and nuclear events, the integral energy spectrum of 
cascades has been determined. 
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Integral spectrum of maximum central density D of cascades produced by 
electromagnetic events. For explanation of the distributions A, B and 
C, see text. The cascade energy EH, is related to D by the relation 
Hy=4000D Bev. The absolute flux of cascades in number/cm?/ster/sec, 


can be obtained by multiplying the number of events shown in curve C, 
by 1-0 x 10-19, 


2.1. Electromagnetic Events 


In determining the densities, and thence the ‘ energy -spectrum z 
of the cascades, it is important to ensure that the effects of scanning loss 
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for the cascades of low energy are made as small as possible. For this 
purpose, the events have been divided into three categories : 


(i) As mentioned previously, we are confident that all events for which 
Dyax>90:7 have been detected in both stacks. This group of events 
will be referred to as selection C. The ‘integral’ distribution of the 


42 events of this class is shown by the open circles, and as steps in the 
distribution C, of fig. 1. 


(i) For somewhat smaller values of D, the events in the Comet 2 stack 
are subject to severe scanning losses, and the measurements of D, for 
reasons given in Part I, are not considered reliable. For the Comet 2 
stack, however, an analysis of the distribution of the measured angles of 
dip and of projected angle of the events found indicates that the scanning 
efficiency was high for angles of dip <30°, and for values of Dmax>0°5. 
It was therefore appropriate to make a second selection of events from the 
Comet 2 stack for which Dmax>0-5, 5<30°. 42 events fell in this. 
selection B, and their distribution is included in fig. 1B. 


(iii) To extend the spectrum to lower densities, only the events found 
in the most favourable conditions of observation in the Comet 2 stack 
(G5 section) are considered, i.e. those in layer 2, see fig. 2 (Part I) for which 
2°<6<11°; together with those in layer 3, where the radiation length is 
shorter, for which 10°<6<30°. Among about 100 events thus chosen, 
60 have measured values of D>0-23. These are plotted as distribution 
A of fig. 1. Conditions for detection are most favourable in layer 3. If 
events in this layer only are considered, the estimated number of events 
with D> 0-23 occurring in both layers is 75, a result which suggests that 
about 15 in layer 2 have been missed. We believe that the selection A is 
not seriously affected by scanning loss for values of D > 0-3. 


The above procedures result in the collecting rate (cm?. ster. sec) for the 
various selections being different. For selection C the whole collecting 
area of both stacks was used, and there was no restriction for the angle of 
dip. In selection B, only the second stack was used, and this stack 
represented 43% of the total collecting rate, when due allowance is made 
for the size of the stacks and the durations of their exposures. The 
further restriction §5<30°, that was imposed, results in an additional 
factor of 0-73 in the collecting rate, which is computed on the basis of 
the value of h/A determined experimentally from the observed distribution 
of projected angles. It follows that selection B corresponds to only 
31% of the total collecting rate of both stacks. In a similar manner, 
selection A is found to correspond to 10% of the whole. These factors 
for the rate of collection of events in the different classes are approximately 
in the ratios 10:3:1, in agreement with the numbers of events greater 
than a common energy that are observed. 

The best estimate of the fluxes over the whole range of energies can now 
be obtained by scaling selections A and B by the appropriate factors. In 
the range 0:50< D<0-72, selection B is used and added to the integral 


2, 
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number in selection C at D=0-72; this gives the results represented by 
the solid circles in distribution C. Similarly, for D< 0-5, selection A is 
added, and the results are represented by the crosses in distribution C. 

The resulting integral spectrum, whether expressed in terms of density, 
D, or energy, corresponds to a power spectrum, N(£)=KE”", with an 
index n= —3-1+0-25, for D> 0-3, corresponding to H>1200 Bev. This 
‘slope’ is very steep compared with the generally accepted index for the 
primary cosmic-ray spectrum (of the order of —1-8). It does not depend 
at all critically on the conversion factors employed above; for it may be 
seen that all three distributions in fig. 1 can be fitted by a line of this slope. 
It may be remarked that the two points of lowest energy fall below this 
line; at least a part of this fall must be ascribed to scanning losses,which 
become progressively more important in approaching lower energies. 
Selection C flattens off below D=0-7, and selection B below D=0-5; a 
comparison with selection A indicates that the fall-off in these two cases 
was due to scanning losses. 

For the spectrum C, fig. 1, the absolute vertical flux of y-rays/cem?/ster/ 
sec may be obtained by multiplying the number of events shown in the 
figure by a factor 1-0 x 10~-°;. further details are given in § 3. 

It should be emphasized that the large value of the index n for the 
integral-spectrum here reported refers only to the energy interval from 
about 1500 Bev. At some lower value of the energy, the spectrum must 
“flatten off’, since it is well known that for secondary particles in the 
energy region of a few Bev, the energy distribution of pions expressed as 
a power spectrum, has an index similar to that of the primary radiation 
producing them. 

Independent measurements of the y-ray spectrum, using emulsions, 
were also reported by the Japanese group at the 1959 Moscow Conference 
(Nishimura 1959). In experiments at mountain altitudes, they observed 
an integral spectrum of individual y-rays with an index — 2-8 in the energy 
range 500-1500Bev. On the other hand, at balloon altitudes, they 
observed an integral y-ray spectrum, for energies below 500 Bev, with an 
index of —1-8. Neither of these results can be directly compared with 


those presented in this paper, for the altitudes of observation were quite 
different. 


2.2. Nuclear Cascades 


An analysis similar to that described above has been applied to the 
cascades associated with nuclear events. The spectra are shown in fig. 2, 
the selection criteria being identical with those for the y-rays. Spectrum 
C is based on 39 events with D>0-72. Above this limit, the spectrum is 
markedly ‘flatter’ than that of the electro-magnetic events, with an index 
n=—1-:9+0-4. Spectrum B is based on 37 events in the Comet 2 stack 
only, and spectrum A on 97 events, of which 70 have measured values of 
D>0-23. For D> 0-3, spectrum A has an index of — 3-0 + 0-5, in contrast 
with the slope (spectrum C) at higher energy. Again, the lowest energy 
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point in spectrum A has two values, the upper one denoting the number 
calculated using only layer 3 of the stack (as in fig. 1). In the case of the 
nuclear events, the normalization factors for the three selections differ 
slightly from those for the electro-magnetic events. This is due to the 
fact that the nuclear events originate nearly uniformly throughout the 


volume of the stack, whereas the electro-magnetic cascades are initiated 
near its surfaces. 
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Integral spectrum of maximum central density D of cascades produced by 
nuclear events. For explanation of the distributions A, B and C, 
see text. The cascade energy EH, is related to D by the relation 
E,=4000D sey. The line D is the estimated flux of nucleons which 
would interact in the stack; the energy of these nucleons, as read off 
from the cascade energy scale, is to be multiplied by a factor 10. Thus, 
we expect 20 nucleons with energy exceeding 80 000 Bev to have inter- 
acted in the stacks. 

The absolute flux of cascades in number/cm?/ster/sec, can be 
obtained by multiplying the number of events shown in curve C, 
by 1:0 x 10-*. 


The ratio of the calculated factors for distributions C and B is independent 
of the assumed value of the interaction length, whereas the observed 
number of events in distribution A depends on the magnitude of the 
shielding provided by layers 1 and 2 of the stack assembly. Assuming 
the ‘ geometrical ’ value for the nuclear interaction length, the normalization 
factors for the distributions C, B and A are estimated to be in the ratio 
9:3: 1,imexcellent agreement with the numbers of events actually observed. 
It is therefore again reasonable to construct a composite spectrum shown 
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in spectrum C, as was done for the electro-magnetic events. This com- 
posite spectrum rather strongly suggests a ‘flattening’ in the slope of 
the spectrum with increasing energy. The best fit to the experimental 
results for the interval 0-3 < D < 2-0 corresponds to an index n= — 2-9 + 0-2. 

A summary of the slopes of the spectra for electro-magnetic and nuclear 
events in different energy intervals is given in table 1. 


Table 1 
Cascade energy E.M. events Nuclear events 
range (Bev) 
Selection C > 3000 —3:7+0°6 —1-9+0-4 
Selection A > 1200 —2:94+0°5 —3:0+0°5 
Composite spectrum 1200-8000 —3-0+0-2 —2-9+0-2 


§ 3. DIscUSSION OF THE ENERGY SPECTRA 
3.1. Introduction 


The principal results of the present measurements on the energy spectra 
of the cascades may be summarized as follows: (i) the integral energy 
spectrum of cascades initiated by individual y-rays incident from the 
atmosphere is of the form H-’, where B~3-0; (ii) the spectrum of the 
cascades associated with nuclear interactions is of the same form for 
cascade energies between 1200 and 8000 Bev, but at higher energies the 
distribution appears to be appreciably flatter; (iii) the proportion of all 
the cascades which are associated with nuclear interactions inside the stack 
is about 50%, except at the highest energies. where events with associated 
local nuclear interactions predominate. 

If we assume that the y-rays are produced in the decay of 7°-mesons 
with a differential energy spectrum of the form 

aN, =aH-4+) dH, 
it may readily be shown that the resultant y-rays have a spectrum 
(B+ 1) 
thus the spectrum of y-rays is similar to that of their parent 7°-mesons, 
but they are diminished in intensity by a factor of 2/(8+1). 

By similar arguments, the energy spectrum of individual electrons 
arising from conversion of these y-rays, and of y-rays arising from 
bremsstrahlung, may also be shown to follow a power law of the same 
‘slope’. The spectrum of ‘electro-magnetic’ events entering the stack 
from outside, will therefore reflect the energy spectrum of individual 
7°-mesons produced in nuclear interactions in the air above the stack. 

In contrast, the energy of the cascade resulting from a nuclear event, as 


determined photometrically, is a measure of the total energy released 
into the ‘soft component’ by the interaction. 


dN, = 


E-O) dE; 
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As the two spectra are very similar between 1200 and 8000 Bev, we 
conclude that there is little change in the ratio of the total energy 
appearing as electro-magnetic radiation as a result of a nuclear inter- 
action to the energy carried by the 7°—meson of highest energy. 


3.2. The Significance of the Spectrum of Individual y-rays 
3.2.1. Comparison of photon and muon fluxes 


From the observed y-ray spectrum, it is possible to calculate the pion 
spectrum in the atmosphere, under the assumption that all the y-rays 
originate as the decay products of neutral pions; if there are other sources 
of y-rays, we shall obtain an upper limit to the pion flux. This flux may 
than be compared with that calculated from observations of the muon 
flux at sea-level and underground. 

Suppose A exp(—2/Ap) represents the source distribution of charged 
pions above a certain energy, H, at vertical incidence, where x is the depth 
in the atmosphere (g/cm?), and Ap the attenuation length of the pion- 
producing component. The total number of pions throughout the 
atmosphere is then A[l—exp(—H/Ap)], where the total depth of the 
atmosphere corresponds to H=1000g/cm?. The source distribution of 
y-rays of energy exceeding H, at depth x, will then be 4/(8 + 1) exp (—2/Ap), 
where f is the index of the integral pion spectrum, under the assumption 
that the number of neutral pions is half the number of charged pions. 
The number of such y-rays (and/or electrons), vertically incident on a stack 
at depth h, is then 


Ai, 
7 exp jp) exp (=A) 
een eM a Aae ig 
where A, is.the attenuation length of the y-rays. It follows that the 
(vertical) pion flux can be obtained from the observed (vertical) y-ray 
flux by multiplying by a factor 
Ap—A 
q=o—) g+1) 


oY, 


(1—exp(—H/Ap) 
[exp (—h/Ap) — exp (—A/A,)] 
For Ap=110g/em?, A,=60g/em?, B=3, h=220g, one obtains Y=30. 
The corresponding pion spectrum is shown by the line B in fig. 3. 

This spectrum may be compared with the pion spectra deduced from 
observation on muons at sea-level or underground. In such a comparison 
it may be assumed tentatively, that pions are the principal sources of both 
y-rays and muons. Magnetic spectrometer measurements of the. muon 
energy spectrum at sea-level have been made by both the Cornell and 
Durham groups (Pine et al. 1959, Ashton et al. 1969). The results extend 
at present up to an energy of several hundred Bev. The pion spectrum 
can be derived from them by calculating the energy loss of the muons in 
the atmosphere, and the probability that the parent pions decay before 
interacting. Both groups of results yield a pion integral spectrum of the 
form: N(£)=9:5x 10-2 #-1%4/cm?/ster/sec, where H is in Bev, for 
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E<400xBev. This spectrum, extrapolated to the energies considered in 
the present experiment, is indicated by the line C1 in fig. 3. Within the 
energy range 1200-4500Bev, the two independent estimates of the 
absolute number of pions agree to within a factor of about 2. 

The Cornell underground measurements (Barrett et al. 1952) of the 
range spectrum of muons yield a pion spectrum, above 400 Bev, somewhat 
steeper (integral index — 2-0) than that deduced from the spectrometer 
results. Normalizing the two spectra at 400 Bev, the line C2 in fig. 3 
shows the integral pion spectrum deduced from these underground 
measurements. 


Fig. 3 
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The energy spectrum of all charged pions produced throughout the atmosphere 

as calculated by three different methods. 

Curve C1 is the extrapolation of the pion spectrum deduced from the 

muon energy spectrum at sea-level, as determined with magnetic 

spectrometers (Ashton et al. 1960, Pine et al. 1959). The muon 

spectrum has been measured reliably up to ~400 Bey, and up to this 
energy, the pion spectrum has an index of —1-64, 

Curve C2 shows the slope of the pion spectrum calculated from the 

Ae range spectrum underground (Barrett et al. 1952). The absolute 
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Although the pion fluxes deduced from the muon measurements on the 
one hand, and the y-ray cascades on the other, agree in order of magnitude, 
the discrepancy in absolute number, a factor 2 at 1500 Bev, cannot easily 
be accounted for in terms of errors of measurements. Much mote serious 
however, is the difference in the ‘slopes’ of the pion spectra. The present 
experiment suggests n= — 3-0 for H between 1200 and 8000 Bev, whereas 
the muon results suggest a value of n between — 1:65 and — 2-0. 

The pion spectrum deduced from the y-rays, B in fig. 3, has been 
calculated assuming the index 8 of the power law to be constant. If in 
fact there is a change in slope of the y-ray spectrum, below 1500 Bey, 
for example, then the conversion factor Q will be reduced over part of the 
spectrum. It should be noted that the attentuation length Neots the 
electron—photon component, depends slightly on the slope of the spectrum : 
thus as B>oo, A,>48 g/cm?, the conversion length (see, for example, 
Rossi 1952). The effect of a moderate reduction in the index of the photon 
spectrum near our lower limit (1500 Bev) will therefore result in an 
appreciable reduction in the estimate of the pion flux at this energy (see 
table 2). 


Table 2 
B A, (g/cm?) Q 
—3:-0 it, 32-6 
—2-5 62 25-4 
— 2-0 He, 18-0 


It does not seem possible, however, to account for the whole of the 
difference in the two estimates of the pion flux in terms of errors of measure- 
ments, or small variations in the index f of the spectrum. It should be 
pointed out that the estimate of the pion production throughout the atmos- 
phere as inferred from observations on the «-mesons is largely determined 
by the rate of pion production at atmospheric depths less than 40 ¢/cm?, 
and depends on the assumed value of the nuclear interaction length of 
the pions, A_. In the present experiment, however, the main contribution 
to the photon flux arises from nuclear interactions at a depth of about 
160g/cm?. There will, therefore, be serious discrepancies between the two 
estimates of pion flux if the source distribution of pions departs appreci- 
ably from the exponential one assumed (Ap =110g/cem?); or if either this 
source distribution, or A,, vary with energy. 


3.2.2. Effects of other sources of photons and muons 

Hitherto, we have assumed that the photons and muons arise only 
through the decay of pions. It is evident, however, that, in the high-energy 
region heavier particles, such as K*-mesons, on account of bnelt aie 
rest-mass and shorter lifetime, may constitute an important source of 
muons. In addition, most of any Ky! particles, because of their very short 
lifetime (10-!° sec), decay before interacting, and of these one-third provide 


y-rays through the decay mode K,!->27°. An approximate formula for 
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the probability that a particle will decay before interacting, assuming an 
interaction length of 70g/cm?, and an exponential source distribution 
corresponding to an attenuation length of 110g/cm? (a value which is 


not critical), is 
fx [1+8 x 10% pr/me)]™, 


where 7 is the mean lifetime in seconds, and p is the momentum. 

Tablo 3 shows the relevant decay processes of mesons and hyperons. 
In many cases the .-meson or y-ray is produced as the result of a succession 
of two, sometimes three, decay processes. The column headed ‘ Weighting 
factor’ gives the number of particles decaying in the appropriate mode ; 
it takes account of the branching ratio in cases where alternate decay 
modes occur; and it gives the relative production rates for the parent 
particles assuming charge-independence to hold at production. Column 
3 gives the probability that the decay chain will be completed in the 
atmosphere. 

Values of émax =H max/H and ¢min= 2 min/E, the maximum and minimum 
values for the ratio of the energy of a secondary particle to that ofits parent, 
are given in columns 4 and 5; they apply to the first decay in cases where 
a chain of decay is involved. 

For a particular decay mode, the ratio, R, of the flux of secondaries above 
a given energy to the flux of the parent particles above the same energy 
which decayed according to the particular mode, is given by the relation 


= n (€max)°t* 7 (€min)? 
B+ 1 €max — €min 


where n is the number of secondaries of the appropriate type produced 
in each decay. 

In cases where the ,.-meson or y-ray is produced as a result of a decay 
chain, the ratio of the flux of 4.-mesons or y-rays, to the flux of particles 
initiating the decay chain which decayed in the manner indicated, will 
be given by the product of appropriate values of R. This ratio is given 
in an approximate form in column 6. 

The final column shows the corresponding contributions made by the 
various decay chains to the flux of u-mesons and y-rays for the particular 
value, B=2. These results, taking account of the various weighting 
factors and decay probabilities, are the product of columns 2, 3 and 6, 
except in the case of the decay of A°-hyperons. There, production resulting 
ie the decay of &° and 3 hyperons is taken account of in the last column 
only. 

Other decay modes, such as K++7+ + 79°>y.++ 2y, contribute very little 
to the flux of »-rays or p-mesons. This follows because the decay 
probability f;, xf, for the production of ,.-mesons compares unfavourably 
with the other processes at high energy; and as a source of y-Tays a 


value of f,, for the decay probability has to be compared with 1 for the 
decay of K,! and 7°. 
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At high energies (> 1000 Bev), f,/f.~7. Thence, for B=2, using the 
values shown in table 3, column 7, it follows that the ratio of the number 
of muons above a given energy arising from K-decay to the number arising 
from -decay is 2-7 (N,+/N_+). This figure stands in sharp contrast 
with the figure for the ratio of the number of y-rays arising from the decay 
of K and z-mesons respectively, viz. 0:21 N pr /N,*; neither coefficient 
is very sensitive to the value of 8. It therefore follows that the contribution 
of K-particles to the muon flux is at least ten times more important than 
its contribution to the photon flux. 

It is also interesting to consider the possible effects of any hyperons. 
The relevant factors are included in table 3. It may be seen that, assuming 
charge-independence to be valid, the relative numbers of muons and 
y-rays above a given energy derived from this source is close to 2f,; almost 
the same ratio as if directly produced pions were the only source. Any 
difference in the spectra of muons and photons cannot therefore be ascribed 
to hyperon production. Hyperon production would be an important 
factor if, for example, in a high-energy interaction, the highest energy 
baryon were equally likely to be a hyperon or a nucleon, and if its energy 
were comparable with that of the primary. In these circumstances, the 
energy transformed into y-rays by decay in the atmosphere would be of 
the order of 5° of the primary energy, a value comparable with that of 
the most energetic photon resulting from direct pion production (see 

3.4). 

: It - interesting to note that the decay of hyperons results in a strong 
negative excess of the muons, the negative to positive ratio being 4:1. 
Accompanying these hyperons, however, there will, of course, be K® 
and K~ particles of strangeness +1; the net charge excess of the muons 
from both hyperon and K-particle decays will depend critically on the 
partition of energy between these particles at production. It may be 
remarked that the positive excess of muons (at all energies) is customarily 
explained as due to the preponderance of protons in the primary co smic-ray 
beam. In fact, if strange-particle production is an important process in 
high-energy collisions, the interpretation of the observed positive/negative 
ratio for the muons will be greatly complicated. 

Finally, direct nuclear production of y-rays may occur at high energies, 
through the process of ‘inner bremsstrahlung ’—this radiation being 
associated with acceleration and creation of the many charged secondary 
particles in the nuclear encounter. At present, it is not possible to assess 
experimentally the importance of this contribution ; general theoretical 
considerations suggest, however, that the contribution to the total y-ray 


flux would be very small. 


3.3. The Cascade Spectrum from Nuclear Interactions in the Stack 


The interpretation of the ‘nuclear cascades’ is complicated by the fact 
that cascades are apparently produced by both nucleons and pions, in 
comparable numbers, and possibly by other particles. This result, it 
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will be remembered, follows from the observed relative numbers of neutral 
and charged primaries of the ‘jets’ recorded in pure emulsion. 


3.3.1. Proportion of pion primaries 

It is of importance to try to estimate the proportion of pion primaries 
from the observed y-ray flux. Assuming all the y-rays to originate in 
7-decay, and N,:=2N,,, the ratio of numbers of charged pions and 
photons above a given energy, vertically incident on the stack from the 
atmosphere, is 

pa b+ —exp (= 2p —A,)[ApAn) - (Ap Ay) - An] 

1—exp (—2a(Ap—A,)/ApA,) - (Ap —Az) «A, 
where, as before, 8 is the exponent of the integral pion spectrum, and Ap 
is the attenuation-length of the pion source; A, is the attenuation length 
for y-rays above a certain energy, and is tabulated in table 3. A, is to 
be identified with the interaction-length of pions, since pion production 
by pions is accounted for in the source distribution. Thus, taking 
Ap = 110 g/cm?, A,=70g/cm?, and A,=57g/cem?, and an effective value of 
a =260g/cm?, one obtains R=5. This value is rather insensitive to the 
choice of parameters employed. 

Virtually all the energetic photons incident on the stack are detected, 
whereas, assuming the interaction mean free path to be geometrical, only 
~ 50% of the pions will interact at points such as to give a detectable 
cascade. Further, only a fraction, F, of the pion energy is transformed 
into y-rays. The relative numbers of cascades due to pions and y-rays 
respectively is then (R/2)(F)’. Upper and lower limits to F can be set 
as follows: The upper limit of (f')’ will be 4; that is, in 4 of the events, 
practically the entire energy will be transformed into y-rays if charge- 
exchange is an important process. A lower limit is F=4. However, 
even if one-third of the energy is radiated on the average, the quantity 
1/(F8) will certainly exceed 4, on account of fluctuations. Hence, the 
ratio of numbers of cascades produced by pions and y-rays should lie 
between 0-1 and 0-9. Thus, we cannot expect the number of cascades 
produced by interactions of pions releasing a given energy into the 
electromagnetic cascade, to be more numerous than the number of such 
cascades initiated by y-rays or electrons. In fact, the observed neutral 


to charged ratio of the primaries interacting in emulsion ($4, Part I) is 
consistent with an effective value of F of 0:5. 


3.2.2, The change in the index of the spectrum; inelasticity of the 
collisions 


One of the salient features of the energy-spectrum is the great pre- 
ponderance of local nuclear origins for the cascades of highest energy. 
Thus, all six cascades of energy > 8000 Bev were associated with nuclear 
events in the stack. This observation leads to the conclusion that the 
nuclear disintegrations responsible for these cascades of highest energy 
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were not produced by pions, but were, very probably, produced by nucleons. 
The integral energy-spectrum in this high-energy region has an index 
n= —1-9+0-4, a value consistent with the assumed index of the primary 
cosmic-ray spectrum. Using the flux values for primary protons given 
by the MIT group (Clark et al. 1957), and assuming that the nucleons are 
attenutated in the atmosphere with a coefficient of 110 g/cm?, the fraction 
of the energy transformed into the soft component in these nuclear inter- 
actions can be estimated. Considering the 21 nuclear events in which 
more than 4000Bev was transformed into y-rays, we find from fig. 2 
curve D, that there are 21 nucleons of H > 80000 Bev expected to interact 
and give visible cascades in the stacks. It follows therefore that typically 
it requires nucleons with energy ~80000Bev to produce cascades with 
energy ~ 4000 Bev, i.e. that the proportion of energy, K,, transformed 
into y-rays is ~5%. Further, the total energy transformed to all pions, 
neutral and charged, and taking fluctuations into account, is only ~ 10% 
of that of the primary particle. At lower energies, the slope of the energy 
spectrum changes appreciably, becoming steeper. This must in part at 
any rate be attributed to the effects of interactions due to pions. 
Empirically, the spectrum can be expressed as: 


N=8D-32414D-*6, 


and the two terms may, tentatively, be attributed to the contributions 
from pions and other particles respectively. Such a distribution for. 
interactions due to pions would follow if, as seems highly probable, the. 
fraction of the primary pion energy transformed into the soft component: 
is independent of primary energy. It then follows from the above formula 

that for D> 0-30 (H > 1200 Bev) ~ 30% of all nuclear interactions should 

be ascribed to pions. The fraction of neutral primaries predicted by the 

formula is not however in good agreement with that directly observed 

(§4.2, Part I). For 3000>#H>1200Bey, neutral primaries should be 

only 9% of the total; yet among 15 primaries observed, five were neutral}. 

For H> 3000 Bev, 25° of the primaries should be neutral; in fact three 

primaries out of a total of 11 were uncharged. 

Because of these apparent inconsistencies between the neutral/charged 
ratio actually observed at the lower energies, and that predicted by the 
foregoing equation, it seems improbable that the change in slope of the 
spectrum can be explained entirely in terms of the contribution of pion 
primaries. ; 

§ 4. Lonerrupinat DEVELOPMENT 


In any dense medium, where decay of mesons into muons is improbable, 
it is clear that practically the whole of the primary energy will eventually 
be transformed into the soft component via successive generations of the 
nuclear cascade. On the other hand discussion in the last section led to 


{ In a similar study in the same energy range, the presence of about 30% 


neutral particles among the primaries of nuclear events has again been 
confirmed. 
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the conclusion that, on the average, only a small fraction of the primary 
energy was transformed into y-radiation in each individual nuclear 
interaction. 

Similar low values of ‘inelasticity’ in nuclear collisions have also been 
suggested by the results of experiments in which the energy of a primary 
particle making a nuclear collision in emulsion was estimated from a 
kinematical analysis of individual interactions (see, for example, Perkins 
1960). Ifthe ‘inelasticity ’ in such collisions is indeed low, the longitudinal 
development of the ensuing nucleonic and electromagnetic cascade should 
show evidence of it. 


Fig. 4 


D 


DENSITY , 


—e Event 133 (nuclear orign) 7 


Event 136 (nuclear origin) 
—-o Event 173 ( ¥-ray) 


DISTANCE FROM THE ORIGIN (radiation lengths) 
Three examples of the variation of the central density D of cascades with 


depth from the origin. Event 136 shows a striking rejuvenation due to 
one or more secondary interactions. 


To investigate such an effect in our present ‘stacks’, we considered 
only the cascades of highest energy, for those of lower energy quickly 
pass below the threshold of reliable measurement. That indeed a high- 
energy secondary particle, capable of rejuvenating the cascade, is sometimes 
present is shown by the observation of ‘double-humped. : cascades. An 
example of such a development is shown in fig. 4. The limited depth of 
the stacks employed made it necessary to limit the study of the cascade 
development to 24r.1. from the origin, or through a distance corresponding 
to one nuclear interaction length in the medium employed (essentially, 
75% by weight Pb, 25% emulsion). In these conditions, any secondary 
particle of great energy emerging from an interaction should have interacted. 
within the emulsion in about 60% of the events studied. 

The longitudinal development of 21 nuclear cascades of high energy, 
and with suitable geometry, was investigated. Figure 5 shows the average 
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behaviour of the central projected density of tracks per micron and this 
may be compared with the longitudinal development of cascades nee 
by y-rays of a similar energy. On the average, as stressed in § 8, Part L, 
for nuclear cascades the maximum of the central projected density is 
delayed by several radiation lengths, as compared with the results for 
single y-rays, and its subsequent decline is, on the average, much less rapid. 
The ratio of the integral under the curve—see fig. 5—for nuclear events to 
the corresponding quantity for y-rays gives a measure of the ratio of energy 
transformed into y-rays in the nuclear cascade, to that liberated in the first 
interaction. More closely it in fact approximates to the ratio of the energy 
transformed into the soft component in about the first 3 r.1. to that liberated 


Fig. 5 
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The graphs show the average behaviour of cascades produced by 12 single 
y-rays, and 21 nuclear interactions respectively. All events above a 
certain threshold energy and of suitable stack geometry have been 
included. The much more rapid attenuation of the cascades initiated 
by individual y-rays is clearly demonstrated. The errors shown are 
deduced from the internal deviations from the average of the individual 
events. In summing the contributions from the different events, the 
mean of the observed densities at 4, 6 and 8r.l. has been normalized 
to unity, so that undue weight is not given to the events of greatest 
energy. Similar results have been obtained by Grigorov et al. (1959), 
who employed ionization chambers in an assembly mainly composed 
of iron, exposed at sea-level. 
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in between about 3 and 20r.1. (the energy liberated after 20r.]. will not 
have passed through enough matter for its cascade to have developed). 

The above approximation should be fairly accurate if all secondaries 
can be assumed to be collinear. Clearly, the effects of divergence will 
result in some underestimate of the total energy transformed into y-rays. 
The integral under curve | (fig. 5) is approximately double that under 
curve 2, so we conclude that the extra energy radiated in these secondary 
nuclear interactions is, on the average, approximately equal to or greater 
than that liberated in the primary nuclear interaction. If one takes a 
value of § for the fractional energy transformed into y-rays in the secondary 
events (probably correct for pion secondaries, but clearly an over-estimate 
for nucleon secondaries), the nucleonic cascade will be attenuated exponen - 
tially with a characteristic length equal to 3 interaction lengths. Thus, 
taking Aint= 160 g/cm? in the composite medium of the stack, one finds 
that, in order to liberate an additional energy in the soft component between 
3 and 20r.1., equal to #,,, the amount liberated in the first interaction, the 
total energy carried by the secondary particles must be at least 6£#,,. 
The energy carried by the charged pions would be expected to be on average 
approximately equal to that in neutral pions from the first interaction 
(#,), and not twice as great, on account of fluctuations, so on its own 
could not provide the degree of rejuvenation that is observed. This 
behaviour is consistent with an inelasticity, in the first interaction, 
K,~20%. This value for K, is certain to be an overestimate, however, 
for if the primary interaction by a nucleon is inefficient at producing y-rays, 
it would seem very likely that the high energy secondary nucleons or 
hyperons were also inefficient. In order to get an estimate of this effect 
a much deeper stack is required. Furthermore, some of the primary 
interactions will themselves have been produced by pions, and in these 
no large secondary contributions are to be expected. Thus, the value 
of K, appropriate to nucleon primaries is certainly less than 20%. 

In the simplified analysis considered so far, no account has been taken 
of the angle of divergence of the individual particles. This clearly reduces 
the central density of the cascade—and is most important at large distances. 
Any formal analysis awaits more detailed knowledge of the transverse 
momenta possessed by the secondaries carrying the highest energy. 
Experiments are now being made with a stack exposed more recently, 
of similar size to the Comet 1 and Comet 2 stack, but of much higher mean 
density (12g/cm2), and this should provide better conditions for investi- 
gating the longitudinal development of nuclear cascades. 


§ 5. CONCLUSIONS 


It is not possible from the present data to give an unambiguous inter- 
pretation of the results relating to the flux and energy spectrum of the 
cascades. However, it seems that the following main conclusions can be 


drawn. 


P.M. 
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(i) The integral pion spectrum at 220g/cm* depth in the atmosphere, 
between the energy limits 1500-4000 Bev, follows a power law. The 
index is —3-1+0-3. The discrepancy between this value and that 
obtained from the muon measurements (—2-0)—appropriate however to 
pions or other particles produced near the top of the atmosphere—remains 
to be explained. 

(ii) The integral energy spectrum of cascades produced by nuclear 
interactions (mainly in lead) has an index — 1-9 + 0-4 for cascade energies 
from 4000-30000Bev. The primaries of these events are very probably 
nucleons. The fractional energy radiated into the soft component in 
these interactions appears to be K,~5%. At lower cascade energies 
(1500-4000 Bev), the spectrum is much steeper, with an index equal to 
that for the y-rays. It seems clear that this change in slope cannot be 
ascribed to pion primaries. 

(iii) The observed longitudinal development of nuclear cascades is 
consistent with the conclusion that in the interactions of nucleons producing 
cascades of energy ~ 4000 Bev, only a very small fraction of the primary 
energy is transformed into y-rays in the first interaction. 


Clearly, much more experimental data will be required before a more 
detailed analysis of such results is possible. In particular, it will be most 
important to extend the y-ray spectrum to both lower and higher energies, 
and to investigate the longitudinal development of nuclear cascades— 
especially those initiated by nucleons—at greater depths from the origin. 
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ABSTRACT 


Measurements of the capacity of the electrical double layer at the interface 
between a polarizable mercury electrode and an aqueous solution of a strong 
electrolyte are interpreted theoretically. Attention is concentrated on the 
layer of charges at the actual phase boundary, whose behaviour has not been 
satisfactorily explained before. 

It is suggested that for anodic polarization an increasing number of mercury 
atoms is situated at adsorbed sites on an otherwise flat surface. These carry 
a charge, and lead to a rise in the measured capacity as the polarization 
becomes anodic. It is thus shown why the capacity as a function of polari- 
zation has a form for fluoride electrolytes, where there is no specific adsorption 
of the anion, so similar to that for other electrolytes, where there is. 

The behaviour of the water molecules in contact with the mercury is con- 
sidered in some detail, and compared with the behaviour of bulk water and of 
a monolayer of water adsorbed to a mercury surface. It is shown that the 
hump, which appears in most capacity curves slightly to the anodic side of the 
electrocapillary maximum, may be due to the high polarizability of the water 
layer between the metal and the outer Helmholtz plane, when the field there 
is small. It is not necessary to assume as Grahame does that there is any 
appreciable chemical or image force polarizing the surface molecules, or that 
the water takes up a rigid ice-like structure for low fields. 

The suggestion that the rise in the capacity for cathodic polarization is due 
to compression of the double layer by electrostatic forces is also considered, 
with the implication of this idea on the mechanism of overvoltage in the dis- 
charge of cations onto mercury. 


§ 1. INTRODUCTION 


MEASUREMENTS of the capacity of the interface between mercury and 
aqueous solutions of strong electrolytes, using a dropping mercury 
electrode, have recently been made with considerable precision: those of 
the late Professor D. C. Grahame are of special importance. They have 
demonstrated that the electrical double layer at the interface falls into 
two parts: a diffuse layer of ions analogous to the long-range screening 
charge of an ion in the solution, and an inner layer, between the metal 
surface and the plane of closest approach of hydrated ions to it, called the 
outer Helmholtz plane (Grahame 1947). This is sketched in fig. 1. The 
capacity of the inner layer is roughly in series with that of the diffuse layer, 
and it is much smaller except in very dilute solutions so that it dominates 
the measured capacity. 

The forms of experimental capacity curves for strong electrolytes, for 
example due to Grahame (1947, 1951, 1954, 1957a), are rather similar 
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and after the capacity of the diffuse layer has been removed by the method 
used by Grahame (1947), they have the form shown in fig. 2. Here o is 
the surface charge on the mercury and C is the differential capacity da/d¢, 
where ¢ is the potential of the mercury relative to a standard electrode in 
the solution. 


Fig. 1 
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Sa atts values of C®, the differential capacity of the inner part of the electrical 

ouble layer, between mercury and aqueous solutions of sodium fluoride 
at several temperatures. Calculations made from data obtained with 
0-8 N solutions. 


Most anions can be ‘specifically ’ adsorbed to the mercury and are then 
partially dehydrated. They have their centres on the inner Helmholtz 
plane in fig. 1. Their concentration in the plane is proportional to their 
concentration in the solution, and since they are charged, o and C depend 
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on th rati i 
a e concenti ation. However, C is found to be practically independent 
a oo alee if the anion is fluoride (Grahame 1954), so fluoride ions 
n S ifi aAleon Odi i i | 
not be specifically adsorbed. Sodium fluoride is therefore a simple 
electrolyte to study. 


re 
§ 2. THe Evecrric Frevp In tHe Meran 


The electric field at the interface penetrates into the metal. Rice (1928) 
applied the Thomas—Fermi approximation to this problem, and showed 
that for mercury, with two free electrons per atom moving in a medium of 
mercuric ions, the depth of penetration is about 14. This medium has a 
low dielectric constant (Pauling 1927) as is pointed out by Rice, so that 
the layer in the metal would be expected to have a capacity of at most 
16 F/cm?, even if the solution behaved like a perfect conductor. How- 
ever, the measured capacity in moderately concentrated aqueous solutions 
of strong electrolytes is never less than this, and it depends strongly on 
the temperature. Temperature-dependence cannot be due to the 
electrons, so the layer contributing to the capacity must extend into the 
solution. It is shown in the appendix that if the positions of all the 
charges in the layer are known, the potential drop across it is the same as 
if they were smeared out in planes parallel to the layer, which can therefore 
be regarded as a number of condensers connected in series. However the 
field on an individual charge is by no means the same as if the other 
charges were smeared out, so the equilibrium distribution of charges 
depends on the potential across the whole layer, not just the drop in the 
metal. It is for this reason that the penetration of the field into the 
metal does not imply that the capacity of the whole layer is always less 
than 16yF/cm?. This is considered further by Mott and Watts-Tobin 
(1961). 


§ 3. ADSORPTION OF IONS 


The capacity rises sharply when the polarization is anodic and it depends 
strongly on temperature there. Measurements with 0-916 molar sodium 
fluoride (Grahame 1954) show that the sharp rise occurs as the potential 
approaches the exchange potential for a normal calomel electrode. This 
suggests that it is associated with the electrochemical reaction which 
would occur on the mercury at even more anodic potentials. The capacity 
curves for sodium fluoride are remarkably similar to those for electrolytes 
whose anions are specifically adsorbed, for example that for sodium 
nitrate (Grahame 1957a). The adsorption to the interface of positively 
charged ions from the mercury, or of negatively charged ones from the 
solution, would explain this, and we shall consider these possibilities. 

Gerischer (1958) has observed the current during sudden changes of 
polarization in electrode reactions. It provides evidence that on other- 
wise flat surfaces of solid metals in contact with electrolytic solutions con- 
taining their own ions, some of the metal atoms are adsorbed as ‘ad-atoms’ 
in a partially charged state. We therefore suggest that mercury atoms 
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are pulled out of its surface by the electric field to form ad-atoms. 
Gerischer and Krause (1958) have shown that the exchange current of the 
charge transfer between a mercury electrode and a solution containing 
mercurous ions is almost unaffected by freezing the mercury. From the 
fact that this current is almost as great on solid mercury as on liquid 
mercury they deduce the existence of ad-atoms on the solid, for if they did 
not exist charge transfer could presumably occur only at kink-sites and the 
rate would be much less. But the fact that the currents are almost the 
same suggests that ad-atoms may be present on the liquid as well. 
Otherwise there would be a considerable change in the number of atoms 
which can participate in the charge transfer when the surface is melted, 
unless the solid mercury surface is practically covered with ad-atoms at the 
melting point. This is unlikely, because the exchange current on mercury 
(Gerischer and Krause 1958) is less than that on silver (Gerischer and 
Tischer 1957), and there is only about one ad-atom per hundred atomic 
sites on silver. 

The idea of a liquid metal surface with ad-atoms of the same metal on it 
seems to be new, but it is possible thermodynamically: if there are few 
ad-atoms their entropy will be that of a two-dimensional gas, and it 
follows from the Sackur—Tetrode formula that the entropy per ad-atom 
increases indefinitely as the degree of coverage tends to zero. This is also 
true if the ad-atoms are immobile on the surface. If this entropy is S, 
and that per atom in the liquid is S, which is a constant, then at constant 
temperature 7’ and pressure we may expect ad-atoms to be adsorbed 
until 7'(S — Sy) is reduced to the value of the heat of adsorption from the 
liquid. If there is no actual energy barrier to cross in pulling an ad-atom 
out of a liquid surface, there will be a rapid fluctuation of atoms in the 
surface in various stages of this process. However, their average contri- 
bution to the potential can be described by means of definite ad-atoms. 

The distance from the metal of the inner Helmholtz plane in which the 
ad-atoms is not known, nor is their charge. Only the product of these, 
the dipole moment of the ad-atom and the charge induced on the metal by 
it, affects the potential drop from the metal to the solution. This is also 
true for anions adsorbed from the solution, although it is more convenient 
to consider the dipole moment of the ions and the charge induced by it on 
the outer Helmholtz plane. There is no reason why the distance of 
such ions from the metal should be the same as the distance of ad-atoms 
from the metal. If both ad-atoms and adsorbed anions are formed, they 
will cluster round each other to form.complexest. 

Suppose the dipole moment of a complex is 4, positive in the direction 
from the metal to the solution, the potential drop from the metal to the 
solution is 4, and the thickness of the double layer, namely the distance 
from the conducting surface of the metal to the outer Helmholtz plane, is 
d, ‘Then if the distance between the complexes is large compared with d, 


} This was pointed out to me by Dr. Gerischer. 
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the electric field in the neighbourhood of one of them. excluding its own 
field, will be approximately constant across the layer, so it will be ¢/d. 
This was first pointed out by Ershler (1946). The number of complexes 
per unit area of the surface will therefore be 


LG d—W 
n= N (2) exp ae ay 


In this equation V is the number of adsorption sites per unit area of the 
metal surface, and W, the heat of adsorption of a complex in the absence 
of an electric field, but including the image-force attraction. Wis 
supposed to be measured when the concentration of anions in the solution is 
cg, their concentration in this experiment being c,, and m is the number of 
anions in a complex. The complexes give a contribution to 4 of magni- 
tude 


= ATTN, fy [Ky Co elt, As wee roe nS (1 ) 
where K, is the dielectric constant of the layer. 


$4. THe Water IN THE INTERFACE 
4.1. Dielectric Theory 

Slightly to the anodic side of the electrocapillary maximum, where the 
mercury is uncharged, there is a hump in most graphs of C against potential 
or charge. We propose to show that this may be due to an increase in the 
effective dielectric constant of the layer, marking the point where the 
average dipole moment of the water in the layer changes sign. Of course, 
for this to be noticeable as a hump it must be possible to polarize the layer 
so strongly that an appreciable proportion of the water dipoles is orien- 
tated by the field. The average field H necessary to orientate a water 
molecule with dipole moment p is H~kT'/u. At 300°K, # is 6-6 x 10° 
volts/em. In capacity measurements the potential can be changed by 
more than a volt and it acts over a distance of only a few angstroms, so 
this field is reached. 

Such fields are much stronger than any which can be applied to bulk 
specimens of water. The water molecules in the interface will also 
behave differently from those in bulk owing to their proximity to the 
mercury, and to the fact that the electric field is applied not through the 
whole solution but by the double layer of ions near the mercury and the 
balancing charge on it. All the same, it is worth considering the theory 
of the dielectric constant of water in bulk to see how it must be modified in 
the interface. 

This theory was given by Kirkwood (1939) and applies generally to 
associated polar liquids. For a substance whose dielectric constant K is 
as high as that of water his formula takes the form 
N3yu-p* 


V3 347° e ° e e e¢ (2) 


K=K,+47 
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In this formula, there are V molecules in a volume V, each having a dipole 
moment p in the absence of an applied field. p. is not quite the same as 
the moment of a molecule in vacuo, owing to its polarization by the ‘reaction 
field’ when it is in the liquid; however, the difference is small. w* is the 
dipole moment of a large sphere of the liquid surrounded by the same 
liquid, when one of its molecules is fixed with a dipole moment wp. The 
interaction energy of the liquid in the sphere with the field C in the cavity 
is —p*-C. If the sphere is large enough to treat its surroundings as a 
continuum of dielectric constant K, and the cavity is filled with a con- 
tinuum of dielectric constant K, (representing the polarizability of the 
liquid when the dipoles in it are fixed in definite directions), then 


C=3KE/(2K + K,)3E/2, 


where E is the field in the surrounding liquid. The factor 3 in the denomi- 
nator of (2) arises from averaging the directions of p. 

The high dielectric constant of water is due to the fact that |z*| is con- 
siderably larger than |u|. The dipole moments of neighbouring water 
molecules are correlated by their hydrogen bonds, and if this correlation 
extends only to nearest neighbours, .*=7u/3. The theory only applies 
when the electric field is too weak to produce dielectric saturation. To 
calculate the saturation of the polarization in stronger fields is much 
harder, as the field in the cavity is no longer 3KE/(2K +K,), and the 
averaging of the magnitude and direction of u* are no longer independent. 
Only the first term or so in an expansion of the dielectric constant in 
even powers of E can reasonably be calculated (Buckingham 1956), but 
they contain all the information which can be measured with fields that 
can in practice be applied to bulk specimens of liquid. 


4.2. Structure 


Pauling (1939) gives the energy of the hydrogen bond in water as 
4:5keal/mole. It is pointed out by Lennard-Jones and Pople (1951) that 
this is so much greater than 7’, that at any instant of time the hydrogen- 
bonding in it is practically complete. The strength of the electrostatic 
field H necessary to break an appreciable number of hydrogen bonds is 
given by NEu~4-5keal where N is Avogadro’s number, so #~5 x 107 
volts/em. The maximum fields applied across a mercury-solution inter- 
face are almost certainly less than this, though perhaps not much. It is 
therefore a reasonable first approximation to suppose that the hydrogen 
bonds between water molecules there are all intact, except for the water 
molecules actually in contact with the mercury, and those in the primary 
hydration shells of ions. Since the dipole moment of the water molecule 
can be resolved into components along the hydrogen bonds, it follows 
that the only water molecules whose moments contribute to the capacity 
lies between the outer Helmholtz plane and the mercury surface. This 
is shown in fig. 3, in which the arrows represent dipole moments. 
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Quite apart from the structure of the water, if the water molecules on 
the solution side of the outer Helmholtz plane have any average dipole 
moment, the ions, which are free to move there, will arrange themselves 
So as to screen out the field, in a reasonably concentrated solution, and the 
dipole moment in this region will not be observed. However, the fact 
that the structure of the water is maintained may help to explain why the 
capacity of the diffuse layer outside the outer Helmholtz plane can be 
calculated accurately with the assumption that the water has its normal 
dielectric constant (Grahame 1954). 


Fig. 3 
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The rate at which the water molecules turn over is of the order of 
magnitude of the reciprocal of the dielectric relaxation time. This is 
much faster than the rate at which the ions diffuse, so that in calculating 
the equilibrium configuration of the water we shall assume the ions to be 
frozen into fixed positions. 

The water molecules whose dipole moments contribute to the capacity, 
between the outer Helmholtz plane and the metal, are probably in contact 
with the metal. For the barrier to movement of an ion towards the metal 
is presumably the high energy needed to remove its primary hydration 
shells, so the distance between the outer Helmholtz plane and the metal 
is probably only the diameter of a water molecule plus half that of an ion, 
and the distance of penetration of the field into the metal. Indeed the 
outer Helmholtz plane may even pass through the first layer of water 
molecules, and if this happens, their contribution to the capacity will be 
‘reduced. 

Kemball (1947) has shown that the heat of adsorption of water molecules 
from the vapour on to a mercury surface is 17-6 kcal/mole. This is so much 
higher then the energy of a hydrogen bond between two water molecules. 
that it seems likely that the water molecules in contact with the mercury 
in an inter-face are in much the same state as if they were adsorbed in a 
close-packed mono-layer on mercury, and are only slightly affected by the 
presence of the rest of the water. This is supported by the fact that the 
interfacial tension of a mercury-water interface (427 dynes/cm) is almost 
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the same as the surface tension of mercury with a monolayer of water on it 
(438 dynes/cm). 


4.3. Hlectrostatics 


To calculate the contribution of these water molecules to the capacity 
we need to know at what angle their dipole moments are held by the 
mercury, whether they have a preferential orientation corresponding to the 
‘natural field’ on the metal surface discussed by Grahame (1957a), and 
to what extent their orientations are correlated with those of other water 
molecules (giving a value of w* for them)t+. None of these questions can 
be answered with any certainty, but we can attempt to answer them on 
the following lines. . 

If the surface of the mercury were an ideal electrical conductor a water 
molecule would be attracted by its electrical image in it, and would there- 
fore form a bond with the mercury very similar to that with another 
water molecule. The surface is not an ideal conductor, but a measure of 
its deviation from one is given by the ratio of the screening distance of an 
electric field in the metal to the distance of the charges we are interested in 
from it. In the case of a water molecule this ratio is less than one for most 
of the charge, so that to treat the mercury as a conductor is reasonable. 

The non-zero screening distance of the field in the metal means that 
image charges are further away than they should be, and image forces 
correspondingly weaker. Therefore if the interaction of a water molecule 
with the mercury is akin to hydrogen-bond formation, we expect that 
three of its bond-forming directions are attached to it, to account for the 
heat of adsorption. This would leave one bonding direction free to link 
with another water molecule. In the case of water adsorbed on mercury, 
Kemball shows that the entropy of adsorption suggests the assumption 
that the water molecules are associated in pairs. Our picture accounts for 
this in terms of bonds, though it is difficult to vizualize geometrically. 

The attachment of three bond-forming directions to the mercury may 
be aided by the atomic structure of the mercury surface, assumed close- 
packed. In such a configuration the dipole moment of the molecule makes 
an angle cos“! (1/4/3) with the normal to the surface, since the bond- 
forming directions of the molecule are directed to the vertices of an 
approximately regular tetrahedron. Just two configurations of the water 
molecule can be realised in this way. The components of its dipole 
moment along the normal to the interface are equal and opposite in them, 
and they are drawn in fig. 4. 

The water molecules will have a preferred sense of orientation on the 
mercury surface, if they bond to a mercury atom with a lower energy 
through a proton (+) bond than through a lone-pair (—) bond, or vice 
versa. If we regard the metal as a conductor, and the bond as due to 
electrostatic forces between rigid charged structures, we can determine 
a a 

} Iam grateful to Dr. A. D. Buckingham for discussion of these questions, 
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this electrostatically, in a very similar manner to that used by Buckingham 
(1957) to explain the difference in hydration energies of positive and 
negative ions of similar sizes. 


Fig. 4 
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We define the dipole and quadrupole moments of a system containing 
charges e (not necessarily all equal) at positions r= (aj, a, 23) to be 


u= Der 


6,5) UA, 

respectively, where the summation is over the charges. We consider two 
systems of charges, each electrically neutral as a whole, having dipole 
and quadrupole moments p, p.’ and 6,4, 0,,’.. The unit vector from the 
origin of the undashed system to the origin of the dashed system is 
A=(A,, As, Az), and the distance between them is R. This is supposed to 
be large compared with the spread of either system about its origin. The 
interaction energy between the two systems is 


W= — Hype (BA Ag —Sqg)/K RF 
— (tug9 gy! — ba’ 9 y)(BAx3 py + BAgDya+ BA,Sap— LBAAgA,)/2K Rt 
+ 0(1/R°), ee PORE a ac Sa Legit aire Hine GR OR ae (ah 


if the charges are in a medium whose dielectric constant is K. If the 
dashed system is the electrical image of the undashed one in a plane con- 
ductor, the energy of the image-force is half W: this can be seen either by 
considering moving the system into position from an infinite distance 
from the conductor, or from the fact that W is the part of the field energy 
due to interaction between every pair of charges, one in each system. If 
the conductor is replaced by the image charges, the contributions to the 
field energy are clearly symmetrical about the plane of the conductor, but 


the field only really exists on one side of it. 


and 
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2 is now the unit normal to the plane of the conductor, and for every 
charge e at r in the undashed system, there is a charge —eatr—2(A.r)aA 
in the dashed system, relative to the image of the origin of the undashed 
system. Therefore 

p= —pt2(A.u)a 
and 
Bg = —9 apt 20 pAAg + 20 yAAg — 46, AAsd Az: 


a. yByo arty 
Substituting these expressions in (3) we obtain 
LW = —[p?+(A.p)?]/2K Be 
—3[A.wlOaparp— aq) +2 aphabip]/2KR 
+0(1/R?). 


If the system of charges is a water molecule, its dipole moment is along 
one principal axis of its quadrupole tensor. We take this as the x-axis, 
and the x,- and x,-axes along the other two principal directions of the 
quadrupole tensor, with the «,-axis parallel to the line joining the protons. 
Writing 


w= (0, 0, HL), 
[0,6] = diag (9;, 6, 63), 
and 
A; = cos 0, 
we obtain 


4W = —p?(1+ cos? @)/2K R8 
— 3y cos 8[ (83 — 44)(1 —A4?) + (83 — 8y)(1 — Ag) ]/2K RA 
+0(1/R>). swith Ne A ge ce RE ge, 


An accurate theory would also take into account the polarization of the 
molecule by its image. 

Given accurate values of the components of the quadrupole tensor, we 
could use (4) to calculate the equilibrium orientation of the molecule. 
However here we shall assume it has one of the orientations shown in fig. 4, 
with three bonds in the direction of the metal (or more exactly, one away 
from it), and use (4) to calculate the relative energies and hence in 
thermal equilibrium the probabilities of the two orientations. In the 
discussion of both these points it should be borne in mind that the series 
for the interaction energy does not converge rapidly when the charges are 
as near the conductor as those we are considering, compared with their 
average spread. However, the next term in the energy comprising the 
dipole-octupole and quadrupole—quadrupole interaction, is complicated, 
and it is of the same order as the van der Waals forces, so to get a better 
approximation they would have to be brought in as well. 
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With the axes used in (4), the components of A are (+ V/ 2,0, — 4/4) 
when a proton (+) in the water molecule points away from the metal, 
and (0, + /3, \/4 when a lone-pair orbital (—) does. The components 
of the dipole moment along the outward normal to the metal are + pl/3 
and —u/1/3 respectively. It follows from (4) that the difference between 
the interaction energies in the two positions is G= 2u[205 — 0, — 6,]/4/3K RS, 
the position with the dipole pointing towards the metal having the lower 
energy if the quantity in square brackets is positive. The first term in (4), 
the dipole-dipole energy, is — 22/3.K R3 in both positions. 

The value of # is uncertain, but it can hardly be less than 34. The 
dipole-dipole energy multiplied by K is then at most 1-2 kcal/mole. 
Buckingham derives an empirical value of 26,—6,—6, for water, namely 
39x 10% e.s.u. Theoretical estimates (Duncan and Pople 1953) are an 
order of magnitude smaller. Buckingham’s value gives K@ as at most 
1-5keal/mole. These values are so small compared with the energy of a 
hydrogen bond, let alone the heat of adsorption of water on to mercury, 
that we are almost compelled to suppose that K =1 in this interaction, 
and we shall do so. The smallness of the dipole-dipole interaction shows 
that our assumption, that in equilibrium the water dipole is not along 
the normal to the metal surface, is not an absurd one. 

The dielectric constant of normal water is so high, that in calculating the 
equilibrium configuration of the water molecules in the double layer, we 
can consider the potential difference across it to be applied between two 
planes, in the metal and in the water. Outside them the material can be 
treated as a conductor to a high degree of approximation. Kirkwood 
(1939) makes a similar approximation in supposing that the field in his 
spherical cavity is 3E/2 instead of 3KE/(2K + Ky). 

We must now refer to fig. 3 again. If each of the water molecules in 
contact with the mercury has three bonding directions occupied, the 
remaining one is bonded to another water molecule. The component of 
its dipole moment in this direction is repeated along a bonding direction 
in the next molecule, and so on until the chain of dipoles ends on an ion, 
or on the surface of the metal again. For the purpose of calculating its 
interaction with the field, however, this chain can be considered to run 
from the plane of the conducting surface of the metal to the plane beyond 
which the water is normal, provided the ions are not too close together. If 
a typical normal to these planes cuts water molecules between them, the 
dipole moment of this chain will have a component .* = + m|+/3 along the 
normal to these planes. The potential drop between them is ¢, the 
potential drop from the metal to the solution, since the material outside 
them is supposed to be conducting. ‘The distance between them is 
approximately equal to nd, so the electric field H is p|nd, a he valor 
action energy of the chain of dipole with the field is —Ly*= ¥ ¥o/ / 3d. 
It is useful that this is independent of n, since it would be very difficult to 
estimate the minimum distance between the planes. However, the 
independence is not really surprising : since the criterion of a suitable 
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plane on the solution side is that the water should be normal beyond 
it, any plane further into the solution should do equally well, as long as 
the distance between the planes remains small compared with the distance 
between the ions in the outer Helmholtz plane. 


§ 5. A CALCULATION OF THE CAPACITY 


We can now make a theoretical estimate of the capacity, taking into 
account the water dipoles and the adsorption of ions. We shall find that 
with a reasonable choice of some parameters, they explain the hump in the 
middle of capacity curves, and the rise on the anodic side, respectively, and 
we shall try to interpret the necessary values of the parameters. 

If there are N free water molecules per unit area of the interface, of 
which L have a component of their dipole moment from the metal to the 
solution, and M the other way, then the energies of the chains of dipoles 
associated with these molecules are 

A=0'—pd)/3d 
and 

B=D' +pd/+/3d 
respectively: C’ and D’ are the adsorption energies of the water molecule 
in contact with the metal and are equal to half the relevant values of W. 
Provided the layer across which the chains of dipoles run is not thick, the 
chains will be approximately independent of one another, so that 
L/M=exp{—(A—B)/kT}. Therefore the contribution of the water 
molecules to ¢ is 


— 4n(L—M)p/>/3Ky = (47N yu /+/ 3K) tanh {(A —B)/2kT! 
= — (4nN yp] /3K,) tanh {u($— V)/V/3dkT} 
where 
V=/3d(C’ —D')/2u=+/3dGj2p.. . . . . . (6) 
If there is a surface charge o on the mercury balancing a charge —o on 
the outer Helmholtz plane, it contributes 47od/K,) to ¢. Adding the 


contributions (1) from the adsorbed complexes, and the contribution from 
the water dipoles, we obtain 


p=4rod/K, 
— (47 N 0" /K geo”) EXP {H44($— Vo) /dkT} 
— (407Np/1/3Ko) tanh {u(d— V)/4/3dkT}. U6) 
In this equation we have put 


W othfity Vg. hs ee Mea 


and we have assumed that the number of adsorption sites per unit area 


for the complexes is the same as NV, the number of water molecules. It 
follows that 


(on Kod/4nd 
+ (Nyye;"/¢o"d) exp {Hi(P — Vo)/ dkT} 
+ (Nyu//3d) tanh {u(¢ — V)//3dkT}. Aiea, (3) 
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The differential capacity C of the layer is da/0d, namely 


C=K,/4rd 
+ (Npy?ey” /cg"PET) exp {u4($— Vo) /dkT? 
+ (Np? 30k) sech? {u(d—V)/+/3dkT%. 5 alte ERED 


In this expression, the first term is the ordinary capacity per unit area 
of a parallel-plate condenser, the next is due to the adsorption of complexes, 
and the last is due to the water. The exponential in the term due to 
adsorption rises extremely rapidly when ¢ is greater than Vy. We shall 
calculate a capacity curve for mercury in sodium fluoride in water from 
(9). The values of the unknown parameters in (9) will be adjusted to fit 
the experimental curve given in fig. 2. 

The following values have been taken in the calculation, and the resulting 
capacity curves at 0°c and at 85°C are shown in fig. 5. 


N = 10? molecules/em!, w= 1-84 x 10-18 e.s.u.; V=0; 
V o=0-7volt; »,=3-31 x 10-8e.s.u.; d=3-454; K,=2. 


The first two values. are standard. Auty and Cole (1952) have found 
that the dielectric constant of ice tends to the value 3-1 as the frequency 
of an incident alternating field is increased above 50 ke/see. This seems to 
be the value corresponding to libration and distortion of the water mole- 
cules fixed in definite orientations (Onsager and Dupuis 1960) and is 
appropriate here. However, part of the thickness of the condenser in 
the interface is the depth of penetration of the field into the metal, where 
the dielectric constant is probably lower, so we have taken K,=2. It 
follows from the calculation in the Appendix that to take an average 
value of K, for the whole layer is not strictly correct. d is adjusted to 
fit fig. 2: its value is reasonable. 

V can be calculated from (5), if we know the quadrupole moment of 
water. The value of 1-5kcal/mole for G gives V as 0-5l volt. If. this 
were correct, the hump due to water molecules would be about half a volt 
to the anodic side of the electrocapillary maximum, in the absence of 
specific adsorption of ions. However, Buckingham’s value of the quad- 
rupole moment is doubtful, and is probably too high (Everett 1957 et 
is noticeable that the hump on most experimental capacity curves is about 
0-1 volt to the anodic side of the electrocapillary maximum. If our inter- 
pretation of the hump is correct, this could be used to find an empirical 
value of the quadrupole moment of water, which would be of the same 
order of magnitude as Duncan and Pople’s. Their value would give V 
as 0:06volt. In the calculation V has been taken as zero for simplicity. 

If we assume that mercury ad-atoms are formed, and that the expo- 
nential term in (9) is due to them, p,/d could be measured by a technique 
due to Gerischer (1958). This consists in estimating the variation in the 
number of ad-atoms in equilibrium on the surface, when the concentration 
of metal ions in the solution is changed. Applied to silver, this method 
yields a value of 0-19 of a proton charge ¢ for u,/d, which we may call ye 


K 
P.M. 
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to conform with Gerischer. It would be interesting to apply this tech- 
nique to mercury: presumably the formation of an ad-atom is a step in 
detaching a mercurous ion from its surface, but since these ions are asso- 
ciated in pairs in solution they are not really comparable with silver ions. 
For the calculation, we have taken y=0-2, giving the value of fy, and 
adjusted V, to fit fig. 2. This fitting can be performed with quite a range 
of values of y, so we should not attach much significance to the value of 
V,, especially as Gerischer’s method of estimating y 1s very rough. 
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W, the heat of formation of an ad-atom, can also be measured, as it is. 
the difference between the heat of dissolution of the metal and the heat of 
dissolution of an ad-atom. ‘The latter can be found from the dependence 
of the number of ad-atoms on temperature (Gerischer 1958). However, 
its connection with V depends on y (eqn. (7)). Very roughly we should. 
expect the heat of formation of an ad-atom on a mercury surface in contact 
with its vapour to be half the latent heat of vapourization of mercury. 
The heat of formation of an ad-atom at a mercury-solution interface should 
be lower, owing to the interaction of its dipole moment with the charges 
in the solution. Our value of V, gives W,=0-l4ev, and this may be. 
compared with the latent heat of evaporation of mercury, which is about 
0-6ev. It is therefore quite reasonable to suppose that the exponential 
term in the capacity is due to ad-atoms. 

An alternative possibility, in a fluoride solution, is that it is due. 
to the adsorption of hydroxyl ions+. The test of this is to measure capa- 


cities in solutions of different px. The contribution of hydroxyl ions to. 
C may be written from (9), 


(Nypy"/@°kT) exp {u4(6 — V)/dkT} 
} This suggestion is due to Dr. J. Wanklyn. 
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where ; 
Vy =Vot(madkT’ 114) In (Wp, ¢9). 

We must remember that 1, is the dipole moment of an adsorbed hydroxyl] 
ion and its induced charge on the outer Helmholtz plane. If we assume 
that in forming a hydroxyl ion from a water molecule in contact with the 
metal, a proton moves through about half the double layer, then p,/d =e 
(or less if the charge on the hydroxyl ion is partly neutralised by its bond 
with the metal). £7’ is about 1/40ev, so that an increase of one unit in 
pH should shift the part of the capacity curve due to hydroxyl ions by at 
least 0-12 volt to the cathodic side, even if m=1. Grahame (1947) has 
measured capacities in normal solutions of sodium hydroxide and per- 
chloric acid. His values of C are plotted against potential measured 
relative to the electrocapillary maximum, whose potential on our scale 
varies with pH if the hydroxyl ion is adsorbed. However, the hump on 
the curves gives a fixed reference point if it is due to the water. A com- 
parison of these curves with that for sodium fluoride shows that the rise: 
on the anodic branch of Grahame’s curves does shift to the cathodic side 
with increasing pH, but only about a third as much as if it were simply due 
to the adsorption of hydroxylions. The very considerable shift in the curve 
for sodium hydroxide suggests that there the rise is due to adsorption of 
hydroxyl ions. The rise in the other cases is more likely to be due to the 
formation of ad-atoms of mercury, or to the formation of complexes on 
the surface between mercury atoms and hydroxyl ions. A careful series 
of capacity measurements in fluoride solutions having different values of 
pH would be very valuable. Acid solutions would have to be contained 
in material unattacked by hydrofluoric acid. 


§ 6. Discussion 


The calculated values of C shown in fig. 5 are in fair agreement with the 
experimental ones for sodium fluoride (Grahame 1957 a) shown in fig. 2, 
except on the far cathodic side. The calculated temperature-dependence- 
is not large enough on the anodic side, and the fact that the hump on the 
experimental curves is always at a small positive value of the surface 
charge, suggests that water molecules tend to have their dipole moments 
pointing towards rather than away from the mercury. 

Adsorbed anions give a very similar contribution to C to that of ad- 
atoms, so we may also compare fig. 5 with many other capacity curves for 
electrolytes whose anions are specifically adsorbed, for example sodium 
nitrate (Grahame 1957a). The similarity between the curves is remark- 
able, and is very strong evidence that an adsorption process also oceurs in 
sodium fluoride solutions. We suggest that capacity measurements with 
mercury amalgams or other liquid metals should be made to alter the 
adsorption energy of ad-atoms, in order to check the hypothesis that they 
contribute to the capacity. If anions are very easily adsorbed the capa- 
city due to their adsorption may rise rapidly when the polarization is only 


K2 
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slightly anodic and this may mask the capacity due to the water molecules 
turning over. This may be the reason for the absence of humps in curves 
for sodium hydroxide, carbonate and iodide (Grahame 1947). 

When humps appear in experimental capacity curves for aqueous 
solutions they generally have much the same height. This supports the 
hypothesis that it is the water which causes them. An important test of 
it is to measure capacity curves in other solvents, or in fused salts. 
Grahame (1955) and Mine and Jastrzebska (1960) measured capacity 
curves for sodium and potassium fluorides dissolved in methyl and ethyl 
alcohol, and these do not show humps. This may be because the alcohol 
molecules have a strong preferential orientation on mercury, so that they 
do not turn over in the usable range of polarization, or because the freezing 
points of the alcohols are lower and the measurements in them should be 
compared with those at higher temperatures in watery. Measurements 
of the capacity in methyl and ethyl alcohols over a range of temperatures 
should decide between these explanations. 

An interesting feature of some capacity curves in water is that speci- 
fically adsorbed anions sometimes seem to intensify the hump rather than 
masking it: for example the curve for potassium nitrate may be compared 
with those for sodium fluoride (Grahame 1957a). A possible explanation 
of this is that the outer Helmholtz plane is further from the metal in a 
nitrate solution and therefore the layer of water contributing to the 
capacity is thicker{. The nitrate ion is so strongly adsorbed when the 
polarization is anodic that on the average there are many more cations 
than anions in the outer Helmholtz plane, or rather in the diffuse layer 
(Grahame and Soderberg 1954). However, the cations will be clustered 
round the specifically adsorbed anions: this is shown in fig. 6. The ions 
facing the rest of the metal surface will mainly be anions. If they are 
nitrate ions they will be further from the metal than fluoride ions, owing 
to their larger size and to the greater number of water molecules in their 
primary hydration shells. The water molecules contributing to the 
capacity will be in the thicker part of the layer, away from the adsorbed 
ions. 


6.1. Compression of the Double Layer 

So far we have not considered the fact that all experimental capacity 
curves rise on the cathodic side. 1t has been suggested by Macdonald 
(1954) that this is due to a compression of the layer in the strong electric 
field. This is supported by the fact that Grahame’s curves for methyl 
alcohol solutions rise roughly symmetrically on both sides, and more 
sharply than those for water solutions on the cathodic side. The layer is 
presumably more compressible in methyl alcohol solutions than in water, 
as Macdonald points out. 


} This suggestion is due to Dr. J. N. Agar. 
¢ Our assumption that just one layer contributes is approximate, and the 
contribution is really proportional to the thickness of the layer. 
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We can obtain a rough estimate of the elastic force, necessary to produce 
the capacity rise observed on cathodically polarized mercury in an aqueous 
sodium fluoride solution, as follows. We suppose that only the first term 
in (9) is important there. This is not a very good approximation with 
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the values of the constants we have used, as it can be seen from fig. 5 that 
the capacity due to the water, the third term in (9), is still decreasing 
throughout the range of measurements on the cathodic side. The thick- 
ness of the layer is d—eH/x, where x is the force necessary to produce a 
decrease of one unit in the thickness of the layer, and 
f = 4ro0(d — eH /x)/K,+constant, 
from which it follows that 
C= K,/4n(d — 87e0/K 9k). 
If the second term in this denominator is small, that is, the layer is not 
much compressed, 
dC |do = 87reC | K kd. 

Using the graph reproduced as fig. 2, dC/do is about § volt~? when C is 

20 uF/em?. With the previous values of K, and d we obtain 
Ke 2ev/A%. 

The interest of this is that it gives a clue to the height of the free energy 
barrier which an ion which is specifically adsorbed, or deposited on the 
metal, has to cross to become free of part of its hydration shell. Of course 
alkali metal ions, which are the ones whose displacement we are estimating, 
are not specifically adsorbed, though in principle they have a low energy 
state on the metal side of the interface, when they form an amalgam with 
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the mercury. If we imagine that there is such a state, and that the ions init 
are at a distance J from the outer Helmholtz plane, then when they are in 
equilibrium with the hydrated ions there, the free energy of an ion in an 
intermediate state at a distance z from the other Helmholtz plane will be 
something like 
F=1A[1—cos (272/l)]=A sin? (72/1), 

where 4 is the height of the energy barrier between the two states. But 
we know that for small values of z, 

F=hx2", 
so that 

A =k? /27. 
The value of J is uncertain, though for a specifically adsorbed ion in the 
inner Helmholtz plane it is given by 

lL=d—d’. 
Here d’ is the distance of the inner Helmholtz plane from the conducting 
surface of the metal. Since a hydrated ion is probably surrounded by a 
close-packed shell of water molecules, we may expect that / is roughly the 
diameter of a water molecule. Taking this as 2-764, the energy barrier 
height is about 0-8ev. It is difficult to say whether this is a reasonable 
value. Vermilyea (1956) shows that an activation energy of this amount 
in the electrolytic deposition of a metal would make the reaction have a 
fairly high overvoltage (which he attributes to a different mechanism), 
but the value of the over-voltage is very sensitive indeed to small changes 
in the activation energy. 

The compression of the double layer at high field strength should 
certainly be studied further, on account of the light it can throw on electro- 
deposition. If it occurs it has the further consequences that the thickness 
of the water layer contributing to the capacity varies, and that d in our 
theory of the capacity is a function of 6. To take this into account would 
complicate the calculation without much improving its accuracy. 
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APPENDIX 
AVERAGE PoTENTIALS 


We set up an orthogonal Cartesian coordinate system with coordinates 
©, Y, 2, with the (x, y) plane parallel to the double layer, and we consider a 
disc-shaped right cylinder in the double layer whose generators pass 
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through a circle C in the (x, y) plane. Its interior we call S, of area A, and 
the disc is bounded by planes z =u and z=», which are cut by the generators 
in circles whose interiors are U, V respectively. These are the pro- 
jections of S on to the planes. The construction is drawn in fig. 7 


Fig. 7 


If the electrostatic potential is ¢, and the field and displacement are E 
and D respectively, generated by a volume charge density p and a polari- 
zation P, then 

D=E+47P= —ogradd+4nP. See ene CAL) 
Applying Gauss’s theorem to the disc whose interior is J and unit outward 
normal n, we obtain 


{| D.ndS+ | | D-ndS=4n|| | eee een 
U V Mi 


provided that the radius of C is so much larger than the thickness of the 
disc (v—w) that we can neglect the flux through its curved surface com- 
pared with that through the flat ones. Substituting for D from (A 1) into 
(A 2) gives 


I, Iz Yu) 4P Avs Y; u) — (FE) y, yt 4rP,(x, Yy, »| dx dy 
=A | dw {| p(x, y, w)da dy theres) 
u Ss 


where P, is the z-component of P. Denoting averages over the projections 
of S on to planes es to the (x, y) plane by bars, for example 


B= 7 [| eu )da dy, 


and denoting derivatives with respect to z by dashes, (A 3) becomes 
b' (u) — 4 P(u) — 6’ (v) +40 P(r) 


=4n[" ployde RM ayaa ee (AS) 


u 
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We now take the plane z=w to the left of all the charges in the layer. 
Subject to this condition the integral on the right of (A 4) is independent 
of u, and p,(w)=0. It follows that ¢'(u) is independent of wu, so it is a 
constant, and if the layer is electrically neutral as a whole, and no field is 
imposed on it externally, it must vanish. Of course we must keep wu so 
that 


ljv—u|</A. 
to apply this equation. On integrating it with respect to v from u to z 
and rearranging it, we obtain 


d(u) — b(z) = —47 i P(v)dvt+ cal ae p(w)dw 


=4n | “[— P,(v) + (2—vplv) dv. 


This is exactly the same potential drop as if the charge were smeared 
out in planes parallel to the double layer, for then the part between the 
planes z=v, z=v+dv could be regarded as a combination of an electrical 
double layer of moment P,(v)dv, and a surface charge of density p(v)dv at 
a distance (z—v) from the field point. In particular, if part of the polari- 
zation of the medium to the right of z=w is due to a perfect dielectric of 
constant Ky, then 


P= —(K,—1)grad¢/47+P, (say) 
and. 


7 5 4n [* 


bu) -$@) =F | [-A,(e) + —v))] dv. 


Equation (6) is of this form; the potential drop ¢ in it is between two 
conductors, so averaging ¢ parallel to the double layer makes no difference. 
A special case of this theorem was proved differently by Grahame (1958). 
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CORRESPONDENCE 


Magnetic Properties of Chromium between 0° and 350°c 


By E. W. Cottrnest, F. T. Hepecock and A. SID DIQrt 


Department of Physics, University of Ottawa, 
Ottawa, Canada 


[Received August 29, 1960] 


‘THE maximum in the susceptibility of chromium observed by Bates and 
Baqi (1936) was attributed by McGuire and Kriessman (1952) to the 
presence of the antiferromagnetic oxide Cr,O, in their samples. Fine 
et al. (1951) and Séchtig (1940) reported no abnormality in the susceptibility 
of chromium at 40°c. Recently Lingelback (1958) observed a maximum 
in the susceptibility of high purity Australian chromium near 37°c in 
a sample where the content of Cr,O; must have been extremely low. 
Neutron diffraction investigations of the possible antiferromagnetic 
transition (see Edwards 1959) have led to confusing results in that the 
diffraction pattern seems extremely sensitive to particle size and purity. 
However recent specific heat measurements by Morrison ef al. (1960) 
indicate that there is an antiferromagnetic transition occurring in 
chromium at 38-5°c. The purpose of the present note is to report 
measurements on the magnetic susceptibility of chromium made on samples 
of various particle size and purity and electron spin resonance measure- 
ments in a sample of high purity chromium. 

Figure 1 shows the susceptibility of chromium as a function of tem- 
perature for the various samples studied. The description of the samples 
and room temperature susceptibility value is included as a footnote in 
the figure. The susceptibility of all of the samples was measured as a 
function of field to a maximum field of 8000 oersted. All of the samples 
showed a field dependent susceptibility and the ferromagnetic correction 
Oc, obtained from H— plots is also included in fig. 1. Asexpected the higher 
purity samples showed a smaller ferromagnetic correction; the correction 
for all samples obtaining a maximum value in the region of the susceptibility 
maximum. All of the samples show a maximum in the magnetic 
susceptibility at a temperature of 39+0-5°c. It seems fairly conclusive 
therefore that the maximum in the susceptibility of chromium is inde- 
pendent of the physical nature and purity of the sample in the range of 
samples studied. In particular the apparent influence of nitrogen content 
on some properties of chromium (see Edwards 1959) ee not influence 
the size or position of the susceptibility maximum at 39°C. 
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The absolute susceptibilities listed in fig. 1 are in reasonable agreement 
with the value of 3:17 x 10-* x 10-e.m.u./g reported by Lingelback (1958). 
The high-temperature slope of the susceptibility versus temperature curve 
is in agreement with Lingelback ; however the slope on the low-temperature 
side of the susceptibility maximum has the opposite sign but is in agreement. 
with the earlier observations of Bates and Baqi (1936). 
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Magnetic susceptibility of chromium and the ferromagnetic correction Co, 
plotted as a function of temperature 


x x10; : ; 
Sample A 0-01 a sacle) N, (weight °%) O, (weight %) 
(1) 3-70 0-007 0-02 
(2) 3-41 0-0230 Estimated 0:5 
(3) 3:08 0-0004 Estimated 0-01 
(4) 3-10 


No analysis made but reported 
higher purity than sample 3 


Chromium sesquioxide is known to be antiferromagnetic with a Neél 
temperature of 38-5°c. Therefore one cannot ignore the possibility that. 
all of the chromium samples contained sufficient oxygen to form a fixed 
concentration of Cr,O; even though the analysis of the samples would 
indicate too low an oxygen content to form sufficient Cr,0, to explain 
the. observed size of the susceptibility bump. The amount of Cr,0, needed 
was estimated from the susceptibility data of Trounson et al. (1950) to be 
4 weight per cent. The electron spin resonance signal from Cr,O, has the 
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characteristic behaviour of an antiferromagnetic resonance in that it 
abruptly diminishes in amplitude near the Neél temperature (‘Trounson 
et al. 1950). An attempt was made to detect a Cr,O, signal in a thin strip 
of freshly etched high purity A.R.L. chromium. Although some absorption 
took place it was of sufticiently low level that it would have to be classified 
as noise (curve A in fig. 2). When 46 micrograms Cr,O3 were placed in 
the cavity along with the metallic chromium (the chromium strip was 
retained in order to keep the cavity Q and hence sensitivity constant) a 
distinct increase in resonance absorption was observed above 33°c (curve 
Bin fig. 2). Allowing for the skin effect in the chromium metal the relative 
weights of Cr,O, and metallic chromium ‘seen’ by the r-f. field was 4%. 
A visual estimate of the sensitivity derived from the results in fig. 2 shows 
that certainly the transition in less than 1°/ Cr,O3 in metallic chromium 
would have been detectable in the background level of the pure metal 
(curve A). Itis therefore apparent that there is not sufficient Cr,O, present 
in the purest chromium to give the susceptibility maximum at 39°c and 
the resonance results support the previous conclusion based on the suscep- 
tibility results alone that there is an antiferromagnetic transition taking 
place in chromium itself at 39 + 0-5°c. 


ee 3 
Fig. 2 


A = Chromium and Equiv. of 4 wot. % Cr, 03. 
B = Etched A.R.L. Chromium Strip. 
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The Mechanism of Surface Effects in Crystal Plasticity 
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THE stress-strain curves of many plastic monocrystals are sensitive to 
conditions at the crystal surfaces. Thin surface films of solid materials 
often harden crystals (the Roscoe Effect, Metzger and Read 1958), and 
surface-active liquids tend to cause softening (the Rehbinder Effect, 
Likhtman et al. 1958). The effects appear schematically as in figs: 
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Schematic effects of surface conditions on the stress-strain curves of plastic 
crystals. 


This figure illustrates a prime feature of the effects which is not explained 
convincingly by previous interpretations. This feature is the fact that. 
the effects do not appear full-blown for small plastic strains ; instead, 
rather substantial strains (> 10~‘) are required before they acquire their 
full magnitude. This means that the effects are not primarily yielding 
phenomena, but depend on conditions that develop sometime after 
yielding has occurred. Therefore, a mechanism that simply involves the 
blocking of dislocations by a surface film does not interpret the observa- 
tions satisfactorily. Something more subtle is required, and I believe 
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that insight gained in recent experiments with LiF crystals (Johnston 
and Gilman 1960) suggests a suitable mechanism. Moreover the mech- 
anism is amenable to direct experimental study using etching methods 
to observe the effects of surface environment on the rate of surface 
cross-glide. 

The interpretation that is proposed is illustrated in fig. 2. At (@) in 
the figure, a screw dislocation, AB, lies on a glide plane that is inclined 
to the surface at an angle, 6. The dislocation’s Burgers vector, b, is 
indicated and the dislocation tends to move in the direction AC. However, 
the dislocation line can shorten its length by moving off the primary 
glide plane CDEF as shown at (6) in fig. 1. This ‘ cross-gliding © motion 
is opposed by the fact that a surface step is created by it ; and, in an 
anisotropic crystal, by the fact that the motion may have to occur on 
a plane other than a primary glide plane. Such cross-glide motions 
have been observed in LiF crystals by Johnston and Gilman (1960). 
Since the cross-gliding is opposed by the surface energy of the step that 
is left behind, the distance over which gliding occurs will increase with 
decreasing surface energy and conversely. A surface film that is harder 
than a crystal, on the other hand, will mechanically resist the formation 
of the step and thus restrict the motion of the screw dislocation. 


(a) 


Cross-gliding of a screw dislocation near the surface of a crystal : 
(a) initial configuration ; (6) after some cross-gliding. 


Cross-gliding of screw dislocations at the surface of a crystal has an 
important effect on the overall plastic behaviour because it leads to 
dislocation multiplication (Johnston and Gilman 1960). This in turn 
leads to strain-hardening because the greater the dislocation density in 
a crystal, the harder it becomes. Since this surface mechanism does not 
involve the initial motions of dislocations on their primary glide planes 
it has the required characteristic of not influencing the initial yield stress, 
but only influencing the progress of a stress-strain curve. . 

A feature of the proposed mechanism that makes it especially attractive 
as a surface mechanism is that the force that is resisted by the surface 
step is presented by a dislocation line with relatively low line-energy. 
This occurs because the shear stresses around a screw dislocation that 
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intersects a surface are relaxed at the surface so the dislocation line has 
no elastic strain energy ; it has only the much smaller energy of its core. 
The core energy of a dislocation line is comparable in magnitude with the 
energy of a surface step, so it is reasonable to expect changes in the 
surface energy to affect the position of the end of the dislocation line. 
On the other hand, if one tries to account for surface effects in terms of 
the blocking of edge dislocations that are trying to emerge from a surface, 
then the surface energy must oppose a dislocation that has elastic as well 
as core energy. In this case it seems unlikely that small changes of the 
surface energy would have much effect. 

The fact that an optimum strain-rate is observed for surface effects 
can be accounted for because the cross-gliding has been found to occur 
at a certain moderate rate. This will be too small to cause a surface 
effect at high rates of plastic deformation. At deformation rates that 
are lower than the optimum, the distance of cross-gliding will be so 
large that little strain-hardening will occur. 

Glide plane orientation effects, irregular glide bands on crystals 
deformed in a surface-active medium, and enhanced distortions in the 
surface layers of crystals deformed with surface films present are all 
observations that seem consistent with the proposed mechanism. 
Moreover the mechanism is amenable to direct experimental study using 
etching methods to observe the effects of surface environment on the 
rate of surface cross-glide. 
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Delayed Alpha Particles from '°N 


By R. E. Szczz, J. W. Oxtnzss and E. L. SprenKEL 
Aeronautical Research Laboratory, Dayton, Ohio, U.S.A. 


[Received October 15, 1960] 


Iv is energetically feasible for the *N B-decay to populate three known 
a-unstable states in *O: the 2~ state at 8-88 Mev, the 1- state at 9-58 Mev, 
and the 2* state at 9:84 Mev (Ajzenberg-Selove and Lauritsen 1959). We 
have measured the alpha particle spectrum resulting from the decay of 
*°N in order to obtain information pertaining to the decay through these 
states. 


1®N was prepared by bombarding nitrogen gas, enriched} to 96% N, 


with 2 Mey deuterons from the Aeronautical Research Laboratory van de 
Graaff generator. The radioactive gas from the target was allowed to 
diffuse into a small cell which was viewed by a solid state detector of p-type 
silicon with a shallow phosphorus diffusion surface layer. The running 
procedure was to bombard for 15 sec, de-energize the accelerator, and then 
record for 15sec. This procedure was repeated several hundred times. 

The silicon detector was calibrated with the aid of radioactive «-sources 
andamercury pulser. The ‘dead’ surface layer of the counter, which was 
measured to be about 0-5 microns, was determined by observing the shift 
in pulse height when the angle of incidence of alpha particles was changed. 
In order to reduce the electron counting efficiency the detector was 
operated at a low bias (~1V). 

A pulse height spectrum is shown in the figure. The energy scale is 
given for alpha particles incident on the face of the detector; corrections 
for the ‘dead’ layer have been taken into account. The energy loss in the 
gas was negligibly small, averaging about 30kev. It can be seen from the 
figure that only one alpha particle group was seen: a broad group centred 
‘at about 1-72Mev. We identify this group with the broad state centred 
at 9-58 Mev (Ajzenberg-Selove and Lauritsen 1959). 

The centre of the observed alpha particle spectrum lies about 100 kev 
below the energy corresponding to the centre of the state as determined by 
alpha particle scattering on C (Hill 1953). A distortion of the spectrum 
from the normal (i.e. Breit-Wigner) shape can be expected because of the 
energy dependence of the f-transition rate (roughly proportional to £,° 
where H, is the 8 end-point energy) and, in addition, the shape of the state 
can be expected to deviate from normal because of the coulomb barrier 


+ The authors are indebted to Dr. J. A. Jacobs of the State University of 
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present to the out-going alpha particles. These two effects tend to work 
in opposite directions with the B-decay favouring high energy electrons and 
therefore low energy alpha particles, while the coulomb barrier energy 
dependence implies a shortened lifetime for high energy alpha particles 
and therefore a broadening of the state on the high side. At the energies 
relevant here, the B-factor is the stronger and therefore the alpha spectrum 
is shifted towards lower energy. 
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Pulse height spectrum from solid state ionization chamber alpha detector 
viewing decaying ‘*N nuclei. The arrows represent the pulse heights 


corresponding to the centres of the various alpha unstable states in 160 
The solid line merely represents an eye smoothing of the data. 


The observed spectrum was also subject to instrumental distortion, the 
primary effect being a smearing of the true spectrum due to alpha particles 
traversing the counter dead layer at various angles. Because of this 
instrumental distortion, it was not possible to make an accurate analysis 
of the spectral shape. However, even without allowing for the energy 
resolution of the detecting system, we find the peak narrower than would 
be predicted for the state having the width ‘= 650 kev (c.m.) as given by 
Hill (1953). The data presented here implies [,=450+100kev. This 
discrepancy is not surprising in view of the difficulty in extracting a width 
in a reaction in which the transition rates can be expected to vary widel 
over the width of the state. s 

From the intensity of the alpha group we find log ft =7-2 + 0-2 with the 
B end-point energy taken as corresponding to the centre of the observed 
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alpha particle distribution. We note that the spin of the 1*N ground state 
is known to be 2~ (Ajzenberg-Selove and Lauritsen 1959) and therefore the 
B-transition to the 1- 9-58mev 160 state is allowed. The observed 
B-transition rate is therefore a few hundred times lower than what would be 
normally expected for an allowed transition. Elliot and Flowers (1957) 
indicate that the 9-58 Mev state should be due to a three nucleon excitation 
and it is gratifying to note that the B-decay to this state is indeed highly 
inhibited, thus supporting the Elliott and Flowers description of 160. 

The 8-88 Mev 2~ state is known to be populated by about 1% of the 
*®N decays (Alburger ef al. 1959). Break-up into #C + wis parity forbidden 
and therefore a search for alpha particles from this state is a test of the 
validity of parity conservation in nuclear interactions. As is shown in the 
figure, the pulse height region corresponding to alpha particles from the 
8-88 Mev state contains pulses from the electron background. However, 
as the state is sharp and heavily populated (some 2 x 108 decays from this 
state were viewed) a meaningful upper limit for the alpha width of this 
state could be deduced. It was determined that as many as about 300 
alpha particles emanating from the 8-88Mev state would have been 
detected. Hence, we find I'",/T,<1-3 x 10-+. 

The gamma width of the 8-88 Mev state has been calculated by Elliot 
and Flowers (1957) to be about 3x10-ev, from which we imply 
[,<3-9x10-‘ev. The expected I’, in the absence of a parity inhibition 
can be estimated by calculating the single particle width (allowing for the 
coulomb barrier) and taking the expected IT, as 1/20 the single particle 
width in accordance with the survey of Wilkinson (1957). This estimate 
predicts a parity allowed [’,=3 x 10% ev and, therefore, using the notation 
of Lee and Yang (1956) we find #*<1-3x10-'°, emphasizing that this 
upper limit can be taken only as a rather crude estimate. 

The 2+ assignment to the 9-84mev 160 state implies a first forbidden 
B-transition. No alpha group was seen corresponding to this state: the 
‘pump’ in the spectrum in the figure at about the right energy for this 
group is not statistically significant. The upper limit to the transition 
rate to the 9-84 state was found to be log ft> 7-4. This value is consistent 
with the 2+ assignment to the 9-84 Mev state. 
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Remark on the Seebeck Coefficient of Polar Semiconductorst 


By J. APPEL 
John Jay Hopkins Laboratory for Pure and Applied Science, _ 
General Atomic, Division of General Dynamics Corporation, San Diego, 
California 


[Received September 23, 1960] 


For a non-degenerate extrinsic semiconductor, the Seebeck coefficient is 
given by the Pisarenko formula (Joffé 1958) 


SE AO Uk wes We ee P(t) 


The — and + signs correspond to electron and hole conduction, respectively 
(—e=electron charge). The number b depends somewhat on the electronic 
band structure and on the pertinent scattering mechanism; usually 0 is 
small compared to ¢/k,7'=log(N,/n), where for an isotropic n-type 
conductor with a parabolic band V,=2(27merrk,T/h?)??, and n and mert 
are the concentration and the scalar effective mass of the negative charge 
carriers, respectively. The validity region of eqn. (1) is discussed else- 
where (Lautz 1953). 

The Seebeck coefficient depends on the effective mass mez via €. Here 
we want to consider the dependence of « on that mer: which is associated 
to the appropriate charge carriers in polar crystals, called polarons (Pekar 
1954), provided the polaron eigenvalues can be described by an effective 
mass concept. Then mer; can be quite different from the bare, free electron 
mass m, and also from the rigid lattice mass m*. The polaron effective 
mass Merp=™m** is defined by the low-lying energy eigenvalues of the 
polaron eigenstates. In order to calculate these eigenstates there exist 
two essentially different approaches. 

For the approach, resulting in what is called the large polaron picture, the 
paper of Frohlich et al. (1950) is characteristic. The theory is based on a 
model in which an excess electron interacts with a dielectric continuum, 
instead of a crystal lattice. In restricting to the ground state of the 
system, electron plus the dielectric continuum which is represented by 
plane polarization waves, the polaron eigenstates are determined by virtual 
emission and reabsorption of quantized longitudinal polarization waves 
(polarization phonons). Thus, the eigenstates of slow polarons at very 
low temperatures are calculated by assuming that the dimension of a 
polaron is large compared to the lattice constant. The effective mass m** 
of a large polaron consisting of an electron plus its associated cloud of 


+ Part of this project research is sponsored by the Advance Research Projects 
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virtual polarization phonons is, of course, temperature independent. The 
briefly described large polaron picture, based on the continuum model, 
does not apply in the case of a strong interaction between excess electron 
and lattice, such that the dimension of a polaron is comparable to a lattice 
constant. 

A different kind of approach, resulting in what is called the small polaron 
picture, has been worked out by Sewell (1958) and Holstein (1959). In 
the smali polaron picture it is taken into account that at T’> 0 the electron 
plus its associated lattice displacement polarizations correspond to an 
electron with a cloud of virtual and real phonons. The method of Sewell 
is essentially that of a Bloch tight-binding approximation, but taking the 
electron-lattice interaction into account when calculationg polaron eigen- 
states. Therefore, the overlap integrals which play a central role in the 
_ tight-binding approximation (Bloch 1928), depend on electron and phonon 
quantum numbers. By taking an appropriate thermal average over the 
eigenvalues, finally the low lying energy eigenvalues of the small polarons 
can be defined in terms of a temperature-dependent effective mass m**. 
It may be expected that, in a certain temperature interval around the 
Debye temperature which corresponds to the frequency of the polarization 
phonons, the small polaron effective mass is given by the simple formula 

m7, ** — 97* i ¢ 

m™ an, ex(z): «Odo See eo) 
Here m* is the rigid lattice mass and 7’, is a parameter represented by a 
rather complicated expression which is given in Sewell’s paper. 

From eqns. (1) and (2) it follows that the Seebeck coefficient of a polar 
crystal to which the small polaron picture applies is given by 

k 3a) 2(2am*k pT')3 2 
o= —F(b+ 57+ log |) Ce Ne etal 
with the above-mentioned restriction on the temperature interval. Thus, 
the Seebeck coefficient can depend linearly on temperature, provided the 
polaron concentration n is constant. 

Now, according to the results of experimental investigations at General 
Atomic (Kurnick, private communication) ; the following are characteristic 
for Ce-S semiconductors : 

1. Astrong, almost linear rise of the Seebeck coefficient with temperature 

between 300 and 1000°K. 

2. A slow, non-exponential decrease of the electrical conductivity with 

increasing temperature in samples with S: Ce ratios 1-3-1-4. 

3. An almost undetectable, but finite, Hall voltage. 

If one assumes that Ce-S semiconductors are characterized by a strong 
polar bonding between the lattice particles, a possible explanation of the 
measured temperature dependence of the Seebeck coefficient is given by 
eqn. (3). Then, the linear rise of « may be consistent with a slight tem- 
perature dependence of the electrical conductivity, provided n=const. 
and provided the polaron mobility—a delicate subject—is independent of 
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m** and is determined by the interaction of polarons with acoustical 
phonons. It turns out, from a sketch of the polaron mobility by Frohlich 
and Sewell (1959), that for the interaction of small polarons with acoustical 
phonons via two phonon scattering processes, the mobility might be 
independent of m** (and of the order of 1 cm?v—!sec~). 
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Beta-decay of '"N: Conservation of Spin and Parity in !°O+ 


By D. E. Ausuresr, R. E. Prxtry, D. H. Wirxinsont 
Brookhaven National Laboratory, Upton, New York, U.S.A. 


and P. Donovan 


Brookhaven National Laboratory and Bell Telephone Laboratories, 
Murray Hill, New Jersey, U.S.A. 


[Received October 15, 1960] 


Tue degree to which parity is conserved in the strong interactions is an 
open question. The early experiments (Tanner 1957, Wilkinson 1958 a, 
Segel e¢ al. 1958) made it unlikely that ¥?, the intensity of the component 
of irregular parity in a typical nuclear wave function, can be greater than 
about 10~*. Subsequent work (see, e.g. Blin-Stoyle 1960) has suggested 
F*<10-8. We report here an experiment designed to sharpen this limit. 

The ground state of ®N is 2-. It is known to decay with a half-life of 
74sec and a total energy release of 10-40 Mev through several states of 16O 
(Ajzenberg-Selove and Lauritsen 1959); that at 8-88 Mev is 2~ and enjoys 
a beta branch of 1-19, (Wilkinson et al. 1956). This state is energetically 
unstable by 1-73 Mev to decay into the ground state of !*C plus an alpha 
particle but such decay is rigorously forbidden if spin and parity are good 
quantum numbers. In the present work we seek this decay. 

We irradiated a thick target of TiN with a beam of deuterons of 1-1 Mev 
which impinged at a grazing angle of 30°. The consequent '*N decay was 
examined by a solid state counter close to the target and at 90° to the 
bombarding beam (which was interrupted during the counting period). 
Owing to the grazing incidence of the bombarding beam and the steeply 
energy dependent cross section for the reaction N(d, p)'®N the induced 
16N activity is rather superficial and the source is not thick to the possible 
alpha-particles from the decay of the 8-88 Mev state. 

The left-hand side of the figure shows a typical pulse-height spectrum 
from the solid-state counter. The exponential rise at low channel 
numbers is due to the great numbers of energetic beta-particles that stream 
through the counter. The solid line on the left-hand side of the figure is 
not drawn through the points but is taken directly from another run in 
which the solid-state counter was covered with a 1 mil. Al foil, sufficiently 
thick to stop alpha-particles from the decay of any state of 160 accessible 
to the beta-decay of 1*N but not thick enough appreciably to affect the 
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beta particles responsible for the small pulses. It is clear by the departure 
of the points from the line at high channel numbers that alpha-particles are 
indeed emitted following the beta decay of *N. That the alpha-particles 
are associated with 1®N-decay was demonstrated by a rough half-life 
measurement on the counts at high channel numbers : this gave 7:6 + 1-Osec. 

These alpha particles come from the decay of the broad 1~ state in 160 of 
nominal excitation 9-58Mev which is accessible to an allowed beta 
transition from 16N. The right-hand side of the figure shows the difference 
between the experimental points (grouped in threes) and the solid line of 
the left-hand side. The solid line of the right-hand side is the alpha- 
particle spectrum computed for the break-up of the 9-58 Mev level following 
its population through beta-decay. This computation allows for the 
considerable width of the level (Hill 1953) for the beta-decay phase-space 
factors and for the various effects of target thickness. The energy 
calibration of the solid state counter was effected with the aid of an 
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Decay of 6N as seen with a solid-state counter. The left-hand half of the 
figure shows as the experimental points the pulse distribution in the un- 
shielded counter and as the full line the distribution in the counter when 
shielded with a 1 mil. Al foil. The right-hand half of the figure shows 
the difference between the points and the line of the left-hand half ; the 
eee is the alpha-particle distribution computed for the decay of 
: e broad 1~ state in 1°O of nominal energy 9-58 Mev. Alpha particles 
rom the decay of a sharp level at 9-58 Mev in 160 would reach to channel 
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alpha-particle source and no adjustment in abscissa is made in comparing 
experimental and theoretical distributions. The 1N source strength was 
monitored using a 5in. x 4in. Nal (Tl) crystal and as a result we can quote 
the branching ratio 3-0 x 10-® for transitions to the broad 9-58 Mev state. 
This corresponds to log ft=7-0. (Using as beta-decay energy the value 
appropriate to the maximum of the computed and observed alpha- 
particle spectrum, viz. an effective level energy of about 9-45 Mev rather 
than 9-58 Mev; the lower energy is favoured by the greater energy released 
in the beta-decay.) The uncertainty in the branching ratio is about 50%. 
This log ft value is interestingly high for an allowed transition. We may 
recollect however that this 1~ state is not generated as belonging to 
1p~*(2s, Id) by the independent particle model in intermediate coupling 
(Elhott and Flowers 1957) which gives an excellent account of most of the 
other odd parity 7’=0 states below 12 Mev or so. It is presumably a state 
of three-quantum excitation and so the transition is theoretically forbidden 
although allowed by spin and parity; the large log ft value is therefore a 
confirmation of the expectation of the model. (Note that the jj-coupling 
model of A=16 (Unna and Talmi 1958) describes ®N as 1p,j.11d,/. and 
the 9-58 Mev state of 160 as 1p,/.12s,.: it therefore seeks to understand 
the large log ft value as due to /-forbiddenness.) 

The absence of alpha particles from the break-up of the 2+ state of 1*O 
at 9-84Mmev implies logft>7-3 for the first forbidden transition to this 
state. 

If the 2- state at 8-88 Mev broke up by alpha-particle emission it would 
affect the four points of lowest energy on the right hand side of the figure. 
The absence of any significant excess in this region above the solid line 
representing the alpha-particles from the 9-58 Mev state corresponds to 
l,/T, < 2-4 x 10~ for the 8-88 mev level (this figure represents one standard 
deviation within which the ratio is zero). 

We must now estimate the lifetime of the 8-88Mev state. The IPM 
calculation successfully predicts its position. For the preferred values of 
the parameters of the model it also predicts satisfactorily not only the 
branching ratios of the transitions from this state to the 7-12, 6-14Mev 
and ground states (experiment: 0:15:1:0-09; theory: 0-06: 1: 0-05), 
but also the M1: E2 mixing ratio in the case where it has been estimated 
namely the transition to the 6-14Mev level. It also predicts accurately 
the properties of the low-lying states of *N (Wilkinson 1957, Freeman 
and Hanna 1957) and the speed of the beta-decay to the 8-88 Mev level 
(experimental log ft=4:4; theoretical log ft=4-2): We may therefore 
have some confidence in its prediction of the absolute width of the 8-88 Mev 
level which is 2 x 10~%ev. ; 

An indication of the plausibility of this lifetime comes from the isotopic 
spin forbidden E1 transition between the 8-88 Mev and 6-91 Mev (27) states 
that is not treated by the IPM. If we apply the usual recipe (Wilkinson 
1956) to estimate the speed of an allowed KI transition of this energy 
(0-032 of a Weisskopf unit on a radius constant of 1:5 x 10-4 em) we find 
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0-18ev, We can now combine the theoretical IPM speed with the experi- 
mental branching ratios of the 8-88 Mev level to deduce a width for the E1 
transition of 1-1 x 10-*ev and so an effective isotopic spin impurity of 
ap2(1)¥6x10-. This in turn, using the recommended procedure, 
(Wilkinson 1958b) gives the value of the Coulomb matrix element of 
H° = a9(1) x Ag, = 0-025 x 3-9 0-10 Mev which is eminently reasonable. 

Combining our experimental limit on T,/I’, with this figure for I’, we 
conclude I, < 4:8 x 10-8ev. In order to give a value for 7* we must now 
fix on a typical value for [’, for a transition allowed by spin and parity. 
We are well placed to do this since in 160 the reduced alpha-particle widths 
are known for nine levels of 7’'=0 (Lane 1960). The logarithmic mean of 
these reduced widths is 0-24 single particle units (the arithmetic mean is 
0-64). If this value of 0-24 is used to predict a typical allowed value of I, 
for a state at 8-88 Mev (/= 2) we find 6-7 x 10° ev (using an interaction radius 
of 1:4 x (4,18 + A,"8) x 10-3cm). This then suggests 
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We may similarly determine a hmit for spin conservation if we assume 
that the spin-non-conserving interaction introduces a small element of 1- 
into the 2~ state. The allowed value of I’, is now (/=1) 2-0 x 10*ev and so 
the limit on the intensity of the admixed J =1 wave-function is 2 x 10-1. 

Although we must voice the usual cautions associated with such 
estimates as these we are probably safe in saying F?< 10-1 with the limit 
on spin conservation somewhat sharper. 

This experiment is susceptible of considerable refinements that are now 
being put in train. 
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REVIEWS OF BOOKS 


Quantum Theory of Atomic Structure, Vol. 1. By Joun C. Suater. (New 
York: McGraw-Hill, 1960.) [Pp. 502.] Price £4 5s. 6d. 


ANOTHER volume from Professor Slater’s fluent pen is always welcome. 
Indeed, we have the promise of more good things to come, for this is just the 
first of a series, adding up to a treatise on the grand scale, covering the applica- 
tion of quantum mechanics to the non-nuclear properties of matter. The first 
200 pages set the stage with a more or less conventional account of basic quan- 
tum theory; the remainder of this volume deals with the structure of atoms 
(except for spin-orbit and hyperfine interactions), up to the latest calculations by 
the M.I.T. team on the multiplet structure of the Iron group. The exposition is, 
as always, lucid in the extreme; students will be grateful for the clear and 
simple explanation of difficult topics, and the profound good sense of the 
comments on the physical interpretation of the mathematical results. In fact, 
for a work on such a scale, the treatment may be a little informal in places; 
we sometimes lack a precise statement of the general theorems and principles so 
that we may know exactly where we stand. Perhaps he has really tried to put 
in too much. There are other books, or lecture courses, where one would have 
already learned trivia such as wavepropagation in strings, or where one could 
try to master highbrow mathematical techniques such as the Calculus of 
Variations. Much better, for example, merely to refer to Kemble’s treatment 
of the principle of Least Action than to try to explain it in a two-page appendix. 
Even Professor Slater cannot make a subtle mathematical theory as simple as 
that. But perhaps this is only a comment on the education of the American 
graduate students at whom this work is chiefly aimed. J: Maze 


Nuclear Spectroscopy. Edited by F AszENBERG-SELOVE. (Academic Press, 
1960.) [Pp. 1147.] Two volumes. Price $16 each. 


Tus excellent two volume work is a readable and comprehensive introduction 
to experimental nuclear spectroscopy. It consists of a series of articles by some 
thirty authors, distinguished for this contributions to the subject and sets out 
to meet the needs both of graduate students entering the experimental field 
and of others working in related areas who wish to familiarize themselves with 
present ideas on nuclear spectroscopy. The first volume contains the experi- 
mental parts of the subject and is largely devoted to articles on the spectro- 
scopy of charged particles, gamma rays and neutrons. In these attention is 
focused on up-to-date instruments and techniques and the whole work is 
liberally illustrated with typical results obtained with them. The second 
volume contains material of more theoretical nature and is divided into a section 
on the analysis of experimental data and a section on nuclear models. The 
former collects together articles which describe the reduction of experimentally 
measured quantities to quantities which can be directly related to theory. The 
latter section contains the theory in the form of articles on the shell model the 
collective model and the complex potential model of the nucleus. _ 

The broad treatment of the whole subject and the general lucidity of pre- 
sentation, are likely to make this book a good starting point for graduate students 
for some years to come. It will also appeal to experimental physicists working 
in this field because of its stimulating and refreshing tone. It is well produced, 
referenced and indexed and can be strongly recommended to those interested in 


low energy nuclear physics. M1. A. G. 
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Progress in Nuclear Energy: Technology, Engineering and Safety, Vol. 2. 
P.N.E. series. Edited by R. Hurst, R. N. Lyon and C. M. NIcHoLts. 
(Pergamon Press, 1960.) [Pp. xvl+800.] Price £5 5s. 0d. 


Tuts book is one of a series based on the Proceedings of the Second International 
Conference on the Peaceful Uses of Atomic Energy, and it contains thirty-four 
papers by various authors on technological aspects of nuclear reactors, in 
particular on reactor chemistry, reactor materials, engineering experience and 
practice, engineering studies, and engineering aspects of reactor safety. As its 
title might suggest, it is something of a hotch-potch of miscellaneous themes and 
lacks the unity that is to be found in other volumes devoted to a single topic. 
Despite the editors’ efforts to deal with this basic problem, and despite their 
skilful ordering of the articles, it does not read as a single book but as a collection 
of separate reviews and research papers which cover such widely different topics 
as, for example, ion exchange properties of hydrous oxides, the design of irradia- 
tion experiments, and hydrodynamic problems of reactor containment. The 
chemical aspects of the behaviour of materials in reactors are dealt with in 
detail, although magnesium is a surprising omission. The physical aspects are 
dealt with rather scantily and there is nothing on Wigner effects in graphite 
despite the fact that the chemistry of the graphite —CO, system is allowed 70 
pages. There are several articles on liquid metal fuels which might have fitted 
better into the companion volumes on Metallurgy and Fuels. The most 
systematically developed themes in the volume are those dealing with methods of 
cooling and with the safety of reactors. The mass of technical data presented 
vividly illustrates the well-known view that a thousand man-years of scientific 
and engineering effort are required to develop a single new types of power- 
producing reactor. ASH: 


Advances in Applied Mechanics, Vol. IV. Edited by H. L. DrypeEn, Tu. v. 
KArmAN and G. Kurrti. (London and New York: Academic Press 1960.) 
[Pp. 294.] Price $9.00. 


THIS is the latest of a series of volumes containing review articles on recent 
work in continuum mechanics and allied fields. It maintains the generally 
high standard set by its predecessors, and all that need be said about it is to list 
the titles and authors of the five articles contained in it. They are : 

The theory of unsteady laminar boundary layers, by K. Stewartson; 


ips at theory with dissociation and ionization, by G. Ludwig and 
. Heil; 


The propagation of shock waves along ducts of varying cross section and 
by W. Chester; 


Similarity and equivalence in compressible flow, by K. Oswatitsch; 
Karman vortex streets, by R. Wille. Gik. 
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G. REMAUT and J. VENNIK Phil. Mag. Ser. 8, Vol. 6, Pl. I. 


Oscilloscope trace of response of sodium chloride to cyclic bending (100 ¢/s). 
One division corresponds to | my. 


Fig. 3 


‘Dual trace oscillogram comparing the signal from the specimen with that from 
the vibrating capacitor. 


G. REMAUT and J. VENNIK Phil. Mag. Ser. 8, Vol.6; Pi-=2: 


A photo of the sawtooth current pulses and the obtained electrical effect. The 
polarity of the second pulse is reversed to that of the first pulse. 


Fig. 10 


Photo of the signal obtaine: 


| by imposing on the specimen a series of pulses of the 
same polarity. 


R. J. STOKES et al. Phil. Mag. Ser. 8, Vol. 6, PI. 3. 
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Photomacrograph showing the distribution of slip after fracture in crystal of 
fig. 1 (a). Slip bands revealed by etching and oblique illumination. (x 5.) 


Fracture surface for crystal of fig. 2 illustrating the characteristic (110) [001] 
step or ‘line’ source. Note that a number of steps are produced 
each corresponding with the structure of the slip band on the left of 


fig. 2. (x 20.) 
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Composite photomicrograph showing the location of the (110) [001] slit fracture 
source of crystal (6), fig. 1, with respect to the slip bands. (x 125.) 


R. JMSTOKES@etral: Phil. Mag. Ser. 8, Vol. 6, PI. 4. 


(a) 
0% 
Elongation 


(d) 
06% 


Selected frames from a motion picture showing the development of slip and 
the formation of cleavage cracks at slip band intersections in magnesium 
oxide. 

The crack located between the slip bands A was responsible for fracture. 
Frames (g) and (i), (j) and (&) were consecutive, i.e., 1/16 sec apart. 
Polarized illumination. (x 5.) 
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RalLSTOKES et al Phil. Mag. Ser. 8, Vol. 6, PI. 5. 


Fig. 9 
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Fig. 6 
Etch pits in an expanding band of single slip in magnesium oxide. 
The direction of growth is from right to left. Note the extremely sharp 
demarcation between slipped and unslipped material. (x 12.) 


Fig. 9 
Two ‘ blocks ’ of slip in a sprinkled magnesium oxide crystal deformed in tension. 


(a) Polarized transmitted illumination. P-plane of polarizer, A-plane of 
analyser. (x 12.) 

(6) Normal transmitted illumination. Etched surface. (x 12.) 

(c) Normal reflected illumination. Enlargement of block boundary in (b) to 


show position of cracks. (x 25.) 
Arrows in (6) and (c) indicate directions for dense slip in each block. Check 


by looking along photograph (c). 

Note: The stress birefringence and etch pit patterns are most marked at 
block boundary and for isolated slip bands penetrating across block 
boundary not the regions of dense single slip. 


R. J. STOKES tet al. Phil. Mag. Ser. 8, Vol. 6, PI. 6. 


Dense intersecting slip in the tension region of a sprinkled crystal deformed 
by bending. The bending axis is perpendicular to the plane of the 
paper. (x 19.) 


Fig. 12 


Tension surface of a sprinkled crystal fractured by bending. Tension direction 
horizontal. The two halves have been butted back together for the 
picture. Note the shattering along the edges of the specimen and the 
longitudinal cracks. (x 4.) 


Fig. 13 
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Anticlastic kinking along the edge of the tension surface in crystal of fig. 12. 
Tension direction horizontal. Fine slip on {110} planes making 45° 


with cube edge can be seen within kinked region only. Oblique illumina- 
tion. (x 56.) 


Reads STOKES etal. Phil. Mag. Ser. 8, Vol. 6, Pl. 7. 


Longitudinal crack in tension surface of a sprinkled magnesium oxide crystal 
fractured by bending. Tension direction vertical. (x 250.) 


Fig. 15 


surface of crystal in fig. 14. Fracture source corresponds with the 


Fracture 
step indicated by arrow in fig. 14. (x 20.) 


R. ILLINGWORTH Phil. Mag. Ser. 8, Vol. 6, PI. 8. 
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Luminescent fragment of KBr: Sn. Mag. x13. 


“Aes 
Fig. 2 


Same field of view as fig. 1, but with the ultra-violet light from the mercury 
lamp filtered out. Mag. x 13. 


Fig. 3 


Transmission photograph showing the surface structure. Mag. x13. 


S. HOWE and C. ELBAUM Phil. Mag. Ser. 8, Vol. 6, PI. 9. 


<100> crystal deformed at 400°c. Arrow indicates region where deformation 
took place on different {111} planes in different regions of the crystal 
(no necking). Approximately x 21. 


Fig. 5 


<110> crystal deformed at 400°c. Approximately x 21. 


Fig. 6 
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<111> crystal deformed at 400°c. Approximately x 23. 


Fig. 7 


<100) erystal deformed at 600°c. Note rounded edges of slip bands. Arrows 

show region where deformation occurred on different {111} planes. 
Approximately x 23. 

Fig. 8 


B®: 


<110> crystal deformed at 600°c. Note sharp-edged slip bands. Approxi- 
mately x 24. 


Fig. 9 


é111> crystal deformed at 600'c. Approximately x 25. 


j-G) DUTHIEvetal: Phil. Mag. Ser. 8, Vol. 6, PI. 10. 


Contact prints of individual emulsion sheets, showing cascades visible to the 
unaided eye. The cascades shown in (A) and (B) are members of a 
‘family ’ of about 30 cascades, distributed through the Comet 1 stack : 
the individual cascade axes are parallel to each other, and most were 
produced by y-rays resulting from a nuclear interaction of primary 
energy ~ 10 Bev occurring in the air about 1km above the aircraft. 
In (C), 7 cascades (all unrelated) can be distinguished in a single 
emulsion. Their positions are indicated by the intersections of the arrows. 


J. G. DUTHIE et al. Phil. Mag. Ser. 8, Vol. 6, PI. I]. 


Fig. 9 
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Photomicrographs of event 192a. A single energetic neutral pion (of energy 
~ 3800 Bev) was generated in a nuclear interaction in the dural end- 
plates clamping the stack. Both y-rays from the 7°-decay convert 
after traversing ~7 cm of emulsion, and the tracks of the two electron 
pairs are shown in (a), with a lateral separation in the emulsion plane 
of ~1-8y. In(b), 11.1. beyond (a), the original pairs are now accompanied 
by other electrons. (c) shows the appearance of the cores 4r.1. beyond (a), 
where their separation is ~2-7 yu. In (d), 67.1. from (a), the two cores 
are no longer distinguishable. 


